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Abstract 



Quantum error-correcting codes will be the ultimate enabler of a future quantum computing or 
quantum communication device. This theory forms the cornerstone of practical quantum informa- 
tion theory. We provide several contributions to the theory of quantum error correction — mainly 
to the theory of "entanglement-assisted" quantum error correction where the sender and receiver 
share entanglement in the form of entangled bits (ebits) before quantum communication begins. 
Our first contribution is an algorithm for encoding and decoding an entanglement-assisted quan- 
tum block code. Wc then give several formulas that determine the optimal number of ebits for an 
entanglement-assisted code. The major contribution of this thesis is the development of the the- 
ory of entanglement-assisted quantum convolutional coding. A convolutional code is one that has 
memory and acts on an incoming stream of qubits. We explicitly show how to encode and decode a 
stream of information qubits with the help of ancilla qubits and ebits. Our entanglement-assisted 
convolutional codes include those with a Calderbank-Shor-Steane structure and those with a more 
general structure. We then formulate convolutional protocols that correct errors in noisy entan- 
glement. Our final contribution is a unification of the theory of quantum error correction — these 
unified convolutional codes exploit all of the known resources for quantum redundancy. 



xi 



Chapter 1 



Introduction 



// computers that you build are quantum, 
Then spies of all factions will want 'em. 

Our codes will all fail, 
And they'll read our email, 
Till we've crypto that's quantum, and daunt 'em. 

— Jennifer and Peter Shor 



QUANTUM computation and communication have enormous potential to revolutionize the 
way we compute and communicate. From communicating by teleportation of quantum 
information [3] to breaking RSA encryption [4, 5], quantum technological breakthroughs will 
change society in a way that is difficult to predict. 

One might say that the field of quantum computation began when Richard Feynman and others 
suggested that a computer operating according to the principles of quantum mechanics would have 
advantages over a computer operating according to the laws of classical physics [6, 7, 8]. They 
suggested that it might be able to simulate quantum-mechanical processes more efficiently than a 
classical computer could. Later work then showed that it was possible to achieve this simulation 
speedup [9]. Other major advances in the theory of quantum computing followed — examples 
are Shor's algorithm for breaking RSA encryption [4, 5] and Grover's algorithm for searching a 
database [10, 11]. 
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The field of modern quantum communication began with the discovery of quantum key distri- 
bution [12]. Major discoveries for quantum communication followed examples are the ability to 
send a quantum bit using two classical bits and a bit of entanglement (quantum teleportation [3]) 
and the ability to send two classical bits by sending a quantum bit and consuming a bit of entan- 
glement (quantum superdense coding [13]). Quantum information theorists then began to think 
more deeply about the ways in which we could combine the resources of classical communication, 
quantum communication, and entanglement to formulate new communication protocols. 

In spite of these spectacular advances for the theories of quantum computation and communi- 
cation, a dark shadow lingered over them. A few authors disputed that reliable quantum compu- 
tation or communication would be possible because small quantum errors would accumulate as the 
computation proceeds or channel errors would destroy the quantum bits as we communicate them 
[14]. Rolf Landauer even urged his colleagues to include the following disclaimer in their papers 
on quantum computation: "This scheme, like all other schemes for quantum computation, relies 
on speculative technology, does not in its current form take into account all possible sources of 
noise, unreliability and manufacturing error, and probably will not work." Various other obstacles 
such as the no-cloning theorem [15] and measurement destroying a quantum state seemed to pose 
an insurmountable barrier to a protocol for quantum error correction. 

Despite the aforementioned obstacles, Shor demonstrated the first quantum error-correcting 
code that reduces the negative effects of decoherence on a quantum bit [16]. Shor's code overcame 
all of the above difficulties and established the basic principles for constructing a general theory 
of quantum error correction [17, 18, 19]. Shor's code exploits many of the signature principles of 
quantum mechanics: superposition, entanglement, unitary evolution, and measurement. Mermin 
proclaims it a "miracle" that quantum error correction is even possible [20]. 

1.1 Quantum Error Correction 

Quantum error correction theory [16, 21, 22, 17, 19] now plays a prominent role in the practical 
realization and engineering of quantum computing and communication devices. The first quantum 
error-correcting codes [16, 21, 22, 19] are strikingly similar to classical block codes [23] in their 
operation and performance. Quantum error-correcting codes restore a noisy, decohered quantum 
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state to a pure quantum state. Any future quantum information processing device will operate 
faithfully only if it employs an error correction scheme. This scheme can be an active scheme [17], 
a passive scheme [24, 25, 26], or a combination of both techniques [27, 28, 29, 30, 31]. 

The field of quantum error correction has rapidly expanded in recent years because there 
are so many ways in which one can correct quantum error. We briefly introduce three notable 
developments in the theory — these three are by no means exhaustive. The developments include 
the stabilizer formalism, the entanglement-assisted stabilizer formalism, and the convolutional 
stabilizer formalism. 

1.1.1 Stabilizer Quantum Coding 

Gottesman formalized the theory of quantum block coding by establishing the stabilizer formalism 
[17]. A stabilizer quantum error-correcting code appends ancilla qubits to information qubits that 
we want to protect. A unitary encoding circuit rotates the Hilbert space of the information 
qubits into a subspace of a larger Hilbert space. This highly entangled, encoded state corrects for 
local noisy errors. A quantum error-correcting code makes quantum computation and quantum 
communication practical by providing a way for a sender and receiver to simulate a noiseless qubit 
channel given a noisy qubit channel that has a particular error model. 

The stabilizer theory of quantum error correction allows one to import some classical binary or 
quaternary codes for use as a quantum code [19]. The Calderbank-Shor-Steane (CSS) construction 
is the name for the method for importing classical binary codes [32]. This idea of importing codes 
is useful because quantum code designers can utilize classical block codes with high performance 
to construct quantum codes with high performance. The only "catch" when importing is that the 
classical code must satisfy the dual-containing or self-orthogonality constraint. Researchers have 
found many examples of classical codes satisfying this constraint [19] , but most classical codes do 
not. 

1.1.2 Entanglement-Assisted Stabilizer Quantum Coding 

Bowen was the first to extend the stabilizer formalism by providing an example of a code that 
exploits entanglement shared between a sender and a receiver [33]. The underlying assumption 
of Bowen's code is that the sender and receiver share a set of noiseless ebits (entangled qubits) 
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before quantum communication begins. Many quantum protocols such as teleportation [3] and 
superdense coding [13] are "entanglement-assisted" protocols because they assume that noiseless 
ebits arc available. 

Brun, Devctak, and Hsich extended the standard stabilizer theory of quantum error correction 
by developing the entanglement-assisted stabilizer formalism [34, 35] . They included entanglement 
as a resource that a sender and receiver can exploit for a quantum error-correcting code. They 
provided a "direct-coding" construction in which a sender and receiver can use ancilla qubits and 
ebits in a quantum code. An ebit is a nonlocal bipartite Bell state 

|$+) = (|00) + |ll))/^/2. 

Gottesman later showed that their construction is optimal [36] — it gives the minimum number of 
ebits required for the entanglement-assisted quantum code. 

The major benefit of the entanglement-assisted stabilizer formalism is that we can construct 
an entanglement-assisted quantum code from two arbitrary classical binary block codes or from an 
arbitrary classical quaternary block code. The rate and error-correcting properties of the classical 
codes translate to the resulting quantum codes. The entanglement-assisted stabilizer formalism 
may be able to reduce the problem of finding high-performance quantum codes approaching the 
quantum capacity [37, 38, 39, 40, 41] to the problem of finding good classical linear codes ap- 
proaching the classical capacity [42]. The entanglement-assisted stabilizer formalism thus is a 
significant and powerful extension of the stabilizer formalism. 

1.1.3 Convolutional Stabilizer Quantum Coding 

Another extension of the theory of quantum error correction protects a potentially-infinite stream 
of quantum information against the corruption induced by a noisy quantum communication chan- 
nel [43, 44, 45, 46, 47, 48, 1, 49, 50, 51, 52, 53]. Quantum convolutional codes have numerous 
benefits. The periodic structure of the encoding and decoding circuits for a quantum convolutional 
code ensures a low complexity for encoding and decoding while also providing higher performance 
than a block code with equivalent encoding complexity [1]. The encoding and decoding circuits 
have the property that the sender Alice and the receiver Bob can respectively send and receive 
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qubits in an "online" fashion. Alice can encode an arbitrary number of information qubits without 
worrying beforehand how many she may want to send over the quantum communication channel. 

We believe that quantum convolutional coding theory is one step along the way to finding 
quantum error-correcting codes that approach the capacity of a noisy quantum communication 
channel for sending quantum information [37, 38, 39, 54, 55, 56, 57]. Poulin et al. have recently 
incorporated some of this well-developed theory of quantum convolutional coding into a theory 
of quantum serial-turbo coding [58] with the goal of designing quantum codes that come close to 
achieving capacity. 

The development of quantum convolutional codes has been brief but successful. Chau was the 
first to construct some quantum convolutional codes [59, 60], though some authors [1] argue that 
his construction does not give a true quantum convolutional code. Several authors have established 
a working theory of quantum convolutional coding based on the stabilizer formalism and classical 
self-orthogonal codes over the finite field GF {4) [43, 44, 48, 1]. Others have also provided a 
practical way for realizing "online" encoding and decoding circuits for quantum convolutional 
codes [43, 44, 46, 45]. 

Forney et al. have determined a method for importing an arbitrary classical self-orthogonal 
quaternary code for use as a quantum convolutional code [48, 1]. The technique is similar to that 
for importing a classical block code as a quantum block code [19]. Forney et al. showed how to 
produce quantum convolutional codes from classical convolutional codes by extending the ideas 
from the CSS construction to the convolutional setting [1]; but again, these imported classical 
convolutional codes have to satisfy the restrictive dual-containing constraint in order to form valid 
quantum codes. The dual-containing constraint is actually quite a bit more restrictive in the 
convolutional case because each generator for the quantum convolutional code must commute not 
only with the other generators, but it must commute also with any arbitrary shift of itself and any 
arbitrary shift of the other generators. Forney et al. performed specialized searches to determine 
classical quaternary codes that satisfy the restrictive self-orthogonality constraint [1] . 
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1.2 Connections to Quantum Shannon Theory 

Quantum error correction theory gives practical ways to build codes for protection of quantum 
information against decoherence. A quantum code typically encodes some number k of information 
qubits into a larger number n of physical qubits. We say that the rate of quantum communication 
for this code is k/n. 

One may ask the question: what is the maximum rate of quantum communication for a given 
noisy quantum communication channel? This question is the fundamental question of quantum 
Shannon theory and the answer to it is a quantity known as the capacity of the quantum channel 
for quantum communication, or more simply, the quantum capacity [37, 38, 39, 61]. Quantum 
Shannon theory has many other questions. What is the rate of classical communication for a 
noisy quantum communication channel? What is the rate if the sender and receiver have access 
to entanglement? What is the rate of simultaneous classical and quantum communication? The 
list of questions goes on and on. 

The goal of quantum Shannon theory is to quantify the amount of quantum communica- 
tion, classical communication, and entanglement required for various information processing tasks 
[37, 38, 39, 61]. Quantum teleportation and superdense coding provided the initial impetus for 
quantum Shannon theory because these protocols demonstrate that we can combine entanglement, 
noiseless quantum communication, and noiseless classical communication to transmit quantum or 
classical information. In practice, the above resources are not noiseless because quantum systems 
decohere by interacting with their surrounding environment. Quantum Shannon theory is a col- 
lection of capacity theorems that determine the ultimate limits for noisy quantum communication 
channels. Essentially all quantum protocols have been unified as special cases of a handful of 
abstract protocols [61]. 

The techniques in quantum Shannon theory determine the asymptotic limits for communi- 
cation, but these techniques do not produce practical ways of realizing these limits. This same 
practical problem exists with classical Shannon theory because random coding techniques deter- 
mine the limit of a noisy classical communication channel [62] . 
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1.2.1 Communication Protocols 

The fundamental theorem of quantum Shannon theory is the celebrated Lloyd-Shor-Devetak quan- 
tum capacity theorem. It determines the capacity of a given noisy quantum communication chan- 
nel for reliable quantum communication. Suppose that a noisy quantum channel J\f connects a 
sender Alice ^4 to a receiver Bob B. Let [q — > q\ denote one qubit of noiseless quantum communi- 
cation. We can state the coding theorem as a resource inequality: 

{M)>Q[q^q], (1.1) 

The above resource inequality states that n uses of the noisy quantum channel J\f are sufficient 
to communicate nQ noiseless qubits in the limit of a large number n of channel uses. The rate Q 
is equal to a quantity known as the coherent information: I{A)B) [32]. The entropic quantity is 
maximized with respect to any resource state |'^)'*^^ associated to the noisy channel M and shared 
between Alice, Bob, and the environment E. The goal of the stabilizer codes and convolutional 
stabilizer codes mentioned respectively in Sections 1.1.1 and 1.1.3 is for their rates to approach 
the maximum capacity given above. 

Another example of an important capacity theorem from quantum Shannon theory results from 
the "father" protocol [61]. The father capacity theorem determines the optimal trade-off between 
the rate ii^ of cbits (entangled qubits in the state |(f>+)'^'^ = (|00)'^^ + |ll)^^)/\/2) and the rate Q 
of qubits in entanglement-assisted quantum communication. This protocol is the "father" protocol 
because it generates many of the protocols in the family tree of quantum information theory [61]. 
The nickname "father" is a useful shorthand for classifying the protocol — there exists a mother, 
grandmother, and grandfather protocol as well. Let [gg'] denote one ebit of entanglement. The 
following resource inequality is a statement of the capability of the father protocol: 

{M) + E[qq]>Q[q^q]. (1.2) 

The above resource inequality states that n uses of the noisy quantum channel M and nE noiseless 
ebits are sufficient to communicate nQ noiseless qubits in the limit of large n. The rates E and 
Q are respectively equal to \l {A; E) and \l {A; B). The entropic quantities are maximized with 
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respect to any resource state associated to the noisy channel M and shared between 

Ahce, Bob, and the environment E (the rate E is different from the environment E). The father 
capacity theorem gives the optimal limits on the resources, but it does not provide a useful quantum 
coding technique for approaching the above limits. The goal of the entanglement-assisted stabilizer 
formalism mentioned in Section 1.1.2 is to develop quantum codes that approach the maximal rates 
given in the father capacity theorem. We spend a significant portion of this thesis (Chapters 5, 
6, and 7) developing entanglement-assisted convolutional stabilizer codes. Future research might 
be able to incorporate these convolutional codes into a framework for designing codes that come 
close to achieving the maximal rates. 

Another important capacity theorem determines the ability of a noisy quantum channel to send 
"classical-quantum" states [63]. Let [c — > c] denote one bit of noiseless classical communication. 
The result of the classical-quantum capacity theorem is also a resource inequality: 

{M)>Q[q^q]+R[c^c]. (1.3) 

The resource inequality states that n uses of the noisy quantum channel M are sufficient to 
communicate nQ noiseless qubits and nR noiseless classical bits in the limit of large n. The 
rates Q and R are respectively equal to I (X; B) and / {A)BX). The entropic quantities are with 
respect to any state resulting from sending the A' system of the following classical-quantum state 

la;) (a;|^0 (<A,|^^' (1.4) 

X 

through the quantum channel M^^^ . X and A are systems that the sender keeps. This theorem 
proves that we can devise clever classical-quantum codes that perform better than time-sharing a 

noisy quantum channel A/" between purely quantum codes and purely classical codes. 

The "grandfather" capacity theorem determines the optimal triple trade-off between qubits, 
ebits, and classical bits for simultaneous transmission of classical and quantum information using 
an entanglement-assisted noisy quantum channel M [64]. The grandfather resource inequality is 
as follows: 

{N) + E[qq]>Q[q^q]+R[c^c]. (1.5) 
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The above resource inequality is again an asymptotic statement and its meaning is similar to that in 
(1.1), (1.2), and (1.3). The optimal rates in the above resource inequality coincide with the father 
inequality (1.2) when R = 0, with the classical-cjuantum inequality (1.3) when E = 0, and with 
the quantum capacity (1.1) when both R. — and E — 0. The optimal strategy for the grandfather 
protocol is not time-sharing the channel between father codes and entanglement-assisted classical 
codes. It remains to be proven whether this optimal strategy outperforms time-sharing [64]. We 
develop grandfather convolutional codes in Chapter 8. Again, future research might be able to 
incorporate these grandfather convolutional codes into a framework for designing codes that come 
close to achieving the maximal limits given in the grandfather capacity theorem. 

1.2.2 Distillation Protocols 

The aforementioned protocols employ a noisy quantum communication channel in their operation. 

We now ask the question: given noisy entanglement, is it possible to produce noiseless entangle- 
ment using classical communication? Indeed, it is possible, and this procedure is entanglement 
distillation [65, 66]. The following resource inequality summarizes the resources consumed and 
produced in an entanglement distillation protocol: 

{pAB) + R[c^c]>Q[qq]. 

The statement of the resource inequality is as follows: given a large number n of noisy bipartite 
states p^^ shared between two parties, it is possible to generate nQ noiseless ebits by consuming 
nR classical bits in the limit of a large number n of copies of the noisy state p"^^ . The optimal 
rates are R — I{A;E) and Q = I{A)B) where the entropic quantities are with respect to a 
state IV')'^^^ that purifies the state p^^ . In Chapter 9, we develop convolutional entanglement 
distillation protocols whose goal is to approach the optimal rates given in the above resource 
inequality. 
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1.3 Organization of this Thesis 

This thesis represents a significant contribution to quantum error correction theory. In it, we 
give several methods for designing quantum codes. These methods should be useful in the future 
when coherent quantum encoding devices are available. The quantum code designer will have a 
plethora of techniques to exploit for protecting quantum information and the techniques in this 
thesis should be part of the arsenal. 

A quick glance over the thesis reveals that much of the mathematics involves only linear algebra 
over matrices of binary numbers or binary polynomials. It is a testament to the strength of the 
theory that it reduces to this simple form. We show in detail especially in Chapters 2 and 4 how 
to reduce the theory to have this simple form. This thesis has contributions from the following 
publications [52, 67, 68, 69, 53, 51, 70]. 

We structure the thesis as follows. The next chapter introduces the stabilizer formalism for 
quantum error correction. We show in this chapter how to represent quantum codes with bi- 
nary matrices. Chapter 3 introduces the entanglement-assisted stabilizer formalism and shows 
how to enhance the stabilizer formalism by assuming that the sender and receiver share entan- 
glement. This chapter gives two contributions: a method for determining the encoding circuit for 
a quantum block code and a method to determine the minimal amount of entanglement that an 
entanglement-assisted code requires. Chapter 4 reviews the convolutional stabilizer formalism. In 
particular, it details the operation of a quantum convolutional code and shows how to simplify 
the mathematics with matrices whose entries are binary polynomials. Chapters 5, 6, and 7 form 
the heart of this thesis. They show how to construct entanglement-assisted quantum codes that 
have a convolutional structure. The techniques developed represent a significant extension of the 
entanglement-assisted stabilizer formaUsm. Chapter 8 presents a unifying structure for quantum 
error correction under which most current schemes for quantum error correction fall. We present a 
method for distilling entanglement in a convolutional manner in Chapter 9. The major benefit of 
all of the above convolutional methods is that we can import arbitrary high-performance classical 
codes to construct high-performance quantum codes. 
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Stabilizer Quantum Codes 



Chapter 2 



"Fight entanglement with entanglement, " you say, 
'Cause entanglement keeps Eve away, 
The hell that she raises, 
She damps and dephases. 
She's no match — like UCLA! 

The stabilizer formalisni is a mathematical framework for quantum error correction [71, 17]. This 
framework has many similarities with classical coding theory, and it is even possible to import a 
classical code for use in quantum error correction by employing the CSS construction [21, 22, 32]. 
We review the stabilizer theory for quantum block codes. 

The stabilizer formalism is versatile because it appears not only in quantum error correction 
but also in other subjects of quantum information theory such as stabilizer or graph states [72], 
cluster-state quantum computation [73], and entanglement theory [74]. The stabilizer formalism 
simplifies the theory of quantum error correction by reducing the mathematics to linear algebra 
over binary matrices. 
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2.1 Review of the Stabilizer Formalism 



The stabilizer formalism exploits elements of the Pauli group 11. The set 11 = {7, X, Y, Z} consists 
of the Pauli operators: 





1 










1 




-i 
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1 = 






, x = 










, z = 









1 




1 







i 




-1 



The above operators act on a single qubit — a vector in a two-dimensional Hilbert space and 
axe the most important in formulating a quantum error-correcting code. Two crucial properties 
of these matrices axe useful: each matrix in 11 has eigenvalues equal to +1 or —1, and any two 
matrices in 11 either commute or anticommute. 

In general, a quantum error-correcting code uses n physical qubits to protect a smaller set of 
information qubits against decoherence or quantum noise. An n-qubit quantum error-correcting 
code employs elements of the Pauli group II". The set 11" consists of n-fold tensor products of 
Pauli operators: 

n" = {e''^^i(g)---(g)A„:Vie{l,...,n}, Aj gU, ^ e {0,7r/2,7r,37r/2}} . (2.1) 

Matrices in 11" act on a 2"-dimensional complex vector, or equivalently, an n-qubit quantum 
register. We omit tensor product symbols in what follows so that Ai • • • A„ = Ai (g) • • • (g) A„. 
The above two crucial properties for the single-qubit Pauli group 11 still hold for the Pauli group 
n" (up to an irrelevant phase for the eigenvalue property). The n-fold Pauli group 11" plays 
an important role in both the encoding circuit and the error-correction procedure of a quantum 
stabilizer code over n qubits. 

We can phrase the theory of quantum error correction in purely mathematical terms using 
elements of 11". Consider a matrix gi G 11" that is not equal to ±/. Matrix gi then has two 
eigenspaces each of dimension 2"^^. We can identify one eigenspace with the eigenvalue +1 and 
the other eigenspace with eigenvalue —1. Consider a matrix 92 G H" different from both gi and 
the identity that commutes with gi. Matrix 92 also has two eigenspaces each of size 2"^^ and 
identified similarly by its eigenvalues ±1. Then gi and (72 have simultaneous eigenspaces because 
they commute. These matrices together have four different eigenspaces, each of size 2"~^ and 
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Figure 2.1: The operation of a stabilizer code. Thin Hnes denote quantum information and thick 
lines denote classical information. Slanted bars denote multiple qubits. A sender encodes a multi- 
qubit state jV') with the help of some ancilla qubits |0). She sends the encoded state over a noisy 
quantum channel. The receiver performs multi-qubit measurements to extract information about 
the errors. He finally performs a recovery operation R to reverse the channel error. 



identified by the eigenvalues ±1,±1 of gi and 32 respectively. We can continue this process of 
adding more commuting and independent matrices to a set S. The matrices in S are independent 
in the sense that no matrix in 5 is a product of two or more other matrices in S. Adding more 
matrices from 11" to S continues to divide the eigenspaces of matrices in S. In general, suppose S 
consists of 71 — fc independent and commuting matrices gi, . . . , gn-k G 11"- These n — k matrices 
then have 2"^'^' different eigenspaces each of size 2*^ and identified by the eigenvalues ±1, . . . , 
±1 of gi, gn^k respectively. Consider that the Hilbert space of k qubits has size 2*^. A 
dimension count immediately suggests that we can encode k qubits into one of the eigenspaces of 
S. We typically encode these k qubits into the simultaneous +l-eigenspace of gi, . . . , gn-k- This 
eigenspace is the codespace. The operators gi, . . . , gn-k function in the same way as a parity check 
matrix does for a classical linear block code. An [n, k] quantum error-correcting code encodes k 
information qubits into the simultaneous +l-eigenspace of n — fc matrices gi, . . . , gn-k G II"- Its 
stabilizer S is an abclian subgroup of the rt-fold Pauli group 11" : S C 11" . Note that it is possible 
to multiply the generators in the stabilizer together and obtain an equivalent representation of 
the stabilizer in much the same way that we can add vectors in a basis together and obtain an 
equivalent basis. 
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The operation of an [n,k] quantum error-correcting code consists of four steps. Figure 2.1 
highlights these steps. 

(i) . A unitary operation U encodes both k information qubits in a general state and n — k 

ancilla qubits in the state jo)®""'' into the simultaneous +l-eigenspa€e of the matrices gi, 

• • • 5 9n—k' 

(ii) . The sender transmits the n encoded qubits by using the noisy quantum communication 

channel n times. 

(iii) . The receiver performs quantum measurements of the n — k matrices g\, . . . , Qn-k in the 

stabilizer S. These measurements learn only about errors that may occur and do not disturb 
the encoded quantum information. Each measurement gives a bit result equal to +1 or — 1, 
and the result of all the measurements is to project the n-qubit quantum register into one 
of the 2'""^^ different eigenspaccs of gi, .... gn-k- Suppose that no error occurs. Then 
the measurements project the n qubits into the simultaneous +l-eigenspace and return a 
bit vector consisting of n — fc ones. This "projection of the error concept" is one of the 
fundamental notions in quantum error correction theory. It suffices to correct a discrete 
error set with support in the Pauli group [16]. Now suppose that a quantum error in an 
error set f C 11" occurs. The error takes the encoded quantum state out of the codespace 
and into one of the other 2""*^ — 1 orthogonal eigenspaces. The measurements can detect 
that an error has occurred because the result of the measurements is a bit vector differing 
from the all ones vector. The receiver can identify uniquely which error in £ has occurred if 
the set £ satisfies the following quantum error correction conditions: 

VK, Ei&£ 3gi€S: {g,, ElEf,} = or eIe^ e S. (2.2) 

The first condition corresponds to the active error-correcting capability of the code, and the 
second condition corresponds to its passive error-correcting capability. 

(iv) . If the receiver can identify which error occurs, he can then apply a unitary operation R 

that is the inverse of the error. He finally performs a decoding unitary that decodes the k 
information qubits. 
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2.2 Clifford Encoding Unitary 

We comment briefly on the encoding operation U. The encoding operation f/ is a special type of 
unitary matrix called a Clifford operation. A Clifford operation U is one that preserves elements 
of the PauU group under conjugation: A e 11" UAW e 11". The CNOT gate, the Hadamard 
gate H, and the phase gate P suffice to implement any unitary matrix in the Clifford group [17]. 
A quantum code with the CSS structure needs only the CNOT and Hadamard gates for encoding 
and decoding. The matrix for the CNOT gate acting on two qubits is 



CNOT 



10 
10 
1 
10 



(2.3) 



the matrix for the Hadamard gate H acting on a single qubit is 



H = 



^/2 



1 1 
1 -1 



(2.4) 



and the matrix for the phase gate P acting on a single qubit is 



1 
i 



(2.5) 



The standard basis for H^ is X and Z because any element of H^ is a product of elements in 
this generating set up to an irrelevant phase. The standard basis for elements of the two-qubit 
Pauli group H^ is as follows 



Z I , I Z , XI, IX, 



(2.6) 



for the same reasons. 
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The Hadamard gate H transforms the standard basis of 11^ under conjugation as foUows: 

Z ^X, X^Z (2.7) 

and the phase gate P transforms the standard basis as follows: 

Z^Z, X^Y. (2.8) 

For the CNOT gate, the first qubit is the "control" qubit and the second qubit is the "target" 
qubit. The CNOT gate transforms the standard basis of 11^ under conjugation as follows 

Z I ^ Z I , I Z ^ Z Z , X I ^ X X , I X ^ I X , (2.9) 

where the first qubit is the control and the second qubit is the target. 

Section 3.5 of Chapter 3 details an algorithm that determines a Clifford encoding circuit for 
any stabilizer code or any entanglement-assisted stabilizer code (we review entanglement-assisted 
codes in the next chapter). 

2.3 The Logical Operators and the Information Qubits 

Another aspect of the theory of quantum error correction is later useful for our purposes in 
quantum convolutional coding. This aspect concerns the information qubits and the operators 
that change them. Consider that the initial uncncodcd state of a quantum error-correcting code 
is a simultaneous -|-l-eigenstatc of the matrices Zk+i, . . . , Zn where Zi has a Z matrix operating 
on qubit i and the identity / on all other qubits. Therefore, the matrices Zk+i, . . . , Zn constitute 
a stabilizer for the unencoded state. The initial unencoded logical operators for the k information 
qubits are Zi,Xi, . . . , Zk,Xh. The encoding operation U rotates the unencoded stabilizer matrices 
Zk+i, . . . , Zn and the unencoded logical operators Zi,Xi,. . . ,Zk,Xi. to the encoded stabilizer 
operators Z^+i, ■ ■ ■ , Zn and the encoded logical operators Zi,Xi, . . . , Zk,Xk respectively. The 
encoded matrices Zk+i, . . . ,Zn are respectively equivalent to the matrices gi, . . . , Qn-k in the 
discussion of the previous section. 
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The encoded operators obey the same conimutation relations as their unencoded counterparts. 
We would violate the uncertainty principle if this invariance did not hold. Therefore, each of the 
encoded logical operators commutes with elements of the stabilizer S. Let A denote an arbitrary 
logical operator from the above set and let Zi denote an arbitrary element of the stabilizer S. The 
operator AZi (or equivalently ZiA) is an equivalent logical operator because AZi and A have the 
same effect on an encoded state |'^): 

ZiAl^p) = AZi\tP) = A\^) . (2.10) 

We make extensive use of the above fact in OTir work. 

The logical operators also provide a useful way to characterize the information qubits. Gottes- 
man showed that the logical operators for the information qubits provide a straightforward way 
to characterize the information qubits as they progress through a quantum circuit [17]. As an 
example of this technique, he develops quantum teleportation in the stabilizer formalism. The 
logical operators at the beginning of the protocol are Xi and Zi and become X3 and Zs at the end 
of the protocol. This transformation implies that the quantum information in qubit one teleports 
to qubit three because the logical operators act on only qubit three at the end of the protocol. We 
use the same idea throughout this thesis to determine if our decoding circuits have truly decoded 
the information qubits. 

2.4 The Pauli-to-Binary Isomorphism 

A simple but useful mapping exists between elements of 11 and the binary vector space (Z2) . This 
mapping gives a simplification of quantum error correction theory. It represents quantum codes 
with binary vectors and binary operations rather than with Pauli operators and matrix operations 

respectively. We first give the mapping for the one-qubit case. Suppose [A] is a set of equivalence 
classes of an operator A that have the same phase: 

[A] = {PA\ (3eC, m = l}. (2.11) 
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Let [n] be the set of phase- free PauH operators where [11] = {[A] \ A & 11}. Define the map 
N : (Zz)^ ^ n as 



n 


I X Y Z 




00 01 11 10 



Suppose u,v G (Z2) . Let us employ the shorthand u = [z\x] and v = [z'\x'] where z, x, z', 
x' G Z2. For example, suppose u= [0|1]. Then A'' (u) = X. The map N induces an isomorphism 
[N] : (^2)^ — > [n] because addition of vectors in (^2)^ is equivalent to multiplication of Paulis up 
to a global phase: 

[N{u + v)] = [N{u)][N{v)]. (2.13) 

We name this isomorphism the P2B isomorphism (for Pauli-to-binary isomorphism). 
Let denote the symplectic product between two elements u,v G (^2)^: 

uQv = zx'-xz'. (2.14) 

The symplectic product gives the commutation relations of elements of H: 

N (w) N {v) = (-1)(''0'') TV {v) N («) . 

The symplectic product and the P2B isomorphism [A''] thus give a useful way to phrase Pauli 
relations in terms of binary algebra. 

The extension of the above definitions and the P2B isomorphism [N] to multiple qubits is 
straightforward. Let A = Ai • • • j4„ denote an arbitrary element of H". We can similarly 
define the phase-free n-qubit Pauli group [H"] = {[A] | A e H"} where 

[A] = {/3A|/3eC, |/3| = 1}. (2.15) 

The group operation * for the above equivalence class is as follows: 

[A] * [B] = [Ai] * [Bi] • • • * [Br,] = [AiBi] • • • [A„B„] = [AB] . (2.16) 
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The equivalence class [11"] forms a commutative group under operation *. Consider the 2n- 
dimensional vector space 

(Z2)'" = {(z,x):z,xG(Z2)"}. (2.17) 

It forms the commutative group ((Z2)^" , +) with operation + defined as binary vector addition. 
We employ the notation u = [z|x] and v = [z'|x'] to represent any two vectors u, v e (^2)^" 
respectively. Each vector z and x has elements {zi,...,Zn) and (si, . . . ,a;„) respectively with 
similar representations for z' and x'. The symplectic product of u and v is 

n n 

u0 V =^ Zjcc'j - Xi^;- = Wi, (2.18) 

i=l i=l 

where Ui = [zi\xi] and Vi = [zHx'^]. Let us define a map N : (Z2)^" H" as follows: 

N (u) = TV (mi) (g) • • • (g) TV (u„) . (2.19) 

Let 

X(x) =X"'i (g)---0X"=", Z(z) = (g) ••• (g)Z^", (2.20) 
so that N (u) and Z (z) X (x) belong to the same equivalence class: 

[N(u)] = [Z(z)X(x)]. (2.21) 

The map [N] : (Z2)^" — > [H"] is an isomorphism for the same reason given in (2.13): 

[N(u + v)] = [N(u)][N(v)], (2.22) 

where u, v e (Z2)^". The symplectic product captures the commutation relations of any operators 
N(u) and N(v): 

N (u) N (v) = (-1)("®^) N (v) N (u) . (2.23) 

The above P2B isomorphism and symplectic algebra are useful in making the relation between 
classical linear error correction and quantum error correction more explicit. 
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2.5 Example 



An example of a stabilizer code is the five qubit [5,1] stabilizer code [75, 66]. It encodes k = 1 
logical qubit into n = 5 physical qubits and protects against an arbitrary single-qubit error. Its 
stabilizer consists of n — A; = 4 Pauli operators: 



91 


= X 


z 


z 


X 


I 


92 


= I 


X 


z 


z 


X 


93 


= X 


I 


X 


z 


z 


9i 


= z 


X 


I 


X 


z 



(2.24) 



The above operators commute. Therefore the codespace is the simultaneous +l-eigenspace of the 
above operators. Suppose a single-qubit error occurs on the encoded quantum register. A single- 
qubit error is in the set {Xj, Yi, Zi} where denotes a Pauli error on qubit i. It is straightforward 
to verify that any arbitrary single-qubit error has a unique syndrome. The receiver corrects any 
single-qubit error by identifying the syndrome and applying a corrective operation. 



2.6 Closing Remarks 

The mathematics developed in this chapter form the basis for the work in later chapters. In 
particular, the binary representation of quantum codes is important. The role of the "twisted" 
symplectic product plays an important role for the entanglement-assisted codes of the next chap- 
ter because it helps in determining the commutation relations of an arbitrary (not necessarily 
commuting) set of Pauli generators. 
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Chapter 3 



Entanglement- Assisted Stabilizer Quantum Codes 



Classical to quantum made dinero, 
But the twisted strange product was zero, 
Little did we know, 
What Entanglement could show. 
And again become quantum's good hero. 

The cntanglcnicnt-assistcd stabilizer formalism is a significant extension of the standard stabilizer 
formalism that incorporates shared entanglement as a resource for protecting quantum information 
[35, 34]. The advantage of entanglement-assisted stabilizer codes is that the sender can exploit 
the error-correcting properties of an arbitrary set of Pauli operators. The sender's PauU operators 
do not necessarily have to form an abelian subgroup of 11". The sender can make clever use of her 
shared ebits so that the global stabilizer is abelian and thus forms a valid quantum error-correcting 
code. Figure 3.1 demonstrates the operation of a generic entanglement-assisted stabilizer code. 

Several references provide a review of this technique and generalizations of the basic theory 
to block entanglement distillation [76], continuous-variable codes [70], and entanglement-assisted 
operator codes for discrete- variable [30, 31] and continuous- variable systems [68]. Chapters 5, 6, 
7, 8, and 9 of this thesis extend the entanglement-assisted stabilizer formalism to many different 
"convolutional" coding scenarios. We explain what we mean by "convolutional" in the next 
chapter. 
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For now, wc provide a review of the basics of entanglement-assisted coding. This chapter also 
gives two original contributions: a method to determine the encoding and decoding circuit for 
an entanglement-assisted code and several formulas that determine the number of ebits that an 
arbitrary entanglement-assisted block code requires. 



3.1 Review of the Entanglement- Assisted Stabilizer For- 
malism 

The fundamental unit of bipartite entanglement is the ebit. We must first understand the commu- 
tativity properties of the ebit before developing the entanglement-assisted formalism. We express 
the state |$+) of an ebit shared between a sender Alice and a receiver Bob as follows: 

. m'^'^^^r. (3.1) 

The two operators that stabilize this ebit state are X^X^ and Z^Z^ . These two operators 
commute, 

{X^X^ ,Z^Z^\ = 0, 
but the local operators (operating only on either party) anticommute, 

{X^,Z^} = {X^,Z^}=0. 

The above commutation relations hint at a way that we can resolve anticommutativity in a 
set of generators. Suppose that we have two generators that anticommute. We can use an ebit 
of entanglement to resolve the anticommutativity in the two generators. We explain this idea in 
more detail in what follows. 

We review the general construction of an entanglement-assisted code. An [n, k; c] entanglement- 
assisted code employs c ebits and a ancilla qubits to encode k information qubits. Suppose that 
there is a nonabelian subgroup <S C 11" of size 2c -|- a. Application of the fundamental theorem of 
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symplcctic geometry* [77] (Lemma 1 in [34]) states that there exists a minimal set of independent 
generators {Zi, . . . , Za_|_c, X^+i, . . . ,Xa+c} for iS with the following commutation relations: 

Vz,j [Z„Z,]=0, Vi,j [X^,Xj]=0, (3.2) 

yij^j [Xi,Zj]=o, yi {Xi,Zi}=:0. 

The decomposition of S into the above minimal generating set determines that the code requires 
a ancilla qubits and c ebits. The parameters a and c generally depend on the set of generators in 
S. 

There exists a symplectic Gram-Schmidt orthogonalization procedure that gives the decompo- 
sition [35, 34, 51, 78]. Specifically, the algorithm performs row operations (multiplication of the 
Pauli generators) that do not affect the code's error-correcting properties and thus gives a set of 
generators that form an equivalent code. The decomposition also minimizes the number of ebits 
required for the code and we prove this optimality in Section 3.6. 

We present a simple stabilizer version of the symplectic Gram-Schmidt algorithm. Suppose we 
have a set of m generators gi, . . . , g^- Consider gi. It either commutes with all other generators 
or it anticommutes with at least one other generator. Remove it from the set if it commutes with 
all other generators. Now suppose that it does not, i.e., it anticommutes with one other generator 
gi. Relabel 52 as gi and vice versa. Recall that we can multiply generators without changing 
the error-correcting properties of the set of generators. So we perform several multiplications 
to change the commutation relations to have the standard commutation relations in (3.2). We 
perform the following manipulations on the generators gs, . . . , gm- 

9j=9j-9(^'''''^-9i^'"''^ VzG{3,...,m} 

where the function / is equal to zero if its two arguments commute and it is equal to one if its 
two arguments anticommute. (This function is the symplectic product.) Remove the generators 
gi and 32 from the set. Repeat the algorithm on the remaining generators. When the algorithm 



*We loosely refer to this theorem as the fundamental theorem of symplectic geometry because of its importance 
in symplectic geometry and in quantum coding theory. 
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Encoding 




Figure 3.1: The operation of an entanglement-assisted quantum error-correcting code. The sender 
encodes quantum information in state \ip) with the help of local ancilla qubits |0) and her half of 
a set of shared ebits | <&"''). She then sends her qubits over a noisy quantum channel. The channel 
does not corrupt the receiver's half of the set of shared ebits. The receiver performs multi-qubit 
measurements on all of the qubits to diagnose the channel error. He performs a recovery unitary 
R to reverse the estimated channel error. 

finishes, all the "removed' generators constitute the generators for the code and have the standard 
commutation relations in (3.2). 

We can partition the nonabelian group S into two subgroups: the isotropic subgroup 5/ 
and the entanglement subgroup Se- The isotropic subgroup Si is a commuting subgroup of 
S and thus corresponds to ancilla qubits: Sj = . . . , Zq}. The elements of the entan- 

glement subgroup Se come in anticommuting pairs and thus correspond to halves of ebits: 
Se ^ {Za+i, ■ ■ ■ , Za+c, ^a+i, ■ ■ ■ , -^a+c}- The two subgroups Sj and Se play a role in the error- 
correcting conditions for the entanglement-assisted stabilizer formalism. An entanglement-assisted 
code corrects errors in a set £ C 11" if 

yEi,E2 e £ e\E2 e Si or e\E2 G H" - Z {{Si.Se)) ■ 

The conditions correspond to error pairs Ei , E2 in an error set £ . The first condition corresponds 
to the passive error-correcting capability of the code, and the second condition corresponds to its 
active error-correcting capability. 
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Figure 3.1 illustrates the operation of an entanglement-assisted stabilizer code. The following 
four steps explain the operation of an [n, k; c] entanglement-assisted quantum code: 

(i). The sender and receiver share c ebits before quantum communication begins and the sender 
has a ancilla qubits. The unencoded state is a simultaneous +l-eigenstate of the following 
operators: 

{Za+l\Zi, . . . , Za+c\Zc, -'^a-l-ll-'^l, • • • , ^a+c\Xc, Zi, . . . , Za} ■ (3-3) 

The operators to the right of the vertical bars indicate the receiver's half of the shared ebits. 
The sender encodes her k information qubits with the help of a ancilla qubits and her half 
of the c ebits. The encoding operation is U in Figure 3.1. The encoding unitary transforms 
the unencoded operators to the following encoded operators: 

{Za+l\Zl, . . . , Za+c\Zc, -^a-|-l|^l) • • • ) ^a+c\Xc, Z\, . . . , Za} ■ (3-4) 



(ii) . The sender sends her n qubits over a noisy quantum communication channel. The noisy 

channel affects these n qubits only and does not affect the receiver's half of the c ebits. 

(iii) . The receiver combines his half of the c ebits with those he receives from the noisy quantum 

channel. He performs measurements on all n + c qubits to diagnose an error that may occur 
on the n qubits. 

(iv) . After estimating which error occurs, the receiver performs a recovery operation that reverses 

the estimated error. 



3.2 The Symplectic Product Matrix 

Let us write the local operators and Z^ on the sender Alice's side as respective symplectic 
vectors hi and /12: 



hi= 1 



2 



1 • 
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The "symplectic product matrix" fl of these two symplectic vectors is as follows: 

n = . (3.5) 



1 

1 



We call this matrix the symplectic product matrix because its entries enumerate the symplectic 
products: 

Mij = hiQhj Vi,ie{l,2}. (3.6) 

The symplectic product matrix in (3.5) is so special for the purposes of this thesis that we give it 
the name J = il. It means that two generators have commutation relations that are equivalent to 
half of an ebit. 

A general set of generators for a quantum block code can have complicated commutation 
relations. For a given set of generators, the commutation relations for the generators resulting from 
the symplectic Gram-Schmidt orthogonalization procedure have a special form. Their symplectic 
product matrix fl has the standard form: 

f2 = 0J©0[O] (3.7) 
i=i j=i 

where the large and small ® correspond to the direct sum operation and [0] is the one-element null 
matrix. The standard form above implies that the commutation relations of the 2c + a generators 
are equivalent to those of c halves of ebits and a ancilla qubits — ^the commutation relations are 
equivalent to those in (3.2). 



3.3 The Rate of an Entanglement- Assisted Quantum Code 

We can interpret the rate of an entanglement-assisted code in three different ways [34, 35, 51]. 
Suppose that an entanglement-assisted quantum code encodes k information qubits into n physical 
qubits with the help of c ebits. 
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(i) . The "entanglement-assisted" rate assumes that entanglement shared between sender and 

receiver is free. Bennett et al. make this assumption when deriving the entanglement- 
assisted capacity of a quantum channel for sending quantum information [40, 41]. The 
entanglement-assisted rate is fc/n for a code with the above parameters. 

(ii) . The "trade-off' rate assumes that entanglement is not free and a rate pair determines perfor- 

mance. The first number in the pair is the number of noiseless qubits generated per channel 
use, and the second number in the pair is the number of ebits consumed per channel use. The 
rate pair is (fc/n, c/n) for a code with the above parameters. Quantum information theorists 
have computed asymptotic trade-off curves that bound the rate region in which achievable 
rate pairs lie [61]. Brun et al.'s construction for an entanglement-assisted quantum block 
code minimizes the number c of ebits given a fixed number k and n of respective information 
qubits and physical qubits [35, 34]. 

(iii) . The "catalytic rate" assumes that bits of entanglement are built up at the expense of trans- 

mitted qubits [35, 34]. A noiseless quantum channel or the encoded use of noisy quantum 
channel are two different ways to build up entanglement between a sender and receiver. The 
catalytic rate of an [n, k; c] code is {k — c) /n. 

Which interpretation is most reasonable depends on the context in which we use the code. In 
any case, the parameters n, k, and c ultimately govern performance, regardless of which definition 
of the rate we use to interpret that performance. 

3.4 Example of an Entanglement- Assisted Code 

We present an example of an entanglement-assisted code that corrects an arbitrary single-qubit 
error [34]. This example highlights the main features of the theory given above. Suppose the 
sender wants to use the quantum error-correcting properties of the following nonabelian subgroup 

of n^: 

Z X Z I 
Z Z I Z 

(3.8) 

X Y X I 
XXIX 
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The first two generators anticommute. We obtain a modified third generator by multiplying the 
third generator by the second. We then multiply the last generator by the first, second, and 
modified third generators. The error-correcting properties of the generators are invariant under 
these operations. The modified generators are as follows: 



51 

92 

9s 
9i 



Z X Z I 

Z Z I z 

Y X X Z 

Z Y Y X 



(3.9) 



The above set of generators have the commutation relations given by the fundamental theorem of 
symplectic geometry: 

{91,92} = [91,93] = [91, 9i] = [92,93] = [92, 9i] = [93, 9i] = 0. 
The above set of generators is unitarily equivalent to the following canonical generators: 



XIII 
Z I I I 
I Z I I 
I I Z I 



(3.10) 



We can add one ebit to resolve the anticommutativity of the first two generators: 



XIII 
Z I I I 
I Z I I 
I I Z I 



X 

z 
I 
I 



(3.11) 



A -|-l-eigenstate of the above stabilizer is the following state 



l$+/^loo)^lV')^ 



(3.12) 
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where 1-0) is a qubit that the sender wants to encode. A local encoding unitary then rotates the 
generators in (3.11) to the following set of globally commuting generators; 



z 


X 


Z 


I 


X 


z 


Z 


I 


Z 


z 


Y 


X 


X 


Z 


I 


Z 


Y 


Y 


X 


I 



(3.13) 



The receiver measures the above generators upon receipt of all qubits to detect and correct errors. 



3.5 Encoding Algorithm 

Suppose we have an arbitrary set of Pauli matrices in 11" whose error-correcting properties we 
would like to exploit. The algorithm in this section (in Ref. [51]) determines the encoding and 
decoding circuit for the set of Pauli generators. 

It has two additional benefits. We do not necessarily know beforehand how many ebits we 

require for the Pauli matrices to form a commuting set. Several methods exist such as the Gram- 
Schmidt algorithm outlined in Section 3.1 and the others in Refs. [35, 34, 30, 31, 51], but the 
algorithm here also determines the optimal number of ebits and the measurements the receiver 
performs to diagnose errors. It "kills three birds with one stone." 

We continue with the example in Section 3.4. We detail an algorithm for determining an 
encoding circuit and the optimal number of ebits for the entanglement-assisted code. The operators 
in (3.8) have the following representation as a binary matrix according to the P2B isomorphism: 



H 



10 10 
110 1 
10 




10 

1110 
110 1 



(3.14) 



Call the matrix to the left of the vertical bar the "Z" matrix and the matrix to the right of the 
vertical bar the "X" matrix. 
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The algorithm consists of row and cohinin operations on the above matrix. Row operations 
do not affect the error-correcting properties of the code but are crucial for arriving at the optimal 
decomposition from the fundamental theorem of symplectic geometry. The operations available 
for manipulating columns of the above matrix are Clifford operations (discussed at the end of 
Section 2.2). The operations have the following effect on entries in the binary matrix: 

(i) . A CNOT gate from qubit i to qubit j adds column i to column j in the X matrix and adds 

column j to column i in the Z matrix. 

(ii) . A Hadamard gate on qubit i swaps column i in the Z matrix with column i in the X matrix 

and vice versa. 

(iii) . A phase gate on qubit i adds column i in the X matrix to column i in the Z matrix. 

(iv) . Three CNOT gates implement a qubit swap operation [32]. The effect of a swap on qubits 

i and j is to swap columns i and j in both the X and Z matrix. 

The algorithm begins by computing the symplectic product between the first row and all other 
rows. Leave the matrix as it is if the first row is not symplectically orthogonal to the second row 
or if the first row is symplectically orthogonal to all other rows. Otherwise, swap the second row 
with the first available row that is not symplectically orthogonal to the first row. In our example, 
the first row is not symplectically orthogonal to the second so we leave all rows as they are. 

Arrange the first row so that the top left entry in the X matrix is one. A CNOT, swap, 
Hadamard, or combinations of these operations can achieve this result. We can have this result 
in our example by swapping qubits one and two. The matrix becomes 



110 
110 1 
10 




10 

1110 
110 1 



(3.15) 



Perform CNOTs to clear the entries in the X matrix in the top row to the right of the leftmost 
entry. These entries are already zero in this example so we need not do anything. Proceed to the 
clear the entries in the first row of the Z matrix. Perform a phase gate to clear the leftmost entry 
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in the first row of the Z matrix if it is equal to one. It is equal to zero in our example so we need 
not do anything. We then use Hadamards and CNOTs to clear the other entries in the first row 
of the Z matrix. 

For our example, perform a Hadamard on qubits two and three. The matrix becomes 




10 1 
1110 
10 



1110 
10 
10 
10 1 



(3.16) 



Perform a CNOT from qubit one to qubit two and from qubit one to qubit three. The matrix 
becomes 




10 1 
1110 
110 



10 
10 
1110 
1111 



(3.17) 



The first row is complete. We now proceed to clear the entries in the second row. Perform a 
Hadamard on qubits one and four. The matrix becomes 



10 

1110 
110 1 




110 1 
1110 
1110 



(3.18) 



Perform a CNOT from qubit one to qubit two and from qubit one to qubit four. The matrix 
becomes 



10 

110 
110 1 





10 

10 11 

10 11 



(3.19) 



The first two rows are now complete. They need one ebit to compensate for their anticommuta- 
tivity or their nonorthogonality with respect to the symplectic product. 
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H 



H 



Figure 3.2: Encoding circuit for the entanglement- assisted code in the example of Section 3.4. 
The "H" gate is a Hadamard gate and the "P" gate is a phase gate. 



Now we perform row operations that are similar to the "symplectic Gram-Schmidt orthogo- 
nalization." Add row one to any other row that has one as the leftmost entry in its Z matrix. Add 
row two to any other row that has one as the leftmost entry in its X matrix. For our example, we 
add row one to row four and we add row two to rows three and four. The matrix becomes 



10 

110 
10 1 





10 

11 

11 



(3.20) 



The first two rows are now symplectically orthogonal to all other rows per the fundamental theorem 
of symplectic geometry. 

We proceed with the same algorithm on the next two rows. We can deal with the next two 
rows individually because they are symplectically orthogonal to each other. Perform a Hadamard 
on qubit two. The matrix becomes 



10 



10 

1 




10 
111 
111 



(3.21) 
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Perform a CNOT from qubit two to qubit three and from qubit two to qubit four. The matrix 
becomes 



(3.22) 



1 



































1 














1 


1 








1 











1 





1 





1 









Perform a phase gate on qubit two: 



1 



































1 

















1 








1 

















1 





1 









(3.23) 



Perform a Hadamard on qubit three followed by a CNOT from qubit two to qubit three: 



10 





1 




10 
10 
110 



(3.24) 



Add row three to row four and perform a Hadamard on qubit two: 



10 



10 

1 





10 



10 



(3.25) 
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Perform a Hadamard on qubit four followed by a CNOT from qubit three to qubit four. End by 
performing a Hadamard on qubit three: 

1000 0000 

0000 1000 

1 

1 

The above matrix now corresponds to the canonical PauUs in (3.10). Adding one half of an ebit 
to the receiver's side gives the canonical stabilizer in (3.11) whose simultaneous +l-eigenstate is 
the state in (3.12). 

Figure 3.2 gives the encoding circuit corresponding to the above operations. The above oper- 
ations in reverse order take the canonical stabilizer (3.11) to the encoded stabilizer (3.13). 

3.6 Optimal Entanglement Formulas 

Entanglement is a valuable resource, and we would like to minimize the amount that a sender and 
receiver need to consume for entanglement-assisted quantum coding. Hsieh, Devetak, and Brun 
first addressed this issue by determining a useful formula that gives the optimal number of ebits 
required by a Calderbank-Shor-Steane (CSS) entanglement-assisted code [31]. In particular, the 
number c of ebits that a CSS entanglement-assisted code requires is 

c = rank {HH'^) , (3.27) 

where H corresponds to the parity check matrix of a classical binary block code that we import 
to correct quantum bit and phase flips. The same authors also determined an upper bound for 
the number of ebits that an entanglement-assisted quantum LDPC code requires [79] . 

In this section, we present several generalizations of the above formula. Our first theorem 
gives a formula for the optimal number of ebits that an arbitrary (non-CSS) entanglement-assisted 
quantum block code requires. We find special cases of this formula that apply to an entanglement- 
assisted quantum block code produced from two arbitrary classical binary block codes, and to 



(3.26) 
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codes from a classical block code over Gi^(4). Our last formula applies to continuous-variable 
entanglement-assisted codes [70]. 

Theorem 3.6.1. Suppose we want to build an entanglement- assisted quantum code from genera- 
tors corresponding to the rows in a quantum check matrix 



H= H 



H 



X 



(3.28) 



where H is an {n — k) x2n- dimensional binary matrix representing the quantum code, and both Hz 
and Hx are {n — k) x n-dimensional binary matrices. Then the resulting code is an [[n, fc + c; c]] 
entanglement- assisted code and requires c ebits, where 

c = rank {HxH^ + HzHI) /2, (3.29) 

and addition is binary. 

Proof. Consider the symplectic product matrix Q.h = HxH^ + HzH^. Its entries are the sym- 
plcctic products between all rows of H so that 



pH]ij = hiQ hj, (3.30) 

where hi is the i**^ row of H. The matrix fijj is a (n — fc) x (n — A;)-dimensional binary matrix. 

Refs. [35, 34] and the algorithm in Section 3.1 outline a symplectic Gram-Schmidt orthogonal- 
ization procedure (SGSOP) that uniquely determines the optimal (i.e., minimal) number of ebits 
that the code requires, and Ref. [36] proves that the SGSOP gives the optimal number of ebits. 

The code construction in Refs. [35, 34] shows that the resulting cntanglcmcnt-assistcd quantum 
code requires at most c ebits. The essence of the argument in Ref. [36] is that the resulting 
entanglement-assisted quantum code requires at least c ebits because any fewer ebits would not 
be able to resolve the anticommutativity of the generators on Alice's side of the code. 

The SGSOP performs row operations that do not change the error-correcting properties of the 
quantum code (because the code is additive), but these row operations do change the symplectic 
product relations. These row operations are either a row swap S {i,j), where S {i,j) is a full-rank 
{n — k) X (n — k) matrix that swaps row i with j, or a row addition A{i,j), where A{i,j) is a 
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full-rank (n — k) x (n — k) matrix that adds row i to row j. These row operations multiply the 
matrix H from the left. The SGSOP then is equivalent to a full-rank (n — k) x (n — k) matrix 
G that contains all of the row operations and produces a new quantum check matrix H' = GH 
with corresponding symplectic product matrix ^Ih' — G (^HxH^ + HzHjr^ G^ . In particular, the 
resulting symplectic product matrix is in the standard form in (3.7) so that 

c n—k—2c 

fJi/' = 0J® [0]. (3.31) 

Each matrix J in the direct sum corresponds to half of an ebit as described in Section 3.2 and 
has rank two. Each matrix [0] has rank zero and corresponds to an ancilla qubit. The optimal 
number of ebits required for the code is rank {Q.h') /2: 

(c 7i—k—2c \ c n—k — 2c 

0J® [0] = ^ rank ( J) + ^ rank ([0]) = 2c. (3.32) 
i=i j=i J i=i j=i 

The second equality follows because the rank of a direct sum is the sum of the individual matrix 
ranks, and the third equality follows from the individual matrix ranks given above. The number c 
of ebits is also equal to rank (CIh) /2 because the matrix G is full rank. The code is an [[n, k + c; c]] 
entanglement-assisted quantum block code by the construction in Refs. [35, 34]. 

□ 

Our formula (3.29) is equivalent to the formula at the top of page 14 in Ref. [35], but it provides 
the quantum code designer with a quick method to determine how many ebits an entanglement- 
assisted code requires, by simply "plugging in" the generators of the code. 

The formula (3.29), like the CSS formula in (3.27), is a measure of how far a set of generators 
is from being a commuting set, or equivalently, how far it is from giving a standarcl stabilizer code. 

Corollary 3.6.1 below gives a formula for the optimal number of ebits required by a CSS 
entanglement-assisted quantum code. It is generally a bit less difficult to compute than the above 
formula in (3.29). This reduction in complexity occurs because of the special form of a CSS 
quantum code and because the size of the matrices involved are generally smaller for a CSS code 
than for a general code with the same number of generators and physical qubits. 
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Corollary 3.6.1. Suppose we import two classical [n,ki,di\ and [n, fc2,d2] binary codes with 
respective parity check m,atrices Hi and H2 to build an entanglement-assisted quantum code. The 
resulting code is an [[n, ki + k2 — n + c, min ((^1,^2) ; c]] entanglement- assisted code, and requires c 
ebits where 

c = rank (HiH^) . (3.33) 
Proof. The quantum check matrix has the following form: 



H = 



The symplectic product matrix Qh is then 



Hi 








H2 



(3.34) 



Hi 




















Hi 


+ 


H2 




Hi 




H2HI 



HIH2 




(3.35) 



The above matrix is equivalent by a full rank permutation matrix to the matrix i/ii/J © H2HI , 
so the rank of ^Ih is 

rank (Q^) = rank (ifiiJj ® H2H{) = rank (HiH^) + rank {H2HI) = 2 rank (HiH^) (3.36) 

The second equality follows because the rank of a direct sum is equivalent to the sum of the 
individual ranks, and the third equality follows because the rank is invariant under matrix trans- 
position. The number of ebits required for the resulting entanglement-assisted quantum code is 



rank{HiH^ 



ig the result of the previous theorem. The construction in Refs. [35, 34] produces 



an [[n, A;i + /c2 — n + c, min (^1,(^2) ; c]] entanglement-assisted quantum block code. 



□ 



Corollary 3.6.2. Suppose we import an [n,k,d\^ classical code over GF (A) with parity check 
matrix H for use as an entanglement- assisted quantum code according to the construction in 
Refs. [35, 34]- Then the resulting quantum code is an [[n,2k — n-\- c;c]] entanglement- assisted 
quantum code where c = rank (hW) and f denotes the conjugate transpose operation over ma- 
trices in GF (4) . 
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Proof. Rcf. [34] shows how to produce an ciitaiigicrncnt-assistcd quantum code from the parity 
check matrix H of a classical code over GF (4) . The resulting quantum parity check matrix Hq 
in symplectic binary form is 

\luH 



Q = l 



ujH 



(3.37) 



where 7 denotes the isomorphism between elements of GF (4) and symplectic binary vectors that 
represent Pauli matrices. We augment the table in (2.12) to include the entries from GF(4): 



n 


/ 


X 


Y 


Z 




00 


01 


11 


10 


GF (4) 





w 


1 


w 



(3.38) 



The isomorphism 7 is the mapping from the third row to the second row in the above table. 
The symplectic product between binary vectors is equivalent to the trace product of their GF (4) 
representations (see, e.g., Ref. [34]): 



hi hj = tr 



(3.39) 



where hi and hj are any two rows of Hq, ■ denotes the inner product, the overbar denotes the 
conjugate operation, and tr{a;} = x + x denotes the trace operation over elements of GF (4). We 
exploit these correspondences to write the symplectic product matrix ^Hq for the quantum check 
matrix Hq as follows: 



Q.Hq = tr 



= tr 









V 






LjjH 




LOH 




> = tr 1 


uH 


u)H 




u)H 


> 




ujH 



hW 


u)HH^ 




1 Co 






H 




LoHH^ 


HW 




CO 1 



1 UJ 




1 UJ 


(g)HH^ 


UJ 1 




U! 1 





(3.40) 
(3.41) 
(3.42) 
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where the "tr" operation above is an element-wise trace operation over GF (4) (it is not the matrix 
trace operation.) The matrix ^Ihq is over the field GF(2), but we can consider it as being over 
the field GF (4) without changing its rank. Therefore, we can multiply it by matrices over the 
field GF (4) . Consider the following full-rank GF (4) matrices: 



1 Uj 

1 



1 
1 1 



(3.43) 



We premultiply and postmultiply the matrix ^Irq as follows and obtain a matrix with the same 
rank as flnQ- 



A2AinHQA\Al 











1 




















1 













HH^ 




HH^®HH^ (3.44) 



Therefore, the rank of is 

rank {Q.Hq) = rank {HH'^ ® HH'^) = rank {HH'^) + rank {HH'^) = 2 rank {HH'^) . (3.45) 



The second equality holds because the rank of a direct sum is the sum of the individual ranks and 
the third holds because the rank is invariant under the matrix transpose operation. Therefore, 
the resulting entanglement-assisted quantum code requires c = rank(i?iJ'l') ebits by applying the 
result of the original theorem. The construction in Refs. [35, 34] produces an [[n, 2k — n + c; c]] 
entanglement-assisted quantum code. 

□ 

Corollary 3.6.3. We can construct a continuous-variable entanglement- assisted quantum code 



H 



X 



where H 



from generators corresponding to the rows in quantum check matrix H 
is {n — k) X 2n-dimensional, H is a real matrix representing the quantum code [70], and both Hz 
and Hx are {n — k) x n- dimensional. The resulting code is an [[n, k + c; c]] continuous-variable 
entanglement-assisted code and requires c entangled modes where 



c = rank {HxH^ - HzHl) /2. 



(3.46) 
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Proof. The proof is similar to the proof of the first theorem but requires manipulations of real 
vectors instead of binary vectors. See Ref. [70] for details of the symplectic geometry required for 
continuous-variable entanglement-assisted codes. □ 

Remark 3.6.1. A similar formula holds for entanglement-assisted qudit codes by replacing the 
subtraction operation above with subtraction modulo d. Specifically, we can construct a qudit 
entanglement-assisted quantum code from generators corresponding to the rows in check matrix 
whose matrix entries are elements of the finite field 1^4. The code requires c 
edits (a ci-dimensional state (^X]f=o K) K)) where 

c = rank {HxH^ Qd HzHJI) /2 

and Qd is subtraction modulo d. We use subtraction modulo d because the symplectic form over 
c?-dimensional variables includes subtraction modulo d. 

3.7 Closing Remarks 

This chapter reviewed the theory of entanglement-assisted coding for block codes. We showed 
how to manipulate these codes with both a Pauli representation and a binary one. Our con- 
tributions included an algorithm for computing an encoding circuit and a method to determine 
the optimal number of ebits for an entanglement-assisted code. Much of this thesis extends the 
entanglement-assisted technique to the domain of quantum convolutional coding. We introduce 
quantum convolutional coding in the next chapter. 



H = 



H 



X 
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Quantum Convolutional Codes 



Chapter 4 



Forney, Viterbi, and Elias, 
To Andy the Trojans are pious, 
Their m,ethods, we want 'em, 
Who'd think they'd go quantum? 
So now they still stupefy us. 

The block codes in the previous two chapters are useful in quantum computing and in quantum 
communication. One of the drawbacks of the block-coding technique for quantum communication 
is that, in general, the sender must have all her qubits ready before encoding takes place. For a 
large block code, this preparation may be a heavy demand on the sender. 

Quantum convolutional coding theory [43, 44, 48, 1] offers a different paradigm for coding 
quantum information. The convolutional structure is useful in a quantum communication scenario 
where a sender possesses a stream of qubits to send to a receiver. The encoding circuit for a 
quantum convolutional code has a much lower complexity than an encoding circuit needed for a 
large block code. It also has a repetitive pattern so that the same physical devices or the same 
routines can encode the stream of quantum information in an online fashion as soon as information 
qubits are available for encoding. 

Quantum convolutional stabilizer codes borrow heavily from the structure of their classical 
counterparts [43, 44, 48, 1]. Quantum convolutional codes are similar because some of the qubits 
feed back into a repeated encoding unitary and give the code a memory structure like that of a 
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classical convolutional code. The quantum codes feature online encoding and decoding of qubits. 
This feature gives quantum convolutional codes both their low encoding and decoding complexity. 

Our techniques for encoding and decoding are an expansion of previous techniques from quan- 
tum convolutional coding theory. Previous techniques for encoding and decoding include finite- 
depth operations only. A finite-depth operation propagates errors to a finite number of neighboring 
qubits in the qubit stream. We introduce an infinite-depth operation to the set of shift-invariant 
Clifford operations and explain it in detail in Section 4.6. We must be delicate when using infinite- 
depth operations because they can propagate errors to an infinite number of neighboring qubits 
in the qubit stream. We explain our assumptions in detail when we include infinite-depth oper- 
ations in entanglement-assisted quantum convolutional codes in Section 5.4 of the next chapter. 
An infinite-depth operation gives more fiexibility when designing encoding circuits — similar to the 
way in which an infinite-impulse response filter gives more fiexibility in the design of classical 
convolutional circuits. 

We structure this chapter as follows. In Section 4.1, we introduce the definition of a quantum 
convolutional code and discuss how it operates. Section 4.2 introduces the Pauli-to-binary isomor- 
phism. This mapping simplifies the mathematics of a quantum convolutional code by showing how 

to represent it with a matrix of binary polynomials. Section 4.3 discusses the important "shifted 
symplectic product" that gives the commutation relations of a set of convolutional generators. We 
discuss in Section 4.4 how to manipulate the binary polynomial matrix that represents a quantum 
convolutional code. In Sections 4.5 and 4.6, we respectively review finite-depth and infinite-depth 
operations for encoding a quantum convolutional code. 

4.1 Review of the Convolutional Stabilizer Formalism 

We review the theory of convolutional stabilizer codes by considering a set of Pauli matrices that 
stabilize a stream of encoded qubits. We first give the mathematical definition of a quantum 
convolutional code and follow by discussing the various steps involved in the operation of it. 
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4.1.1 Quantum Convolutional Code Definition 

A quantum convolutional stabilizer code acts on a Hilbert space H that is a countably infinite 
tensor product of two-dimensional qubit Hilbert spaces {^1}-^^+ where 

00 

n = (^Hi. (4.1) 

and Z+ = {0, 1, . . .}. A sequence A of Pauli matrices {Ai}^^^+, where 

00 

A = (g) ^i, (4.2) 

can act on states in H. Let 11^^ denote the set of all Pauli sequences. The support supp(A) of a 
Pauli sequence A is the set of indices of the entries in A that are not equal to the identity. The 

weight of a sequence A is the size |supp (A)| of its support. The delay del(A) of a sequence A is 
the smallest index for an entry not equal to the identity. The degree deg(A) of a sequence A is 
the largest index for an entry not equal to the identity. E.g., the following Pauli sequence 

I X I Y Z I I ■■■ , (4.3) 

has support {1,3, 4}, weight three, delay one, and degree four. A sequence has finite support if its 
weight is finite. Let F(n^^) denote the set of Pauli sequences with finite support. The following 
definition for a quantum convolutional code utilizes the set F{I1^^) in its description. 

Definition 4.1.1. A rate k/n- convolutional stabilizer code with < k < n is specified by a 
commuting set Q of all n-qubit shifts of a basic generator set Qo- The commutativity requirement 
is necessary for the same reason that standard stabilizer codes require it [32]. The basic generator 
set Qo has n — k Pauli sequences of finite support: 

^?o = {GieF(n^''):l<i<n-fc}. (4.4) 
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The constraint length v the code is the maximum degree of the generators inQo- ^ frame of the 
code consists of n qubits. The definition of a quantum convolutional code as n-qubit shifts of the 
basic set Go is what gives the code its periodic structure. 

A quantum convolutional code admits an equivalent definition in terms of the delay transform 
or £)-transform. The Z)-transform captures shifts of the basic generator set Qo- Let us define the 
?>qubit delay operator D acting on any Pauli sequence A e 11^^ as follows: 

D (A) = 7®" O A. (4.5) 

We can write j repeated applications of £> as a power of D: 

(A) = 7®^" A. (4.6) 

Let {Qq) be the set of shifts of elements of Qq by j. Then the full stabilizer Q for the convolutional 
stabilizer code is 

Q= [j {Qo) . (4.7) 

Example 4.1.1. Forney et al. provided an example of a rate-1/3 quantum convolutional code 
by importing a particular classical quaternary convolutional code [48, 1]. Grassl and Rotteler 
determined a noncatastrophic encoding circuit for Forney et al.'s rate-1/3 quantum convolutional 
code [45]. The basic stabilizer and its first shift are as follows: 
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XXX 


XZY 


III 


III 


777 


zzz 


ZYX 


III 


III 


777 


III 


XXX 


XZY 


III 


777 


III 


ZZZ 


ZYX 


III 



(4.8) 



The code consists of all three-qubit shifts of the above generators. The vertical bars are a visual 
aid to illustrate the three-qubit shifts of the basic generators. The code can correct for an arbitrary 
single-qubit error in every other frame. 
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Figure 4.1: The operation of a quantum convolutional code. The sender appHes the same unitary 
successively to a stream of information qubits and ancilla qubits. The convolutional structure 
implies that the unitary overlaps some of the same qubits. The sender transmits her qubits as 
soon as the unitary finishes processing them. The noisy quantum channel corrupts the transmitted 
qubits. The receiver performs overlapping multi-qubit measurements to diagnose channel errors 
and corrects for them. The receiver performs an online decoding circuit to recover the sender's 
original stream of information qubits. 



4.1.2 Quantum Convolutional Code Operation 

Figure 4.1 illustrates the basic operation of a quantum convolutional code. The operation of a 
rate-Zc/rt quantum convolutional code consists of several steps: 

(i) . The protocol begins with the sender encoding a stream of information qubits with an online 

encoding circuit. The sender encodes n^k ancilla qubits and fc information qubits per frame 
[47, 45]. The encoding circuit is "online" if it acts on a few frames at a time. 

(ii) . The sender transmits a set of initial qubits as soon as the first unitary finishes processing 

them and continues to transmit later qubits after the next part of the online encoding circuit 
finishes processing them. 

(iii) . The above basic set Qq and all of its n-qubit shifts act like a parity check matrix for the 

quantum convolutional code. The receiver measures the generators in the stabilizer to de- 
termine an error syndrome. It is important that the generators in Qq have finite weight so 
that the receiver can perform the measurements and produce an error syndrome. It is also 
important that the generators have a block-band form so that the receiver can perform the 
measurements online as the noisy encoded qubits arrive. 

(iv) . The receiver processes the error syndromes with an online classical error estimation algorithm 

such as the Viterbi algorithm [2] or any other decoding algorithm [80] to determine the 
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most likely error for each frame of quantum data. Syndrome-based versions of the Viterbi 
algorithm that are appropriate for quantum coding are available in Refs. [43, 44, 81]. 

(v) . The receiver performs unitary recovery operations that reverse the estimated errors. 

(vi) . He finally processes the encoded qubits with a decoding circuit to recover the original stream 

of information qubits. The qubits decoded from this convolutional procedure should be error 
free and ready for quantum computation at the receiving end. 

In the above description, each step is online — ^the sender and receiver do not have to employ 
the above procedure sequentially. Processing of one step can begin once the previous step has 
finished some amount of processing. All steps above are always active during processing once the 
initial rounds have been completed. 

4.2 The Pauli-to-Binary Isomorphism for Quantum Convo- 
lutional Codes 

This section contains the most important part of this quantum convolutional coding review — the 
isomorphism from the set of Pauli sequences to the module over the ring of binary polynomials 
[44, 45, 1]. We also name it the Pauli-to-binary (P2B) isomorphism (whether we are considering 
the P2B isomorphism for block codes or convolutional codes should be clear from the context). 
The P2B isomorphism is important because it is easier to perform manipulations with vectors of 
binary polynomials than with Pauli sequences. 

We first define the phase- free Pauli group [ll^] on a sequence of qubits. Recall that the delay 
transform D in (4.5) shifts a Pauli sequence to the right by n. Let us assume for now that n = 1. 
Let denote the set of all countably infinite Pauli sequences. The set is equivalent to the 
set of all one-qubit shifts of arbitrary Pauli operators: 




(4.9) 
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We remark that {Ai) — £>* {Ai ® /^^). Wc make this same abuse of notation in what follows. 
We can define the equivalence class [ll^] of phase-free Pauli sequences: 

[n^] = {/3A I Ae n^,/3 e c,|/3| = 1} . (4.10) 

We develop the Pauli-to-binary (P2B) isomorphism between binary polynomials and Pauli 
sequences. The P2B isomorphism is useful for representing the shifting nature of quantum convo- 
lutional codes. Suppose z [D) and x [D) are arbitrary finite-degree and finite-delay polynomials 
in D over Z2: 

z{D) = ^z^D\ x{D) = ^XiD\ G Z2 G Z, (4.11) 

i i 

where del(2; del(x deg{z deg(x {D)) < oo. Suppose 

u {D) = {z {D) , X {D)) G (Z2 {D)f , (4.12) 

where (Z2 (D))^ indicates the Cartesian product Z2 (D) x Z2 (D). Let us employ the following 
shorthand: 

u{D) = [z{D)\x{D)]. (4.13) 
Let N he a map from the binary polynomials to the Pauli sequences, N : (Z2 {D)y n^, where 

N {u (£))) = {Z'^'X'"') . (4.14) 

i 

Let V {D) = [z' {D) \x' {D)] where v (£>) e (Z2 (£>)). The map N induces an isomorphism 

[N]:{Z,{D)f^[U^], (4.15) 

because addition of binary polynomials is equivalent to multiplication of Pauli elements up to a 
global phase: 

[N {u {D) + V {D))] = [N {u {D))] [N {v {D))] . (4.16) 
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The above isomorphism is a powerful way to capture the infiniteness and shifting nature of con- 
volutional codes with finite-degree and finite-delay polynomials over the binary field Z2. 

A quantum convolutional code in general c;onsists of generators with n rjTibits per frame. 
Therefore, we consider the n-qubit extension of the definitions and isomorphism given above. Let 
the delay transform D now shift a Pauli sequence to the right by an arbitrary positive integer n. 
Consider a 2n-dimensional vector u {D) of binary polynomials where u (D) e (Z2 (£>))^". Let us 
write u {D) as follows 

u(^) = [z(^)|x(D)]= [ ^l(D) ••• Zn{D)\x^{D) ■■■ Xn{D)], 

where z (D) ,x(D) e (Z2 (D))". Suppose 

Zi (D) = Y^Zi^jD\ Xi (D) = 

j j 

Define a map N : (Z2 (£1))^" ^ H^: 

N (u (£))) = l\D^ (Z^i.^X^i.J) {1 Z^^.^X^^.i) ...j^j ^j®n-i _ 

j 

N is equivalent to the following map (up to a global phase) 

N (u {D)) = N (ui {D)) (/ ® TV {U2 {D))) ■ ■ ■ ® N K {D))) , 

where 

Ui {D) = [zi {D) \xi (£>)] . (4.17) 

Suppose 

v(D) = [z'(^)|x'(D)], (4.18) 

where v (£») e (Z2 (£>))^". The map N induces an isomorphism [N] : (Z2 (D))^" ^ [n^] for the 
same reasons given in (4.16): 

[N (u {D) + V m] = [N (u {D))] [N (V {D))] . (4.19) 
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The P2B isomorphism [N] is again useful because it allows us to perform binary calculations 
instead of Pauli calculations. 



4.3 The Shifted Symplectic Product 



Wc consider the commutative properties of quantum convolutional codes in this section. We 
develop some mathematics for the important "shifted symplectic product." The shifted symplectic 
product reveals the commutation relations of an arbitrary number of shifts of a set of Pauli 
sequences. 

Recall from Definition 4.1.1 that a commuting set comprising a basic set of Paulis and all their 
shifts specifies a quantum convolutional code. How can we capture the commutation relations of 
a Pauli sequence and all of its shifts? The shifted symplectic product 0, where 



is an elegant way to do so. The shifted symplectic product maps two vectors u {D) = [z {D) \x {D)] 
and V {D) = [z' {D) \x' {D)] to a binary polynomial with finite delay and finite degree: 



The shifted symplectic product is not a proper symplectic product because it fails to be alternating 
[77]. The alternating property requires that 



: (Z2 {D)f X (Z2 {D)f ^ Z2 (£>) 



(4.20) 



{u v) {D) = z (£>) x' (D-^) - X {D) z' {D'^) . 



(4.21) 



(u v) {D) = -{vQu) (£>) 



(4.22) 



but we find instead that the following holds: 



{uQv) (£>) = -{vQu) (D-i) . 



(4.23) 



Every vector u (D) G Z2 (D) is self time-reversal antisymmetric with respect to 0: 



{u u) {D) = -{uQu) (£)-!) \/u {D) G Z2 {Df . 



(4.24) 
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Every binary vector is also self time-reversal symmetric with respect to because addition and 
subtraction are the same over Z2. We employ the addition convention from now on and drop the 
minus signs. The shifted symplectic product is a binary polynomial in D. We write its coefficients 
as follows: 

{u v) {D) = ^ (u v)^ D\ (4.25) 

iGZ 

The coefficient (m w)^ captures the commutation relations of two Pauli sequences for i one-qubit 
shifts of one of the sequences: 

N {u (£»)) {N {v {D})) = D' {N {v (D))) N {u {D)) . (4.26) 

Thus two Pauli sequences A'' {u (D)) and A'' {v {D)) commute for all shifts if and only if the shifted 

symplectic product (u v) (D) vanishes. 

The next example highlights the main features of the shifted symplectic product and further 
emphasizes the relationship between Pauli commutation and orthogonality of the shifted symplec- 
tic product. 

Example 4.3.1. Consider two vectors of binary polynomials: 

u{D)=[ D \ 1 + ], v{D)=[i + D \ 

The P2B isomorphism maps the above polynomials to the following Pauli sequences: 

N {u {D)) = {... \I\X\Z\I\X\I\ ...), N{v {D)) = (• • • \I\Z\Z\I\X\I\ • • • ) • (4.27) 

The vertical bars between every Pauli in the sequence indicate that we are considering one-qubit 
shifts. We determine the commutation relations of the above sequences by inspection. N {u {D)) 
anticommutes with a shift of itself by one or two to the left or right and commutes with all other 
shifts of itself. N {v (£>)) anticommutes with a shift of itself by two or three to the left or right 
and commutes with all other shifts of itself. N{u{D)) anticommutes with N{v{D)) shifted to 
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the left by one or two, with the zero-shifted N {v (D)), and with {v (D)) shifted to the right by 
two or three. The following shifted symplectic products give us the same information; 

(tt u) (D) = D-^ + D-^ + D + D^, {vQ v) {D) = _^ jj-2 ^ jj2 ^ jj3^ 
{vQu){D) = D-^ + D-'^ + l + D + D^. (4.28) 

The nonzero coefficients indicate the commutation relations just as (4.26) claims. 

We can again define a shifted symplectic product for the case of n-qubits per frame. Let 
denote the shifted symplectic product between vectors of binary polynomials: 

: (Z2 (£>))'" X (Z2 (£>))'" ^ Z2 (D) . (4.29) 

It maps 2n-dimensional vectors u {D) and v {D) of binary polynomials to a finite-degree and 
finite-delay binary polynomial 

n 

(u0v)(£>) =;^(ui0t;i)(D), (4.30) 

i=l 

where 

Ui {D) = [zi {D) \xi {D)] , V, {D) = [zl {D) \x'^ {D)] . 
The standard inner product gives an alternative way to define the shifted symplectic product: 

(u v) (D) =z{D)- x' (D-i) - X (D) • z' (D-^) . (4.31) 

Every vector u{D) G %i (D)^" is self time-reversal symmetric with respect to 0: 

(u0u)(Li) = (u0u) (D-i) yu{D) eZ2{Df" . (4.32) 

The shifted symplectic product for vectors of binary polynomials is a binary polynomial in D. We 

write its coefficients as follows: 

(u0v)(^) = ^(u0v). D\ (4.33) 
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The coefficient (u0 v),- captures the commutation relations of two Pauli sequences for i n-qubit 
shifts of one of the sequences: 



N (u (D)) (N (v {D))) = (-1)("®"'- D' (N (v (D))) N (u (D)) . 

The shifted symplectic product between two vectors of binary polynomials vanishes if and only if 

their corresponding Pauli sequences commute. 

Example 4.3.2. We consider the case where the frame size n = 4. Consider the following vectors 
of polynomials: 



(4.34) 



" u{D) ' 




1 + D D 1 D 


10 


_y{D) _ 




10 


1 + D 1 + D 1 D 



The P2B isomorphism maps u {D) and v (D) to the following Pauli sequences: 



N (u (D)) = (• • • \IIII\ZXZI\ZZIZ\IIII\ ■■■), 
N (v (D)) = (• • • \IIII\XYXI\XXIX\IIII\ ■■■). 



(4.35) 



Wc can determine the commutation relations by inspection of the above Pauli sequences. N (u (D)) 
anticommutes with itself shifted by one to the left or right, N(v(£')) anticommutes with itself 
shifted by one to the left or right, and N (u (£))) anticommutes with N (v (D)) shifted by one to 
the left. The following shifted symplectic products confirm the above commutation relations: 

(u © u) (£)) = + £», (v©v)(L») = + {uQv){D) = D'K 

We note two useful properties of the shifted symplectic product 0. Suppose / (D) G Z2 (D) 
with deg (/) > 0. Let us denote scalar polynomial multiplication as follows: 



(/ u) {D) 



f{D)u, {D) 



f{D) Un {D) 



(4.36) 
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The following identities hold. 



((/ u)0v) {D)^f{D) (u0v) (D), 



(4.37) 



(u0(/ v))(D) = /(l?-i)(u©v) {D). 



(4.38) 



We also remark that 



(u v) {D) = (v u) {D) 



iff 



(u0v) (I?) = (u0v) (Z?-i). 



The shifted symplectic product vanishes for all generators in a quantum convolutional code 
because of the commutativity requirement in Definition 4.1.1. The cases where the shifted sym- 
plectic product does not vanish (where the two Pauli sequences anticommute for one or more 
shifts) are important for constructing entanglement-assisted quantum convolutional codes. 

In general, we represent a rate-fc/n quantum convolutional code with an (n — A;) x2n-dimensional 
quantum check matrix H (D) according to the P2B isomorphism. The entries of H (D) are binary 
polynomials where 



and Z (D) and X (D) are both (n — fc) x n-dimensional binary polynomial matrices. The following 
matrix Cl (D) captures the commutation relations of the generators in H (D): 



The reader can verify that the matrix elements [Cl {D)]-j of Cl {D) are the shifted symplectic 
products between the i^^ and j*^ respective rows hi {D) and hj {D) of H (D): 



We call the matrix Cl (D) the shifted symplectic product matrix because it encodes all of the shifted 
symplectic products or the commutation relations of the code. This matrix is equal to the null 



H{D) 



Z{D) X{D) 



fl{D) = Z (D) + X (D) Z'^ . 



[nmii = ihiQhj)iD). 
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matrix for a valid quantum convolutional code because all generators commute with themselves 
and with all n-qubit shifts of themselves and each other. For a general set of generators, f2 (D) is 
not equal to the null matrix and obeys the symmetry: Q (D) = fi^ (i)"^). 

4.4 Row and Column Operations 

We can perform row operations on binary polynomial matrices for quantum convolutional codes. 
A row operation is merely a "mental" operation that has no effect on the states in the codespace 
or on the error-correcting properties of the code. It just changes the rows of the check matrix for 
a code. We have three types of row operations: 

(i) . An elementary row operation multiplies a row times an arbitrary binary polynomial and 

adds the result to another row. This additive invariance holds for any code that admits a 
description within the stabilizer formalism. Additive codes are invariant under multiplication 
of the stabilizer generators in the "Pauli picture" or under row addition in the "binary- 
polynomial picture." 

(ii) . Another type of row operation is to multiply a row by an arbitrary power of D. OUivier 

and Tillich discuss such row operations as "multiplication of a line by £>" and use them to 
find encoding operations for their quantum convolutional codes [44]. Grassl and Rotteler 
use this type of operation to find a subcode of a given quantum convolutional code with an 
equivalent asymptotic rate and equivalent error-correcting properties [45]. 

(iii) . We also employ row operations that multiply a row by an arbitrary polynomial (not nec- 

essarily a power of D). This type of row operation occurs when we have generators with 
infinite weight that we would like to reduce to finite weight so that the receiver can perform 
measurements in an online fashion as qubits arrive from the noisy channel. 

We can encode all of the above types of row operations in a full-rank matrix R (D) with rational 
polynomial entries. Let H' (D) denote the resulting check matrix after performing a set of row 
operations in the matrix R (D) where H' (D) = R (D) H (D). The resulting effect on the shifted 
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symplectic product matrix f2 (D) is to change it to another shifted symplectic product matrix 
n' {D) related to n {D) by 



Row operations do not change the commutation relations of a valid quantum convolutional code 
because its shifted symplectic product matrix is equal to the null matrix. But row operations do 
change the commutation relations of a set of generators whose corresponding shifted symplectic 
product matrix is not equal to the null matrix. This ability to change the commutation relations 

through row operations is crucial for constructing entanglement-assisted quantum convolutional 
codes from an arbitrary set of generators. We use entanglement to resolve any anticommutativity 
in the generators. 

We can also perform column operations on binary polynomial matrices for quantum convolu- 
tional codes. Column operations do change the error-correcting properties of the code and are 
important for realizing a periodic encoding circuit for the code. We have two types of column 
operations: 

(i) . An elementary column operation multiplies one column by an arbitrary binary polynomial 

and adds the result to another column. We implement elementary column operations with 
gates from the shift-invariant Clifford group [47, 45]. 

(ii) . Another column operation is to multiply column i in both the "X" and "Z" matrix by 

where I G Z. We perform this operation by delaying or advancing the processing of qubit i 
by I frames relative to the original frame. 

(iii) . An infinite-depth column operation multiplies one column in the "X" matrix by a rational 

polynomial whose numerator is one and multiplies the corresponding column in the "Z" 
matrix by a corresponding finite polynomial. 

A column operation implemented on the "X" side of the binary polynomial matrix has a 

corresponding effect on the "Z" side of the binary polynomial matrix. This corresponding effect 
is a manifestation of the Hciscnbcrg uncertainty principle because commutation relations remain 
invariant with respect to the action of unitary quantum gates. The shifted symplectic product is 




(4.39) 
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therefore invariant with respect to cohimn operations from the shift-invariant Chfford group. The 
next two sections describe possible column operations for implementing encoding circuits. 

4.5 Finite-Depth Clifford Operations 

One of the main advantages of a quantum convolutional code is that its encoding circuit has a 
periodic form. We can encode a stream of quantum information with the same physical routines 
or devices and therefore reduce encoding and decoding complexity. 

Ollivier and Tillich were the first to discuss encoding circuits for quantum convolutional codes 
[43, 44]. They provided a set of necessary and sufiicient conditions to determine when an encoding 
circuit is noncatastrophic. A noncatastrophic encoding circuit does not propagate uncorrected 
errors infinitely through the decoded information qubit stream. Classical convolutional coding 
theory has a well developed theory of noncatastrophic encoding circuits [80] . 

Grassl and Rotteler later showed that Ollivier and Tillich's conditions for a circuit to be 
noncatastrophic are too restrictive [45, 46, 47]. They found subcodes of quantum convolutional 
codes that admit noncatastrophic encoders. The noncatastrophic encoders are a sequence of 
Clifford circuits with finite depth. They developed a formalism for encapsulating the periodic 
structure of an encoding circuit by decomposing the encoding circuit as a set of elementary column 
operations. Their decoding circuits are exact inverses of their encoding circuits because their 
decoding circuits perform the encoding operations in reverse order. 

Definition 4.5.1. A finite-depth operation transforms every finite-weight stabilizer generator to 
one with finite weight. 

This property is important because we do not want the decoding circuit to propagate uncor- 
rected errors into the information qubit stream [80]. A finite-depth decoding circuit corresponding 
to any stabihzer for a quantum convohitional code exists by the algorithm given in Ref . [45] . 

We review the finite-depth operations in the shift-invariant Clifford group [45, 46, 47]. The 
shift-invariant Clifford group is an extension of the Clifford group operations mentioned in Sec- 
tion 2.2. We describe how finite-depth operations in the shift-invariant Clifford group affect the 
binary polynomial matrix for a code. Each of the following operations acts on every frame of a 
quantum convolutional code. 
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(i) . The sender performs a CNOT from qubit i to qubit j in every frame where qubit j is in a 

frame delayed by k. The effect on the binary polynomial matrix is to multiply column i by 
and add the result to column j in the "X" matrix and to multiply column j by and 
add the result to column i in the "Z" matrix. 

(ii) . A Hadamard on qubit i swaps column i in the "X" matrix with column i in the "Z" matrix. 

(iii) . A phase gate on qubit i adds column i from the "X" matrix to column i in the "Z" matrix. 

(iv) . A controlled-phase gate from qubit i to qubit j in a frame delayed by k multiplies column i in 

the "X" matrix by D'' and adds the result to column j in the "Z" matrix. It also multiplies 
column j in the "X" matrix by D~'^ and adds the result to column i in the "Z" matrix. 

(v) . A controlled-phase gate from qubit i to qubit i in a frame delayed by k multiplies column i 

in the "X" matrix by D'' + D~'^ and adds the result to column i in the "Z" matrix. 

We use finite-depth operations extensively in the next few chapters. Figure 5.2 gives an ex- 
ample of an entanglement-assisted quantum convolutional code that employs several finite-depth 
operations. The circuit encodes a stream of information qubits with the help of ebits shared 
between sender and receiver. 

Multiple CNOT gates can realize an elementary column operation as described at the end of 
the previous section. Suppose the elementary column operation multiplies column i in the "X" 
matrix by / (£>) and adds the result to column j. Polynomial / {D) is a summation of some finite 
set {/i, . . . , In} of powers of D: 

f{D) =£)'i +••• + £>'". 

We can realize / [D] by performing a CNOT gate from qubit i to qubit j in a frame delayed by 
li for each i e {!,..., n}. 

4.6 Infinite-Depth Clifford Operations 

We now introduce an infinite-depth operation to the set of operations in the shift-invariant Clifford 
group available for encoding and decoding quantum convolutional codes. 
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Definition 4.6.1. An infinite- depth operation can transform a finite-weight stabilizer generator 
to one with infinite weight (but does not necessarily do so to every finite-weight generator). 

A decoding circuit with infinite-depth operations on qubits sent over the noisy channel is 
undesirable because it spreads uncorrected errors infinitely into the decoded information qubit 
stream. But an encoding circuit with infinite-depth operations is acceptable if we assume a 
communication paradigm in which the only noisy process is the noisy quantum channel. 

The next chapter shows several examples of circuits that include infinite-depth operations. 
Infinite-depth operations expand the possibilities for quantum convolutional circuits in much the 
same way that incorporating feedback expands the possibilities for classical convolutional circuits. 

We illustrate the details of several infinite-depth operations by first providing some specific 
examples of infinite-depth operations and then show how to realize an arbitrary infinite-depth 
operation. 

We consider both the stabilizer and the logical operators for the information qubits in our 
analysis. Tracking both of these sets of generators is necessary for determining the proper decoding 
circuit when including infinite-depth operations. 

4.6.1 Examples of Infinite-Depth Operations 

Our first example of an infinite-depth operation involves a stream of information qubits and ancilla 
qubits. We divide the stream into frames of three qubits where each frame has two ancilla qubits 
and one information qubit. The following two generators and each of their three-qubit shifts 
stabilize the initial qubit stream: 



III 
III 



Z I I 
I Z I 



III 
III 



(4.40) 



The binary polynomial matrix corresponding to this stabilizer is as follows: 



1 




















1 















(4.41) 
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Wc obtain any Pauli sequence in the stabilizer by multiplying the above rows by a power of D and 
applying the inverse of the P2B isomorphism. The logical operators for the information qubits 
are as follows: 

' ' ' / / / 



III 

III 



I I X 

I I z 



III 

They also admit a description with a binary polynomial matrix: 


















1 








1 












(4.42) 



We refer to the above matrix as the "information-qubit matrix." 



Encoding 

Suppose we would hke to encode the above stream so that the following generators stabilize it: 



III 
III 



XXX 
Z Z I 



XXI 

III 



or equivalently, the following binary polynomial matrix stabilizes it: 




110 



D+1 D+1 1 




(4.43) 



We encode the above stabilizer using a combination of finite-depth operations and an infinite- 
depth operation. We perform a Hadamard on the first qubit in each frame and follow with a CNOT 
from the first qubit to the second and third qubits in each frame. These operations transform the 
matrix in (4.41) to the following matrix 












1 


1 


1 


1 


1 
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or equivalently transform the generators in (4.40) to the foUowiiig generators: 



III 
III 

The information-qubit matrix becomes 



XXX 
Z Z I 



III 
III 


















1 


1 





1 












Wc now perform an infinite-depth operation: a CNOT from the third qubit in one frame to the 
third qubit in a delayed frame and repeat this operation for all following frames. Figure 4.2 shows 
this operation acting on our stream of qubits with three qubits per frame. The effect of this 
operation is to translate the above stabilizer generators as follows: 



III 
III 



XXX 
Z Z I 



I I X 
III 



I I X 

III 



The first generator above and each of its three-qubits shifts is an infinite-weight generator if the 
above sequence of CNOT gates acts on the entire countably-infinite qubit stream. We represent 
the above stabilizer with the binary rational polynomial matrix 




110 



1 1 1/(1 + £>) 




(4.44) 



where 1/(1 + £>) = 1 + D + D'^ + . . . is & repeating fraction. The operation is infinite-depth 
because it translates the original finite- weight stabilizer generator to one with infinite weight. 

It is possible to perform a row operation that multiplies the first row by + 1. This operation 
gives a stabilizer matrix that is equivalent to the desired stabilizer in (4.43). The receiver of 
the encoded qubits measures the finite-weight stabilizer generators in (4.43) to diagnose errors. 
These measurements do not disturb the information qubits because they also stabilize the encoded 
stream. 
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Frame 1 



Frame 2 



Frame 3 



Frame 4 



Frame 5 



Frame 6 



Figure 4.2: An example of an infinite-depth operation. A sequence of CNOT gates acts on the 
third qubit of every frame. This infinite-depth operation effectively multiplies the third column of 
the "X" side of the binary polynomial matrix by the rational polynomial 1/(14-1?) and multiplies 
the third column of the "Z" side of the binary polynomial matrix by 1 -f . 
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The above encoding operations transform the information-qubit matrix as foUows; 





1 1 + 



1/(1 + £>) 




(4.45) 



The infinite-depth operation on the third qubit has an effect on the "Z" or left side of the 
information-qubit matrix as illustrated in the second row of the above matrix. The effect is 
to multiply the third column of the "Z" matrix by / if the operation multiplies the third 

column of the "X" matrix by 1// {D) . This corresponding action on the "Z" side occurs because 
the commutation relations of the Pauli operators remain invariant under quantum gates (due to 
the Heisenberg uncertainty principle), or equivalently, the shifted symplectic product remains in- 
variant under column operations. The original shifted symplectic product for the logical operators 
is one, and it remains as one because f{{D^^)^^)/f{D) = 1. 



Decoding 

We perform finite-depth operations to decode the stream of information qubits. Begin with the 
stabilizer and information-qubit matrix in (4.44) and (4.45) respectively. Perform a CNOT from 
the first qubit to the second qubit. The stabilizer becomes 




10 



1 1/(1 + 1?) 




and the information-qubit matrix does not change. Perform a CNOT from the third qubit to the 
first qubit in the same frame and in a delayed frame. These gates multiply column three in the 
"X" matrix by 1 + and add the result to column one. The gates also multiply column one in 
the "Z" matrix by 1 -|- and add the result to column three. The effect is as follows on both 
the stabilizer 




10 



1/(1 + 1?) 




(4.46) 
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and the information-qubit matrix 


10 

We can multiply the logical operators by any element of the stabilizer and obtain an equivalent 
logical operator [17]. We perform this multiplication in the "binary-polynomial picture" by adding 
the first row of the stabilizer in (4.46) to the first row of (4.47). The information-qubit matrix 
becomes 












1 








1 


















(4.48) 



so that the resulting logical operators act only on the first qubit of every frame. We have success- 
fully decoded the information qubits with finite-depth operations. The information qubits teleport 
coherently [82, 83] from being the third qubit of each frame as in (4.42) to being the first qubit 

of each frame as in (4.48). Wc exploit the above method of encoding with infinite-depth oper- 
ations and decoding with finite-depth operations for the class of entanglement-assisted quantum 
convolutional codes in Section 5.4.1. 

4.6.2 General Infinite-Depth Operations 

We discuss the action of a general infinite-depth operation on two weight-one "X" and "Z" Pauli 
sequences where each frame has one Pauli matrix. Our analysis then determines the effect of an 
infinite-depth operation on an arbitrary stabilizer or information-qubit matrix. The generators in 
the "Pauli picture" are as follows: 

/ 
I 

or as follows in the "binary-polynomial picture" : 






1 


1 






1 1/(1 + £>) 




(4.47) 



X 


I 


z 


I 



(4.49) 
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An infinite-depth 1// (D) operation, where / (D) is an arbitrary polynomial, transforms the above 
matrix to the following one: 

l/f{D) 

A circuit that performs this transformation preserves the shifted symplectic product because 
/ {D~^) ■ 1// = 1. The circuit operates on a few qubits at a time and is shift-invariant — 

the same device or physical routines implement it. 

First perform the long division expansion of binary rational polynomial l/f{D). This ex- 
pansion has a particular repeating pattern with period I. For example, suppose that / {D) = 
1 + D + D^. Its long-division expansion is 1 + D + + + D'' + + + D^^ + . . . and 
exhibits a repeating pattern with period seven. We want a circuit that realizes the following Pauli 
generators 

I T T X X X r X T 

(4.50) 





I 


I 


X 


X 


X 


I 


X 


I 


z 


I 


z 


z 


I 


I 


I 


I 


I 



where the pattern in the X matrices is the same as the repeating polynomial 1// (D) and continues 
infinitely to the right, and the pattern on the Z matrices is the same as that in / (i*^^) and 
terminates at the left. The above Pauli sequence is equivalent to the following binary rational 
polynomial matrix: 



1 + 



D- 



1/(1 + D + D^) 




We now discuss a method that realizes an arbitrary rational polynomial 1// {D) as an infinite- 
depth operation. Our method for encoding the generators in (4.50) from those in (4.49) consists 
of a "sliding- window" technique that determines transformation rules for the circuit. The circuit 
is an additive, shift-invariant filtering operation. It resembles an infinite-impulse response filter 
because the resulting sequence extends infinitely. In general, the number A'' of qubits that the 
encoding unitary operates on is as follows 



AT = deg (/(D)) -del (/(£))) + !, 
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where deg (/(!?)) and del { f{D)) arc the respective highest and lowest powers of polynomial 
f{D). Therefore, our exemplary encoding unitary operates on four qubits at a time. We delay 
the original sequence in (4.49) by three frames. These initial frames are "scratch" frames that 
give the encoding unitary enough "room" to generate the desired Paulis in (4.50). The first set of 
transformation rules is as follows 



I 


I 


I 


X 


I 


I 


I 


X 


I 


I 


I 


z 


z 


I 


z 


z 



(4.51) 



and generates the first four elements of the pattern in (4.50). Now that the encoding unitary 
has acted on the first four frames, we need to shift our eyes to the right by one frame in the 
sequence in (4.50) to determine the next set of rules. So we shift the above outputs by one frame 
to the left (assuming that only identity matrices lie to the right) and determine the next set of 
transformation rules that generate the next elements of the sequence in (4.50): 





X 


I 


I 


I 


X 


z 


z 


I 


I 


z 


z 



X 

I 



Shift the above outputs to the left by one frame to determine the next set of transformation rules: 



/ 
Z 



X 
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I 
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z 


z 


I 



X 

I 



We obtain the rest of the transformation rules by continuing this sliding process, and we stop 
when the pattern in the sequence in (4.50) begins to repeat: 
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X 


I 
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X 


I 


X 


I 
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I 


X 
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The above set of rules determines the encoding unitary and only a few of them are actually 
necessary. We can multiply the rules together to form equivalent rules because the circuit obeys 
additivity (in the "binary-polynomial picture"). The rules become as follows after rearranging 
into a standard form: 
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I 
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I 


X 



A CNOT from qubit one to qubit four and a CNOT from qubit three to qubit four suffice to 
implement this circuit. We repeatedly apply these operations shifting by one frame at a time 
to implement the infinite-depth operation. We could have observed that these gates suffice to 
implement the "Z" transformation in the first set of transformation rules in (4.51). but we wanted 
to show how this method generates the full periodic "X" sequence in (4.50). Figure 4.3 shows how 
the above encoding unitary acts on a stream of quantum information. 

We can determine the encoding unitary for an arbitrary rational polynomial 1// (D) using a 
similar method. Suppose that del (/ {D)) = n and suppose n ^ as in the above case. First delay 
or advance the frames if n > or if n < respectively. Determine the CNOT gates that transform 
the "Z" Pauli sequence 

^ 



to 







These CNOT gates form the encoding circuit that transform both the "X" and "Z" Pauli sequences. 
We perform the encoding unitary, shift by one frame, perform it again, and keep repeating. Our 
method encodes any arbitrary polynomial l/f{D) on the "X" side and / (D^^) on the "Z" side. 
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Frame 1 { 
Frame 2 { 
Frame 3 { 
Frame 4 { 
Frame 5 { 
Frame 6 { 
Frame 7 { 
Frame 8 { 



Figure 4.3: Another example of an infinite-depth operation. An infinite-depth operation acts on 
qubit i in every frame. This particular infinite-depth operation multiplies column i on the "X" 
side of the binary polynomial matrix by 1/ [l + D + D^^ and multiplies column i on the "Z" side 
of the binary polynomial matrix by 1 + D^^ + D^^ . 



We can implement the "time-reversed" polynomial \/ f {D ^) on the "X" side by first delaying 
the frames by m = deg(/(I?)) — del(/(_D)) frames and performing the circuit corresponding to 
1/D'" (/ (D-^)). These operations implement the circuit D™/!)™ (/ (D-^)) ^ 1/f (D-^). 

4.6.3 Infinite-Depth Operations in Practice 

We assume above that each of the infinite-depth operations acts on the entire countably-infinite 
stream of qubits. In practice, each infinite-depth operation acts on a finite number of qubits at a 
time so that the encoding and decoding circuits operate in an "online" manner. Therefore, each 
infinite-depth operation approximates its corresponding rational polynomial. This approximation 
does not pose a barrier to implementation. We can implement each of the above infinite-depth 
operations by padding the initial qubits of the information qubit stream with some "scratch" 
qubits. We first transmit these "scratch" qubits that contain no useful quantum information so 
that the later information qubits enjoy the full protection of the code. These scratch qubits do 
not affect the asymptotic rate of the code and merely serve as a convenience for implementing the 
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infinite-depth operations. From now on, we adhere to describing infinite-depth operations with 
binary rational polynomials because it is more convenient to do so mathematically. 



4.7 Closing Remarks 

This chapter reviewed in detail the method of quantum convolutional coding. We discussed how 
to represent a quantum convolutional code with a matrix of binary polynomials. The shifted 
symplectic product is crucial in several of the forthcoming chapters because it determines the 
commutation relations for an arbitrary set of convolutional generators. In general, a set of genera- 
tors do not form a commuting set and we can use entanglement to resolve their anticommutativity. 
In this chapter, we also developed many of the operations that manipulate quantum convolutional 
codes. We repeatedly use these operations in the next few chapters. 
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Chapter 5 



Entanglement- Assisted Quantum Convolutional 
Coding: The CSS Case 



O Calderbank, good Shor, and bright Steane, 
Methinks your construction is keen, 
But it makes me tense, 
this Pauli sequence, 
I just want our quantum bits clean. 

In this chapter, wc develop a theory of entanglement-assisted quantum convolutional coding for a 
broad class of codes. Our major result is that we can produce an entanglement-assisted quantum 
convolutional code from two arbitrary classical binary convolutional codes. The resulting quantum 
convolutional codes admit a Calderbank-Shor-Steane (CSS) structure [21. 22. 32]. The rates and 
error-correcting properties of the two binary classical convolutional codes directly determine the 
corresponding properties of the entanglement-assisted quantum convolutional code. 

Our CSS entanglement-assisted quantum convolutional codes divide into two classes based on 
certain properties of the classical codes from which we produce them. These properties of the 
classical codes determine the structure of the encoding and decoding circuit for the code, and the 
structure of the encoding and decoding circuit in turn determines the class of the entanglement- 
assisted quantum convolutional code. 

(i). Codes in the first class admit both a finite-depth encoding and decoding circuit. 
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(ii). Codes in the second class have an encoding circuit that employs both finite-depth and 
infinite-depth operations. Their decoding circuits have finite-depth operations only. 

We structure this chapter as follows. We outline the operation of an entanglement-assisted 
quantum convolutional code and present our main theorem in Section 5.1. This theorem shows 
how to produce a CSS entanglement-assisted quantum convolutional code from two arbitrary 
classical binary convolutional codes. The theorem gives the rate and error-correcting properties 
of a CSS entanglement-assisted quantum convolutional code as a function of the parameters of 
the classical convohitional codes. Section 5.3 completes the proof of the theorem for our first class 
of entanglement-assisted quantum convolutional codes. Section 5.4 completes the proof of our 
theorem for the second class of entangiement-assisted quantum convolutional codes. We discuss 
the implications of the assumptions for the different classes of entanglement-assisted quantum 
convolutional codes while developing the constructions. Our theory produces high-performance 
quantum convolutional codes by importing high-performance classical convolutional codes. 

5.1 Operation of an Entanglement- Assisted Quantum Con- 
volutional Code 

An entanglement-assisted quantum convolutional code operates similarly to a standard quantum 

convolutional code. The main difference is that the sender and receiver share entanglement in the 
form of ebits before quantum communication begins. An [[n, k\ c]] entangiement-assisted quan- 
tum convolutional code encodes k information qubits per frame with the help of c ebits and 
n — k — c ancilla qubits per frame. Figure 5.1 highlights the main features of the operation of an 
entanglement-assisted quantum convolutional code. An entanglement-assisted quantum convolu- 
tional code operates as follows: 

(i). The sender encodes a stream of quantum information using both additional ancillas and 

ebits. The sender performs the encoding operations on her qubits only (i.e., not including 
the halves of the ebits in possession of the receiver). The encoding operations have a periodic 
structure so that the same operations act on qubits in different frames and give the code a 
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Frame 1 



Frame 2 



Frame 3 



Encoding 



0)-" 

i'l) 



_ U 



Noise 



# 



# 



Measurement and Recovery 



Decoding 



Figure 5.1: (Color online) An entanglement-assisted quantum convolutional code operates on a 
stream of qubits partitioned into a countable number of frames. The sender encodes the frames of 
information qubits, ancilla qubits, and half of shared ebits with a repeated, overlapping encoding 
circuit U. The noisy channel affects the sender's encoded qubits but does not affect the receiver's 
half of the shared ebits. The receiver performs overlapping measurements on both the encoded 
qubits and his half of the shared ebits. These measurements produce an error syndrome which 
the receiver can process to determine the most likely error. The receiver reverses the errors on 
the noisy qubits from the sender. The final decoding circuit operates on all qubits in a frame and 
recovers the original stream of information qubits. 
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memory structure. The sender can perform these encoding operations in an online manner 
as she places more qubits in the unencoded qubit stream. 

(ii) . The sender transmits her encoded qubits over the noisy quantum communication channel. 

The noisy channel does not affect the receiver's half of the shared ebits. 

(iii) . The receiver combines the received noisy qubits with his half of the ebits and performs 

measurements to diagnose errors that may occur. These measurements may overlap on 
some of the same qubits. 

(iv) . The receiver then diagnoses errors using a classical error estimation algorithm such as Viterbi 

error estimation [2] and reverses the estimates of the errors. 

(v) . The receiver finally performs an online decoding circuit that outputs the original information 

qubit stream. 

5.2 The CSS Construction 

Our main theorem below allows us to import two arbitrary classical convolutional codes for use 
as a CSS entanglement-assisted quantum convolutional code. Grassl and Rotteler were the first 
to construct CSS quantum convolutional codes from two classical binary convolutional codes that 
satisfy an orthogonality constraint — the polynomial parity check matrices Hi (D) and H2 {D) of 
the two classical codes are orthogonal with respect to the shifted symplectic product: 

Hi{D)H^{D-^)=0. (5.1) 

The resulting symplectic code has a self-orthogonal parity-check matrix when we join them to- 
gether using the CSS construction. Our theorem generalizes the work of Grassl and Rotteler 
because we can import two arbitrary classical binary convolutional codes — ^the codes do not nec- 
essarily have to obey the self-orthogonality constraint. 

The theorem gives a direct way to compute the amount of entanglement that the code requires. 
The number of ebits required is equal to the rank of a particular matrix derived from the check 
matrices of the two classical codes. It generalizes the earlier theorems discussed in Section 3.6 
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that determine the amount of entanglement required for an entanglement-assisted quantum block 
code. 

Theorem 5.2.1 also provides a formula to compute the performance parameters of the entanglement- 
assisted quantum convolutional code from the performance parameters of the two classical codes. 
This formula ensures that high-rate classical convolutional codes produce high-rate entanglement- 
assisted quantum convolutional codes. Our constructions also ensure high performance for the 
"trade-off' and "catalytic" rates by minimizing the number of ebits that the codes require. 

We begin the proof of the theorem in this section and complete it in different ways for each 
of our two classes of entanglement-assisted quantum convolutional codes in Sections 5.3 and 5.4. 
The proofs detail how to encode a stream of information qubits, ancilla qubits, and shared ebits 
into a code that has the CSS structure. 

Theorem 5.2.1. Let Hi [D] and H2 (D) be the respective check matrices corresponding to non- 
catastrophic, delay-free encoders for classical binary convolutional codes Ci and C2. Suppose that 
classical code Ci encodes k-i information bits with n bits per frame where i = 1.2. The respective 
dimensions of Hi (D) and H2 {D) are thus {n — ki) x n and (n — ^2) x n. Then the resulting 
entanglement- assisted quantum convolutional code encodes ki -\- k2 — n -\- c information qubits per 
frame and is an [[n, fci + ^2 — n + c; c]] entanglement- assisted quantum convolutional code. The 
code requires c ebits per frame where c is equal to the rank of Hi (D) H2 (r>~^). 

Let us begin the proof of the above theorem by constructing an entanglement-assisted quantum 
convolutional code. Consider the following quantum check matrix in CSS form; 



HiiD) 








H2 (D) _ 



(5.2) 



We label the above matrix as a "quantum check matrix" for now because it does not necessar- 
ily correspond to a commuting stabilizer. The quantum check matrix corresponds to a set of 
Pauli sequences whose error-correcting properties are desirable. The Pauli sequences inherit these 
desirable properties from the classical codes Ci and C2. 
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The following k;iniiia begins the proof of the above theorem. It details an initial decomposition 
of the above quantum check matrix for each of our two classes of entanglement-assisted quantum 
convolutional codes. 

Lemma 5.2.1. Elementary row and column operations relate the quantum check matrix in (5.2) 
to the following matrix 

E{D) F{D) 
0/0 

where E (D) is dimension (n — ki) x (n — k2), F (D) is (n — ki) x k2, the identity matrix is 
(n — ^2) X (n — k2), and the null matrix on the right is {n — k2) x ^2. We give a definition of 
E [D) and F [D) in the following proof. 



Proof The Smith form [80] of Hi {D) for each j = 1, 2 is 



Hi{D)=Ai{D) I Bi{D) 



(5.3) 



where (D) is (n — fc,) x (n — fcj), the matrix in brackets is (n — fcj) x n, and Bj {D) is nxn. The 
Smith decomposition is somewhat analogous to the SVD decomposition of a matrix over the reals. 
The binary polynomials along the diagonal of the center matrix of the Smith decomposition are 
the invariant factors. The matrices Ai (D) and (D) are a product of a sequence of elementary 
row and column operations respectively [80] . Let Bia {D) be the first n — ki rows of Bi (D) and 
let Bib {D) be the last ki rows of Bi (D): 



Bi{D) 



Bia {D) 

Bib {D) 



The (n — ki) x (n — ki) identity matrix in brackets in (5.3) indicates that the invariant factors of 
Hi (D) for each i = 1,2 are all equal to one [80]. The invariant factors are all unity for both check 
matrices because the check matrices correspond to noncatastrophic, delay-free encoders [80] . 

Premultiplying Hi (D) by Ar^ (£>) gives a check matrix H^(D) for each i = 1,2. Matrix 
H[ (D) is a check matrix for code C, with equivalent error-correcting properties as Hi (D) because 
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row operations relate the two matrices. This new check matrix (D) is equal to the first n — ki 
rows of matrix Bi {D): 

Hl{D) = Bi, {D). 

The invariant factors of Hi (D) H2 {D~^) are equivalent to those of H[ (D) H!^ {^~^) because 
they are related by elementary row and column operations [80] : 

Hi {D) Hi (D-i) = Ai {D) H[ (D) H'^ {D-') [D-') . (5.4) 

We now decompose the above quantum check matrix into a basic form using elementary row 
and column operations. The row operations have no effect on the error-correcting properties of 
the code, and the column operations correspond to elements of an encoding circuit. We later show 
how to incorporate ebits so that the quantum check matrix forms a valid commuting stabilizer. 

Perform the row operations in matrices A^^ (D) for both check matrices Hi (D). The quantum 
check matrix becomes 











B2a (D) _ 



(5.5) 



The error-correcting properties of the above generators are equivalent to those of the generators 
in (5.2) because row operations relate the two sets of generators. The matrix B2 (D) corresponds 
to a sequence {-B2,i of elementary column operations: 

B2 {D) = B2,i (£>)••• B2,i {D) = H B2,i {D) . 

i=l 

The inverse matrix (D) is therefore equal to the above sequence of operations in reverse order: 

1 

B^^ (£») = B2,i (£>)••• B2,i {D) = H B2,i (£>) . 

i=l 

Perform the elementary column operations in B^^ (D) with CNOT and SWAP gates [45]. The 
effect of each elementary column operation i?2,i (D) is to postmultiply the "X" matrix by i32,i {D) 



75 



and to postmultiply the "Z" matrix by Sj^ (_D ^) . Therefore the effect of all elementary operations 
is to postmultiply the "Z" matrix by because 



The quantum check matrix in (5.5) becomes 



Bi„ {D)Bi{D-^) 






1 



Let E [D) be equal to the first n — ki rows and n — k2 columns of the "Z" matrix: 



and let F (D) be equal to the first n — k\ rows and last k2 columns of the "Z" matrix: 



F{D)=B,,^{D)Bl, (D-i). 



The quantum check matrix in (5.6) is then equivalent to the following matrix 



E{D) F{D) 





/ 



where each matrix above has the dimensions stated in the theorem above. 



(5.6) 



(5.7) 



□ 



The above operations end the initial set of operations that each of our two classes of entanglement- 
assisted quantum convolutional codes employs. We outline the remaining operations for each class 
of codes in what follows. 
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5.3 Codes with Finite-Depth Encoding and Decoding Cir- 
cuits 



This section details entanglement-assisted quantum convolutional codes in our first class. Codes 
in the first class admit an encoding and decoding circuit that employ finite-depth operations only. 
The check matrices for codes in this class have a property that allows this type of encoding and 
decoding. The following lemma gives the details of this property, and the proof outlines how to 
encode and decode this class of entanglement-assisted quantum convolutional codes. 

Lemma 5.3.1. Suppose the Smith form of Hi [D) H2 [D^^) is 



Hi{D)Hj{D-')=A{D) 



T{D) 




B{D), 



where A (D) is an {n — ki) x (n — ki) matrix, B (D) is an {n — ^2) x (n — ^2) matrix, T (£)) is a 
diagonal c x c matrix whose entries are powers of D, and the matrix in brackets has dimension 
(n — fci) X (n — ^2). Then the resulting entanglement- assisted quantum convolutional code has both 
a finite- depth encoding and decoding circuit. 

Proof. We begin the proof of this lemma by continuing where the proof of Lemma 5.2.1 ends. The 
crucial assumption for the above lemma is that the invariant factors of Hi (D) if J {D~^) are all 
powers of D. The Smith form of E (D) in (5.7) therefore becomes 



Ai'iD)A{D) 



T{D) 




B{D)A^'{D), 



by employing the hypothesis of Lemma 5.3.1 and (5.4). The rank of both Hi [D) H2 {D~^) and 
E (D) is equal to c. 

Perform the inverse of the row operations in A^^ (D) A (D) on the first n — ki rows of the quan- 
tum check matrix in (5.7). Perform the inverse of the column operations in matrix B (D) A^^ [D) 
on the first n — k2 columns of the quantum check matrix in (5.7). We execute these column oper- 
ations with Hadamard, CNOT, and SWAP gates. These column operations have a corresponding 
effect on columns in the "X" matrix, but we can exploit the identity matrix in the last n—k2 rows 
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of the "X" matrix to counteract this effect. We perform row operations on the last n — k2 rows 

of the matrix that act as the inverse of the column operations, and therefore the quantum check 
matrix in (5.7) becomes 

^ T{D) Fi{D) 

F2{D) 
7 
/0 

where Fi (D) and F2 {D) are the first c and n — k\ — c respective rows of (D) Ai (D) F (D). 
We perform Hadamard and CNOT gates to clear the entries in Fi (D) in the "Z" matrix above. 
The quantum check matrix becomes 



r (D) 

F2{D) 


0/0 







/ 





(5., 



The Smith form of F2 (D) is 



F2 (£)) = (£») TpiD) BpiD) 



where Tf{D) is a diagonal matrix whose entries are powers of D, Ap {D) is (n — fci — c) x 
(n — fci — c), and Bp {D) is k2 x k2. The Smith form of F2 {D) takes this particular form be- 
cause the original check matrix H2 {D) is noncatastrophic and column operations with Laurent 
polynomials change the invariant factors only up to powers of D. 

Perform row operations corresponding to Ap^ (D) on the second set of n — fci — c rows with 
F2 {D) in (5.8). Perform column operations corresponding to Bp^ (D) on columns n — ^2 + 1, . . . , n 
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with Hadamard, CNOT, and SWAP gates. The resulting quantum check matrix has the foUowing 
form; 

'' T{D) 0000 
r^^p) 
0/000 
/ 00000 



(5.9) 



We have now completed the decomposition of the original quantum check matrix in (5.2) for 
this class of entanglement-assisted quantum convolutional codes. It is not possible to perform 
row or column operations to decompose the above matrix any further. The problem with the 
above quantum check matrix is that it does not form a valid quantum convolutional code. The 
first set of rows with matrix F (D) are not orthogonal under the shifted symplectic product to 
the third set of rows with the identity matrix on the "X" side. Equivalently, the set of Pauli 
sequences corresponding to the above quantum check matrix do not form a commuting stabilizer. 
We can use entanglement shared between sender and receiver to solve this problem. Entanglement 
adds columns to the above quantum check matrix to resolve the issue. The additional columns 
correspond to qubits on the receiver's side. We next show in detail how to incorporate ancilla 
qubits, cbits, and information qubits to obtain a valid stabilizer code. The result is that we can 
exploit the error-correcting properties of the original code to protect the sender's qubits. 

Consider the following check matrix corresponding to a commuting stabilizer; 



/ / 

/0 



0/00 







/ / 





(5.10) 



where the identity matrices in the first and third sets of rows each have dimension c x c, the identity 
matrix in the second set of rows has dimension (n — fci — c) x (n — fci — c), and the identity matrix 
in the fourth set of rows has dimension (n — ^2 — c) x [n — — c)- The first and third sets of c 
rows stabilize a set of c ebits shared between Alice and Bob. Bob possesses the "left" c qubits 
and Alice possesses the "right" n qubits. The second and fourth sets of rows stabilize a set of 
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2 {n — c) — ki — k2 ancilla qubits that Alice possesses. The stabihzer therefore stabilizes a set of c 
ebits, 2 (n ~ c) — ki — k2 ancilla qubits, and ki + k2 — n + c information qubits. 

Observe that the last n columns of the "Z" and "X" matrices in the above stabilizer are 
similar in their layout to the entries in (5.9). We can delay the rows of the above stabilizer by 
an arbitrary amount to obtain the desired stabilizer. So the above stabilizer is a subcode of the 
following stabilizer in the sense of Ref. [45]: 



T{D) T{D) 

ri.(D) 



J 







/ / 





The stabilizer in (5.10) has equivalent error-correcting properties to and the same asymptotic rate 
as the above desired stabilizer. The above stabilizer matrix is an augmented version of the quantum 
check matrix in (5.9) that includes entanglement. The sender performs all of the encoding column 
operations detailed in the proofs of this lemma and Lemma 5.2.1 in reverse order. The result of 
these operations is an [[n, ki + k2 — n + c; c]] entanglement-assisted quantum convolutional code 
with the same error-correcting properties as the quantum check matrix in (5.2). The receiver 
decodes the original information-qubit stream by performing the column operations in the order 
presented. The information qubits appear as the last fci + ^2 — n + c in each frame of the stream 
(corresponding to the ki + k2 — n+c columns of zeros in both the "Z" and "X" matrices above). □ 

Example 5.3.1. Consider a classical convolutional code with the following check matrix: 



HiD) 



1 + 1 + D + D'^ 



We can use H (D) in an entanglement-assisted quantum convolutional code to correct for both 
bit-flip errors and phase-flip errors. We form the following quantum check matrix: 



1 + l + D + D"^ 





1 £)2 1 + D + D^ 



(5.11) 
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Frame 1 



Frame 2 



Frame 3 



Frame 4 



Frame 5 



Frame 6 



1^2)^ - 



{ 
{ 
{ 
{ 



l^sf - 

\<^+)BA ^ 
1^4)^ - 



4 — 4 — [h]- 

4 — [h] — i- 
-[Hh 



Figure 5.2: (Color online) The finite-depth encoding circuit for the entanglement-assisted quantum 
convolutional code in Example 5.3.1. The above operations in reverse order give a valid decoding 
circuit. 



This code falls in the first class of entanglement-assisted quantum convolutional codes because 
H{D)H^ (D-i) =1. 

We do not show the decomposition of the above check matrix as outlined in Lemma 5.3.1, but 
instead show how to encode it starting from a stream of information qubits and ebits. Each frame 
has one ebit and one information qubit. Let us begin with a polynomial matrix that stabilizes the 
unencoded state: 



1 


1 























1 


1 






Alice possesses the two qubits on the "right" and Bob possesses the qubit on the "left." We label 
the middle qubit as "qubit one" and the rightmost qubit as "qubit two." Alice performs a CNOT 



81 



from qubit one to qubit two in a delayed frame and a CNOT from qubit one to qubit two in a 
frame delayed by two. The stabilizer becomes 



110 






1 1 D + D^ 



Alice performs Hadamard gates on both of her qubits. The stabilizer becomes 



10 
1 £) + £)2 



10 
10 



Alice performs a CNOT from qubit one to qubit two in a delayed frame. The stabilizer becomes 



10 
D D + D"^ 



1 D 
10 



Alice performs a CNOT from qubit two to qubit one in a delayed frame. The stabilizer becomes 



10 

D l + D + D"^ 



1 + £)2 D 
10 



Alice performs a CNOT from qubit one to qubit two. The stabilizer becomes 



10 

l + £>2 i + D + D^ 



l + £)2 l + £) + £)2 

10 



A row operation that switches the first row with the second row gives the following stabilizer: 



l + £>2 i + D + D^ 
10 



10 

l + £)2 l + £) + £)2 



The entries on Alice's side of the above stabilizer have equivalent error-correcting properties to 
the quantum check matrix in (5.11). Figure 5.2 illustrates how the above operations encode a 
stream of ebits and information qubits for our example. 
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Discussion 

Codes in the first class are more useful in practice than those in the second because their encoding 
and decoding circuits are finite depth. An uncorrected error propagates only to a finite number of 
information qubits in the decoded qubit stream. Codes in the first class therefore do not require 
any assumptions about noiseless encoding or decoding. 

The assumption about the invariant factors in the Smith form of Hi (D) {^~^) holds only 
for some classical check matrices. Only a subclass of classical codes satisfy this assumption, but it 
still expands the set of available quantum codes beyond those whose check matrices H\ (D) and 
H2 (-D) are orthogonal. We need further techniques to handle the classical codes for which this 
assumption does not hold. The following sections provide these further techniques to handle a 
larger class of entanglement-assisted quantum convolutional codes. 

5.4 Codes with Infinite-Depth Encoding and Finite-Depth 
Decoding Circuits 

This section details codes whose encoding circuits have both infinite-depth and finite-depth op- 
erations. We therefore assume that encoding is noiseless to eliminate the possibility of encoding 
errors spreading infinitely into the encoded qubit stream. We later briefly discuss the effects of 
relaxing this assumption in a realistic system. Their decoding circuits require finite-depth opera- 
tions only. Some of the decoding circuits are not the exact inverse of their corresponding encoding 
circuits, but the decoding circuits invert the effect of the encoding circuits because they produce 
the original stream of information qubits at their output. 

Just as with the previous class, this class of codes is determined by the properties of their 
corresponding classical check matrices, as described in the following lemma. 

Lemma 5.4.1. Suppose the Smith form of E (D) does not admit the form from Lemma 5.3.1. 
Then the entanglement-assisted quantum convolutional code has an encoding circuit with both 
infinite- depth and finite-depth operations. Its decoding circuit has finite-depth operations. 
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Proof. Wc perform all of the operations from Lemma 5.2.1. The Smith form of E (D) is in general 
as follows 

Ti {D) 
Ae{D) T^iD) BsiD). 


where Ae {D) is (n — ki) x (n — ki), Fi (D) is an s x s diagonal matrix whose entries are powers 
of D, r2 (-D) is a (c — s) X (c — s) diagonal matrix whose entries are arbitrary polynomials, and 
Be {D) is (n — k2) x (n — ^2). Perform the row operations in A'^^ (D) and the column operations 
in B^^ (D) on the quantum check matrix in (5.7). Counteract the effect of the column operations 
on the identity matrix in the "X" matrix by performing row operations. The quantum check 
matrix in (5.7) becomes 



ri(D) Fi(£>) 

r2(£>) F2(£>) 

F3{D) 











/ 



where Fi (D), F2 {D), and F3 (D) are the respective s, c—s, and n—ki — c rows of A^^ (D) F (D). 
The Smith form of F3 (D) is as follows 



Fs{D)=Af,{D) TeAD) BeAD) 



where Aes {D) is (n — fci — c) x (n — fci — c), Te^ {D) is an (n — fci — c) x {n — ki — c) diagonal 
matrix whose entries are powers of D, and Be^ {D) is k2 x k2- The entries of Fe^ (-D) are powers 
of D because the original check matrix H2 {D) is noncatastrophic and column and row operations 
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with Laurent polynomials change the invariant factors only by a power of D. Perform the row 
operations in [D) and the column operations in Bp^ {D). The quantum check matrix becomes 



T,{D) F,JD) F,,iD) 

T2{D) F^AD) F'^,{D) 

TfAD) 











/ 



where F^^ (D), F[f^ (D), (D), F^^ (D) are the matrices resulting from the column operations 
in Bp^ {D). Perform row operations from the entries in T p.^ [D) to the rows above it to clear the 
entries in F^^ (D) and (D). Use Hadamard and CNOT gates to clear the entries in F-^f^ (D). 
The quantum check matrix becomes 



Pip) 

T^iD) F^.p) 

TfAD) 











/ 



We can reduce F21J (D) to a lower triangular form with an algorithm consisting of column operations 
only. The algorithm operates on the last k2+ki—n+c columns. It is similar to the Smith algorithm 

but does not involve row operations. Consider the first row of (D). Perform column operations 
between the different elements of the row to reduce it to one non-zero entry. Swap this non-zero 
entry to the leftmost position. Perform the same algorithm on elements 2, . . . , A:2 + ^ + c of 
the second row. Continue on for all rows of (-^^) to reduce it to a matrix of the following form 



c-s ki + k2-n+s 



L{D) 
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where L (D) is a lower triangular matrix. The above quantum check matrix becomes 

ri(£») 

r2(£') L (LI) 
Tp^iD) 

00/0_ 

We have completed decomposition of the first set of s rows with Fi (£>), the third set of n — fci — c 
rows with Tp^ (D), and rows n — ki + 1, . . . ,n — ki + s with the identity matrix on the "X" side. 

We now consider an algorithm with infinite-depth operations to encode the following submatrix 
of the above quantum check matrix: 

T2{D) L{D) 


We begin with a set of c — s ebits and c — s information qubits. The following matrix stabilizes 
the ebits 



/ 


/ 























/ 


/ 






and the following matrix represents the information qubits 









I 


























I 



where all matrices have dimension (c — s) x (c — s) and Bob possesses the c — s qubits on the 
"left" and Alice possesses the 2 (c — s) qubits on the "right." We track both the stabilizer and 
the information qubits as they progress through some encoding operations. Alice performs CNOT 
and Hadamard gates on her 2 (c — s) qubits. These gates multiply the middle c — s columns of 
the "Z" matrix by L (D) and add the result to the last c — s columns and multiply the last c — s 
columns of the "X" matrix by (^~^) and add the result to the middle c — s columns. The 
stabilizer becomes 

I I L{D) 





/ 



(5.12) 




7/0 
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and the information-qubit matrix becomes 



0/ 





(£»-!) / 



Alice performs infinite-depth operations on her first c — s qubits corresponding to the rational 
polynomials . . ., 72^c_s {D~^) in r2 ^ (Z)~^). The stabilizer matrix becomes 



I T2{D) L{D) 




and the information-qubit matrix becomes 





1 r^'{D-^) 



0/ 






L'^ {D-^)T^^ (D-^) I 



Alice's part of the above stabilizer matrix is equivalent to the quantum check matrix in (5.12) 
by row operations (premultiplying the second set of rows in the stabilizer by r2 {D).) Bob can 
therefore make stabilizer measurements that have finite weight and that are equivalent to the 
desired stabilizer. 

We now describe a method to decode the above encoded stabilizer and information-qubit 
matrix so that the information qubits appear at the output of the decoding circuit. Bob performs 
Hadamard gates on his first and third sets of c — s qubits, performs CNOT gates from the first 
set of qubits to the third set of qubits corresponding to the entries in L{D), and performs the 
Hadamard gates again. The stabilizer becomes 



/ T2{D) 




and the information-qubit matrix becomes 











(5.13) 



0/ 
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Bob finishes decoding at this point because we can equivalently express the information-qubit 
matrix as follows 









I 


























I 



by multiplying the last c — s rows of the stabilizer by and adding to the last c — s rows 

of the information-qubit matrix. 

The overall procedure for encoding is to begin with a set of c ebits, 2 (n — c) — fci — A;2 ancilla 
qubits, and k\ + k2 — n+c information qubits. Alice performs the infinite-depth operations detailed 
in the paragraph with (5.12) for c— s of the ebits. Alice then performs the finite-depth operations 
detailed in the proofs of this lemma and Lemma 5.2.1 in reverse order. The resulting stabilizer 
has equivalent error-correcting properties to the quantum check matrix in (5.2). 

The receiver decodes by first performing all of the finite-depth operations in the encoding 
circuit in reverse order. The receiver then decodes the infinite-depth operations by the procedure 
listed in the paragraph with (5.13) so that the original k\ + k2 — n + c information qubits per 
frame are available for processing at the receiving end. □ 



5.4.1 Special Case with Coherent Teleportation Decoding 

We now detail a special case of the above codes in this final section. These codes are interesting 
because the information qubits teleport coherently to other physical qubits when encoding and 
decoding is complete. 

Lemma 5.4.2. Suppose that the Smith form of F [D) in (5.7) is 

F {D) = Af (D) [tf {D) o] Bp {D) , 

where Ap {D) is {n — ki) x{n — ki), Tf (D) is an (n — fci) x{n — ki) diagonal matrix whose entries 
are powers of D, and Bf (D) is k^y-k^- Then the resulting entanglement-assisted code admits an 
encoding circuit with both infinite- depth and finite-depth operations and admits a decoding circuit 
with finite- depth operations only. The information qubits also teleport coherently to other physical 
qubits for this special case of codes. 
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Proof. Wc perform all the operations in Lemma 5.2.1 to obtain the quantum check matrix in 
(5.7). Then perform the row operations in Ap^ [D) and the column operations in Bp^ (D). The 
quantum check matrix becomes 



E' (D) Tf (D) 





7 



where E' (D) = Ap'^ {D) E (£>). The Smith form of E' (D) is 



E' {D) = Ae' {D) 



Ti (D) 

r2 (£>) 





Be' (D) 



where Ae' (-D) is (n — ki) x (n — ki), Fi (D) is an s x s diagonal matrix whose entries are powers 
of D, r2 (-D) is a (c — s) X (c — s) diagonal matrix whose entries are arbitrary polynomials, and 
Be' (-D) is (n - ^2) x (n - ^2). 

Now perform the row operations in Apj (D) and the column operations in Bp} (D). It is 
possible to counteract the effect of the row operations on Fir (D) by performing column operations, 
and it is possible to counteract the effect of the column operations on the identity matrix in the 
"X" matrix by performing row operations. The quantum check matrix becomes 



Fi(£>) T[{D) 

T2{D) T'^iD) 

T'siD) 











7 



where r[{D), r'2{D), and T'siD) represent the respective sxs, {c — s)x{c — s), and (n — fci — c) x 
(n — fci — c) diagonal matrices resulting from counteracting the effect of row operations Apj (D) 



89 



on Fj- (D). Wc use Hadamard and CNOT gates to clear the entries in T[ (D). The quantum check 
matrix becomes 

Ti{D) 00 0000 

T2{D) r'2(£>) 

r'3(D) 

00 0700 

The first s rows with Ti{D) and rows n — ki — c + 1, . . . ,n — ki — c + s with the identity matrix 
on the "X" side stabiUze a set of s ebits. The n — ki — c rows with T'^{D) and the n — k2 — c 
rows with identity in the "X" matrix stabiUze a set of 2 (n — c) — fci — k2 ancilla qubits (up to 
Hadamard gates). The s and k2 — n + ki columns with zeros in both the "Z" and "X" matrices 
correspond to information qubits. The decomposition of these rows is now complete. 

We need to finish processing the c — s rows with r2 (D) and Ff, (D) as entries and the c — s 
rows of the identity in the "X" matrix. We construct a submatrix of the above quantum check 
matrix; 



We describe a procedure to encode the above entries with c — s ebits and c — s information 
qubits using infinite-depth operations. Consider the following stabilizer matrix 



/ 


/ 























/ 


/ 






where all identity and null matrices are {c — s) x {c — s). The above matrix stabilizes a set of c — s 
ebits and c— s information qubits. Bob's half of the ebits are the c—s columns on the left in both 
the "Z" and "X" matrices and Alice's half are the next c—s columns. We also track the logical 
operators for the information cjubits to verify that the circuit encodes and decodes properly. The 
information-qubit matrix is as follows 









I 


























I 



(5.16) 




/ 



(5.14) 
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where all matrices are again (c — s) x (c — s). Alice perforins Hadaniard gates on her first c — s 
qubits and then performs CNOT gates from her first c — s qubits to her last c — s qubits to 
transform (5.15) to the following stabilizer 



7 

0/0 



/ T'^{D) 
10 



The information-qubit matrix in (5.16) becomes 



r'2 (£>-!) / 






0/ 



Alice then performs infinite-depth operations on her last c — s qubits. These infinite-depth opera- 
tions correspond to the elements of ^ {D). She finally performs Hadamard gates on her 2 (c — s) 
qubits. The stabilizer becomes 








/ / 



(5.17) 



and the information-qubit matrix becomes 




r^^{D) 



(/?-!) T^iD-^) 





(5.18) 



The stabilizer in (5.17) is equivalent to the following stabilizer by row operations (premultiplying 
the first c — s rows by r2 {D)): 



T2{D) T2{D) T^{D) 





/ / 



(5.19) 



The measurements that Bob performs have finite weight because the row operations are multipli- 
cations of the rows by the arbitrary polynomials in r2 (D) . Alice thus encodes a code equivalent 
to the desired quantum check matrix in (5.14) using c — s ebits and c — s information qubits. 
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Wc now discuss decoding the stabilizer in (5.17) and information qubits. Bob performs CNOTs 
from the first c — s qubits to the next c — s qubits. The stabihzer becomes 



/ T^'iD)r^{D) 





7 



(5.20) 



and the information-qubit matrix does not change. Bob uses Hadamard and finite-depth CNOT 
gates to multiply the last c — s columns in the "Z" matrix by r2 {D~^) T2 (D) and add the result 
to the middle c — s columns. It is possible to use finite-depth operations because the entries of 
r'2 {D) are all powers of D so that T'^ {D'^) = T^'^ {D). The stabilizer in (5.20) becomes 



T^'{D)T'2{D) 





7 



and the information-qubit matrix in (5.18) becomes 




r'2(7)-i) T^^D) 



r'2(D-i) 





We prcmultiply the first c — s rows of the stabilizer by r2 (7?^^) and add the result to the second 
c — s rows of the information-qubit matrix. These row operations from the stabilizer to the 
information-qubit matrix result in the information-qubit matrix having pure logical operators for 
the middle c — s qubits. Perform Hadamard gates on the second set of c — s qubits. The resulting 
information-qubit matrix is as follows 



r'2(7?-i) 
000 



000 
r^(7)-i) 



(5.21) 



so that the information qubits are available at the end of decoding. Processing may delay or 
advance them with respect to their initial locations because the matrix {D~^) is diagonal with 
powers of D. We can determine that the information qubits teleport coherently from the last set 
of c — s qubits to the second set of c — s qubits in every frame by comparing (5.21) to (5.16). 
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The overall procedure for encoding is to begin with a set of c ebits, 2 (n — c) — fci — fc2 ancilla 
qubits, and fci + fc2 — n+c information qubits. Alice performs the infinite-depth operations detailed 
in (5.14-5.19) for c — s of the ebits. Alice then performs the finite-depth operations detailed in the 
proofs of this lemma and Lemma 5.2.1 in reverse order. The resulting stabilizer has equivalent 
error-correcting properties to the quantum check matrix in (5.2). 

The receiver decodes by first performing all of the finite-depth operations in reverse order. 
The receiver then decodes the infinite-depth operations by the procedure listed in (5.20-5.21) so 
that the original ki + k2 — n + c information qubits per frame are available for processing at the 
receiving end. □ 

Example 5.4.1. Consider a classical convolutional code with the following check matrix: 



H{D) 



1 1 + D 



We can use the above check matrix in an entanglement-assisted quantum convolutional code to 
correct for both bit flips and phase flips. We form the following quantum check matrix: 



1 1 + D 






1 1 + D 



(5.22) 



We first perform some manipulations to put the above quantum check matrix into a standard 
form. Perform a CNOT from qubit one to qubit two in the same frame and in the next frame. 
The above matrix becomes 



D'^ + l + D 1 + D 






1 



Perform a Hadamard gate on qubits one and two. The matrix becomes 





1 



D-^ + l + D 1 + D 




93 



Perform a CNOT from qubit one to qubit two. The matrix becomes 





1 



D-'^ + l + D D-^ 




Perform a row operation that delays the first row by D. Perform a Hadamard on both qubits. 
The stabiUzer becomes 

" ' 



1 + £) + £>2 1 




1 



The above matrix is now in standard form. The matrix F (D) = 1 as in (5.7) so that its only 
invariant factor is equal to one. The code falls into the second class of entanglement-assisted 
quantum convolutional codes. We begin encoding with one ebit and one information qubit per 
frame. The stabilizer matrix for the unencoded stream is as follows: 



and the information-qubit matrix is as follows: 



110 














1 


1 





as follows: 

















1 


1 












Perform a Hadamard on qubit two and a CNOT from qubit two to qubit three so that the above 
stabilizer becomes 



1 











1 


1 





1 





1 









and the information-qubit matrix becomes 


















1 





1 


1 
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Perform an infinite-depth operation corresponding to the rational polynomial 1/ (l + D + _D^) on 
qubit three. Follow with a Hadamard gate on qubits two and three. The stabilizer matrix becomes 



1 1 1/ {1 + D + D^) 





110 



(5.23) 



and the information-qubit matrix becomes 



l/(l + £> + £)2) 






1 l + £)-i+£)-2 



(5.24) 



Perform the finite-depth operations above in reverse order so that the stabilizer becomes 



D 



1+-D 



l+D+D^ 1+D+D^ 







1 1 1 + D 



and the information-qubit matrix becomes 







l+D+P^ 1+P+P^ 











The above stabilizer is equivalent to the desired quantum check matrix in (5.22) by a row opera- 
tion that multiplies its first row by 1 + D + D"^. 



The receiver decodes by performing the finite-depth encoding operations in reverse order and 
gets the stabilizer in (5.23) and the information-qubit matrix in (5.24). The receiver performs a 
CNOT from qubit one to qubit two and follows with a CNOT from qubit two to qubit three in the 
same frame, in an advanced frame, and in a twice-advanced frame. Finally perform a Hadamard 
gate on qubits two and three. The stabilizer becomes 








1/{1 + D + D^) 
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and the information-qubit matrix becomes 




10 



1 l/(l + £> + £»2 




The receiver decodes the information qubits successfully because a row operation from the first row 
of the stabilizer to the first row of the information-qubit matrix gives the proper logical operators 
for the information qubits. Figure 7 of Ref. [52] details the above encoding and decoding operations 
for this entanglement-assisted quantum convolutional code. 



5.4.2 Discussion 

This second class of codes assumes that noiseless encoding is available. We require this assumption 
because the encoding circuit employs infinite-depth encoding operations. 

If an error does occur during the encoding process, it can propagate infinitely through the 
encoded qubit stream. The result of a single encoding error can distort both the encoded quantum 
information, the syndromes that result from measurements, and the final recovery operations based 
on the syndromes. 

We may be able to relax the noiseless encoding assumption if nearly noiseless encoding is 
available. The probability of an error would have to be negligible in order to ensure that the 
probability for a catastrophic failure is negligible. One way to lower the probability of an encoding 
error is to encode first with a quantum block code and then further encode with our quantum 
convolutional coding method. Many classical coding systems exploit this technique, the most 
popular of which is a Reed-Solomon encoder followed by a convolutional encoder. 



5.5 Closing Remarks 

This chapter develops the theory of entanglement-assisted quantum convolutional coding for CSS 
codes. We show several methods for importing two arbitrary classical binary convolutional codes 

for use in an entanglement-assisted quantum convolutional code. Our methods outline different 
ways for encoding and decoding our entanglement-assisted quantum convolutional codes. 
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Our first class of codes employs only finite-depth operations in their encoding and decoding 
procedures. These codes are the most useful in practice because they do not have the risk of 
catastrophic error propagation. An error that occurs during encoding, measurement, recovery, or 
decoding propagates only to a finite number of neighboring qubits. 

Our second class of codes uses infinite-depth operations during encoding. This assumption is 
reasonable only if noiseless encoding is available. The method of concatenated coding is one way 
to approach nearly noiseless encoding in practice. 

We suggest several lines of inquiry from here. Our codes are not only useful for quantum 
communication, but also should be useful for private classical communication because of the 
well-known connection between a quantum channel and private classical channel [39]. It may 
make sense from a practical standpoint to begin investigating the performance of our codes for 
encoding secret classical messages. The commercial success of quantum key distribution [12] for 
the generation of a private shared secret key motivates this investigation. It is also interesting to 
determine which entanglement-assisted codes can correct for errors on the receiver's side. Codes 
that possess this property will be more useful in practice. 

We hope that our theory of entanglement-assisted quantum convolutional coding provides 
a step in the direction of finding quantum codes that approach the ultimate capacity of an 
entanglement-assisted quantum channel. 
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Chapter 6 



Entanglement- Assisted Quantum Convolutional 
Coding: The General Case 



These patterns that anticommute, 
They laugh and they think it's so cute, 
Expand until fit, 
Then ole Gram and Schmidt, 
They'll line up with a grand salute. 

In this chapter, wo show how to encode and decode an entanglement-assisted quantum convo- 
hitional code that does not necessarily possess the CSS structure. The methods in this chapter 
represent a significant extension of our work on CSS entanglement-assisted quantum convolutional 
codes in the previous chapter. 

In particular, we develop a set of techniques that make the commutation relations for the 
generators of a general code the same as those of entangled qubits (ebits) and ancilla qubits. 
This procedure first "expands" the check matrix for a quantum convolutional code and applies 
an extended version of the symplectic Gram-Schmidt orthogonalization procedure [35] that incor- 
porates the binary polynomials representing quantum convolutional codes. We then show how 
to encode a stream of information qubits, ancilla qubits, and ebits to have the error-correcting 
properties of the desired code. We follow by detailing the decoding circuit that Bob employs at 
the receiving end of the channel. The algorithms for encoding and decoding use similar techniques 
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to those outlined in the previous chapter. The encoding circuits incorporate both finite-depth and 
infinite-depth operations as discussed in respective Sections 4.5 and 4.6 and the decoding circuits 
incorporate finite-depth operations only. 

One benefit of the techniques developed in this chapter is that the quantum code designer 
can produce an entanglement-assisted quantum convolutional code from a classical quaternary 
convolutional code. More generally, quantum convolutional code designers now have the freedom 
to design quantum convolutional codes to have desirable error-correcting properties without having 
to search for codes that satisfy a restrictive commutativity constraint. 

We structure this chapter as follows. In Section 6.1, we present an example that demonstrates 
how to expand the check matrix of a set of quantum convolutional generators and then show 
how to expand an arbitrary check matrix. We show in Section 6.2 how to compute the symplectic 
Gram-Schmidt orthogonalization procedure for the example in Section 6.1 and then generalize this 
procedure to work for an arbitrary set of quantum convolutional generators. The Gram-Schmidt 
orthogonalization technique reduces the quantum check matrix to have the same commutation 
relations as a set of ebits and ancilla qubits. This technique is essential for determining an encoding 
circuit that encodes a stream of information qubits, ancilla qubits, and ebits. Section 6.3 gives 
the algorithms for computing the encoding and decoding circuits for an arbitrary entangienient- 
assisted quantum convolutional code. We present a detailed example in Section 6.4 that illustrates 
all of the procedures in this chapter. We finish with a discussion of the practical issues involved 
in using these entanglement-assisted quantum convolutional codes and present some concluding 
remarks. 

The focus of the next two sections is to elucidate techniques for reducing a set of quantum 
convolutional generators to have the standard commutation relations in (3.2), or equivalently, the 
standard symplectic relations in (3.7). The difference between the current techniques and the 
former techniques [35, 34, 51, 78] is that the current techniques operate on generators that have a 
convolutional form rather than a block form. Section 6.1 shows how to expand a set of quantum 
convolutional generators to simplify their commutation relations. Section 6.2 outlines a symplectic 
Gram-Schmidt orthogonalization algorithm that uses binary polynomial operations to reduce the 
commutation relations of the expanded generators to have the standard form in (3.7). 
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6.1 The Expansion of Quantum Convolutional Generators 

We begin this section by demonstrating how to expand a particular generator that we eventually 
incorporate in an entanglement-assisted quantum convolutional code. We later generalize this 
example and the expansion technique to an arbitrary set of generators. This technique is important 
for determining how to utilize entanglement in the form of ebits in a quantum convolutional code. 

6.1.1 Example of the Expansion 

Let us first suppose that we have one convolutional generator: 



X Z I 



(6.1) 



This generator does not yet represent a valid quantum convolutional code because it anticommutes 
with a shift of itself to the left or to the right by one qubit. We are not concerned with the error- 
correcting properties of this generator but merely want to illustrate the technique of expanding 
it. 

Let us for now consider a block version of this code that operates on eight physical qubits with 
six total generators. The generators are as follows: 



X 


z 


I 


I 


I 


I 


I 


I 


I 


X 


z 


I 


I 


I 


I 


I 


I 


I 


X 


z 


I 


I 


I 


I 


I 


I 


I 


X 


z 


I 


I 


I 


I 


I 


I 


I 


X 


z 


I 


I 


I 


I 


I 


I 


I 


X 


z 


I 



(6.2) 
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We still use the vertical bars in the above block code to denote that the frame size of the code 
is one. Observe that we can view the frame size of the code as two without changing any of the 
error-correcting properties of the block code: 
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I 



(6.3) 



The frame is merely a way to organize our qubits so that we send one frame at a time over 
the channel after the online encoding circuit has finished processing them. We can extend the 
above block code with frame size two to have a convolutional structure with the following two 
convolutional generators; 



/ / 
/ / 



X Z 
I X 



I I 

Z I 



I I 
I I 



The above two generators with frame size two have equivalent error-correcting properties to the 
original generator in (6.1) by the arguments above. We say that we have expanded the original 
generator by a factor of two or that the above code is a two-expanded version of the original 
generator. We can also extend the original generator in (6.1) to have a frame size of three and we 
require three convolutional generators so that they have equivalent error-correcting properties to 
the original generator: 



III 
III 
III 



X Z I 
I X z 
I I X 



III 
III 

Z I I 



III 
III 
III 
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The representation of the original generator in (6.1) as a quantum check matrix in the poly- 
nomial formalism is as follows; 



9{D) 



D 



1 



(6.4) 



The two-expanded check matrix is as follows 

G2 (D) 

and the three-expanded check matrix is as follows: 






1 


1 





D 








1 



(6.5) 








1 





1 








Ga {D) = 








1 





1 







D 














1 



An alternative method for obtaining the polynomial representation of the two-expanded check 
matrix consists of two steps. We first multiply g {D) as follows: 



1 








D 





1 

1 



We then "plug in" the fractional delay operator and apply the flooring operation [-J that 
nulls the coefficients of any fractional power of D: 



G2{d)=[g',(d'/')\ 



A similar technique applies to find the check matrix of the three-expanded matrix. We first 
multiply g {D) as follows: 



Gs {D) 



1 

D 
£>2 



[gim 



1 5 




1 





1 ^ 



1 

£)2 
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We then "plug in" the fractional delay operator D^l^ and apply the flooring operation [-J : 



We discuss the general method for expanding an arbitrary check matrix in the next subsection. 
6.1.2 General Technique for Expansion 

We generalize the above example to determine how to expand an arbitrary quantum check matrix 
by a factor of I. Suppose that we have an n — A; x 2n quantum check matrix H (D) where 



H{D) 



Z{D) X{D) . 



Let D denote a diagonal matrix whose diagonal entries are the delay operator D. We take the 
convention that D° is the identity matrix and D™ is the matrix D multiplied m times so that its 
diagonal entries are D"*. Let Ri {D) and Ci {D) denote the following matrices: 



Ri {D) = 

Ci {D) ^ 



DO Di 
DO 

••• 



... D'-i 

D-C-i) 
DO 





D-C-i) 



where the diagonal D matrices in Ri (D) and C; {D) have respective dimensions n — k x n — k 
and n X n. We premultiply and postmultiply the matrix H (D) by respective matrices Ri {D) 
and Ci {D), evaluate the resulting matrix at a fractional power l/l of the delay operator D, and 
perform the flooring operation [.J to null the coefficients of any fractional power of D. The 
Z-expanded check matrix Hi (D) is as follows: 



Ht{D)=[Ri[D'/^)H[D'/^)Q[D'/') 



The Z-expanded quantum check matrix Hi (D) has equivalent error-correcting properties to the 
original check matrix. 
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6.2 Polynomial Symplectic Gram-Schmidt Orthogonaliza- 
tion Procedure 

In general, a given set of generators may have complicated commutation relations. We have to 
simplify the commutation relations so that we can encode information qubits with the help of 
ancilla qubits and halves of ebits shared with the receiver. In this section, we begin with an 
arbitrary set of convolutional generators. We show how to determine a set of generators with 
equivalent error-correcting properties and commutation relations that are the same as those of 
halves of ebits and ancilla qubits. We first show an example of the technique by illustrating it for 
Pauli sequences, for the quantum check matrix, and with the shifted symplectic product matrix. 
We then state a polynomial symplectic Gram-Schmidt orthogonalization algorithm that performs 
this action for an arbitrary set of quantum convolutional generators. 

6.2.1 Example of the Procedure 
Pauli Picture 

Let us consider again our example from the previous section. Specifically, consider the respective 
expressions in (6.1) and (6.3) for the convolutional generator and the block code. Recall that we 
can multiply the generators in a block code without changing the error-correcting properties of 
the code [32]. Therefore, we can multiply the sixth generator in (6.3) to the fourth. We can then 
multiply the modified fourth to the second to yield the following equivalent code: 
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(6.6) 
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Wc have manipulated the two-expanded matrix rather than the code in (6.2) because the com- 
mutation relations of the above code are equivalent to the commutation relations of the following 
operators: 
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We can use three ebits to encode the set of generators in (6.6) because they have the same 
commutation relations as the above operators and the above operators correspond to halves of 
three ebits. We resolve the anticommutation relations by using the following entanglement-assisted 
code: 
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The convention above is that the first qubit of each frame belongs to Bob and corresponds to 
half of an ebit. The second two qubits of each frame belong to Alice. The overall code forms 
a commuting stabilizer so that it corresponds to a valid quantum code. Bob could measure the 



105 



above operators to diagnose errors or he could measure the following operators that are equivalent 
by row operations; 
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One can check that the operators corresponding to the second two qubits of each frame are 
equivalent to the desired generators in (6.3). 

Polynomial Picture 

Let us consider the convolutional generator in (6.1). We now use the polynomial formalism because 
it is easier to perform the manipulations for general codes in this picture rather than in the Pauli 
picture. 

We extend the row operations from the above block code to the polynomial picture. Each row 
operation multiplied the even-numbered generators by themselves shifted by two qubits. Extending 
this operation to act on an infinite Pauli sequence corresponds to multiplying the second generator 
in (6.5) by the rational polynomial 1/ (1 + D). Consider the two-expanded check matrix with the 
operation described above applied to the second generator: 

91 (D) 
_ 92 (D) 

The shifted symplcctic products (gi Q gi) (D) = 0, (,gi 52) (-D) = 1, and (.92 .92) (-D) = 
capture the commutation relations of the resulting code for all frames. We can therefore resolve 
this anticommutativity by prepending a column that corresponds to Bob possessing half of an 
ebit: 

1 1 

. tId 



D 
1+D 



1 

1+D 



1 

1 1^5 . 
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Bob possesses the qubit corresponding to column one of both the "Z" and "X" matrix and Ahce 
possesses the two qubits corresponding to the second and third columns. The code above has 
equivalent error-correcting properties to those of the desired generators. 

The generators in the above polynomial setting correspond to Pauli sequences with infinite 
weight. Infinite- weight generators are undesirable because Bob cannot measure an infinite number 
of qubits. There is a simple solution to this problem and it is similar to what we did at the end 
of the previous subsection. We multiply the second row of the above check matrix by 1 + f to 
obtain the following check matrix that has equivalent error-correcting properties to the above one: 



1 1 
£> 



1 
1 + D 1 



Bob can measure operators that have finite weight because the above check matrix corresponds 
to Pauli sequences with finite weight. 

Shifted Symplectic Product Matrix Picture 

We can also determine the row operations that simplify the commutation relations by looking 
only at the shifted symplectic product matrix. It is easier to determine the row operations that 
resolve anticommutativity by considering the shifted symplectic product matrix. We would like to 
perform row operations to reduce the shifted symplectic product matrix to the standard form in 
(3.7) so that it has commutation relations equivalent to those of halves of ebits and ancilla qubits. 

The shifted symplectic product matrix corresponding to the check matrix in (6.4) is a one- 
element matrix Clg (D) where 

n,iD)=[D-' + D]. 

It is clear that Qg {D) is not reducible by row operations to the 2x2 matrix J in (3.5) because 
fig {D) is a one-element matrix. It is also not reducible to the one-element null matrix [0] because 
there is no valid row operation that can zero out the term + D. 
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Wc therefore consider the two-expanded check matrix in (6.5) to determine if we can reduce it 
with row operations to the standard form in (3.7). The shifted symplectic product matrix ilc^ 
for the two-expanded check matrix G2 {D) is as follows: 



l + 

l + D 



We can formulate the row operation of multiplying the second generator by 1/ (1 + D) as a matrix 
R {D) where 

1 

1/{1 + D) 

The effect on the matrix ficj (-D) is to change it to 



R{D) 



R{D) Qg, {D)R^{D-') = J 



as described in (4.39). The above matrix J has equivalent commutation relations to half of an 
ebit. We can therefore use one ebit per two qubits to encode this code. 

In a later section, we show how to devise encoding circuits beginning from c ebits per frame 
and a ancilla qubits per frame and also give the algorithm for their online decoding circuits. It 
is important for us right now to devise a procedure to reduce the commutation relations of any 
check matrix to those of ebits and ancilla qubits. 



6.2.2 The Procedure for General Codes 

We detail a polynomial version of the symplectic Gram-Schmidt orthogonalization procedure in 
this section. It is a generalized version of the procedure we developed for the above example. 
Before detailing the algorithm, we first prove a lemma that shows how to determine the shifted 
symplectic product matrix for an /-expanded version of the generators by starting from the shifted 
symplectic product matrix of the original generators. 
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The Shifted Symplectic Product Matrix for an /-Expanded Code 



Lemma 6.2.1. Suppose the shifted symplectic product matrix fl {D) of a given check matrix H (D) 
is as follows: 

Q{D) = Z {D) + X {D) Z'^ (D'^) . 

The shifted symplectic product matrix fli {D) of the l-expanded check matrix Hi (£>) is as follows: 



where the flooring operation [-J nulls the coefficients of any fractional power of D. 

Proof. Consider that the "X" matrix Xi (D) of the /-expanded check matrix Hi (D) is as follows: 



Xi {D) = [Ri (dV; j X j c'l {d'^')\ 



where 



C'l (d^/^) ^ 



DO D-i ••• D-('-i) 



and each diagonal D matrix is n x n-dimensional. It is also then true that the "Z" matrix of 
Hi {D) is as follows: 

Zi {D) = [Ri (d^/') Z (i)!/') C'l (i?!/') . 

The matrix transpose operation, the time reversal operation (substituting for D), and matrix 
addition are all invariant under the flooring operation [-J for arbitrary matrices M (D) and N (D): 



D=D-^ 



Additionally, the following property holds for two arbitrary binary polynomials / {D) and g {D) : 



i=0 



(6.7) 
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Now consider the product Xi [D) Zf i^D ^): 



X,{D)Zf {D-') = 



Ri\D 



D=D- 



= [Ri (d^/') X (d^') C[ {d^")\ [c[^ (d-^') Z'^ (d-^/') Rf 
= ^ [d-^/'Ri X (l)i/')J [d'/'Z^ Rf {d-'^')\ 

i=0 

= [Ri (z)!/') X Z^ Rf 



where the second Une uses the invariance of the flooring operation with respect to matrix transpo- 
sition and time reversal, the third line expands the matrix multiplications using the matrix C; (D) 
defined above, and the last line uses the matrix generalization of the multiplication property de- 
fined in (6.7). Our final step is to use the invariance of the flooring operation with respect to 
matrix addition: 



ni {D) = Xi {D) zf {D-^)+ Zi (D) Xf (D-^) 

= Ri (^DT^ X (^D^^ Z'^ (^D^^ Rf (^D^^ + Ri (^Dt^ Z (^D^^ X^ (^D^^ Rf (^D^^ 

= Ri (£)t) X (dt j Z^ (^^) {^^) ^ {^^) ^ {^^) {^^) {^^)\ 

= Ri (d^") n (z)!/') Rf 



□ 



The Gram-Schmidt Procedure 

Wc now present the polynomial symplectic Gram-Schmidt orthogonalization procedure that re- 
duces the commutation relations of a given set of convolutional generators to have the standard 
form in (3.7). 

Consider the following n — k x 2n-dimensional quantum check matrix H (D): 



H{D)= Z{D) X{D) . 



Label each row as hi (D) 



z, [D) Xi {D) 



for alH e {1, . . . , n — k}. 
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We state the Gram-Schmidt procedure in terms of its effect on the shifted symplectic product 
matrix. It is easier to see how the algorithm proceeds by observing the shifted symplectic product 
matrix rather than by tracking the generators in the check matrix. 

Let I denote the amount by which we expand the check matrix H (D). Suppose first that I — 1 
(we do not expand the check matrix). Let us say that the check matrix has r generators (we take 
a number difi^erent from n — k because the number of generators may not be equal to n — A; for 
future iterations). There are the following possibilities: 

(i). There is a generator hi (D) such that {hi hj) (D) = for all j G {1, • • • ,r}. In this case, 
the generator is already decoupled from all the others and corresponds to an ancilla qubit 
because an ancilla and it share the same commutation relations. Swap hi (D) to be the first 
row of the matrix so that it is hi (D). The shifted symplectic product matrix then has the 
form: 







/l2,2 



T-2,r 



h 



r,2 



hr 



i2,2 



1,r 



h. 



r,2 



hr 



where [0] is the one-element zero matrix and we use the shorthand hij = {hi hj) {D). We 
remove generator hi {D) from check matrix H (D) and continue to step two below for the 
remaining generators in the matrix. 

(ii). There are two generators hi (D) and hj (D) for which {hi hj) (D) = D"^ for some integer 
m and {hi hi) {D) = {hj hj) {D) = 0. In this case these generators correspond exactly 
to half of an ebit. Multiply generator hj {D) by D"^. This row operation has the effect of 
delaying (or advancing) the generator by an amount m and changes the shifted symplectic 
product to be {hi hj) {D) = 1. These two generators considered by themselves now have 
the commutation relations of half of an ebit. Swap the generators hi {D) and hj {D) to be 
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the first and second respective rows of the check matrix H (D). Call them hi {D) and /12 (D) 
respectively. The shifted symplectic product matrix is then as follows: 






1 


hi,3 ■ 




1 





h2,3 ■ 


/l2,r 


3,1 


h3,2 


h3,3 ■ 


h3,r 


r,2 


hr,2 


hr,3 • 





We use the following row operations to decouple the other generators from these two gener- 
ators: 

h', (D) = h (D) + [{h, © h2) (D)] hi (D) + [{h, © hi) {D)] h2 (D) for alH G {3, . . . , r} . 

The shifted symplectic product matrix becomes as follows under these row operations: 

1 ••• 

1 ••• 

/l'3,3 ••• h's^, 

/i;_3 ••• h'^^^ 

The first two generators are now decoupled from the other generators and have the commu- 
tation relations of half of an ebit. We remove the first two generators from the check matrix 
so that it now consists of generators /13 {D), . . ., h'^ {D). We check to see if the conditions at 
the beginning of the previous step or this step hold for any other generators. If so, repeat 
the previous step or this step on the remaining generators. If not, see if the conditions for 
step three hold. 

There are two generators hi {D) and hj {D) for which (/i, hi) {D) = {hj hj) {D) = 
but {hi hj) {D) ^ £)™ for all m and {hi hj) {D) ^ 0. Multiply generator hj {D) by 



J 



''■3, 3 ■ ■ ■ 



'3,r 



^r,3 ' " ^'r, 
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1/ {hj hi) (D). Generator hj (D) becomes infinite weight because {hj hi) (D) is a poly- 
nomial with two or more powers of D with non-zero coefficients. Now the shifted symplectic 
product relations are as follows: {hi hi) (D) = {hj hj) (D) = and {hi hj) {D) = 1. 
We handle this case as we did the previous case after the two generators there had the 
commutation relations of half of an ebit. 

(iv). None of these conditions hold. In this case, we stop this iteration of the algorithm and 
expand the check matrix by the next factor I := I + 1 and repeat the above steps. 

We have not proven that this procedure converges on all codes; however, it does converge on 
all the codes we have tried. We conjecture that this procedure converges for all codes, but even 
if this is true, in principle it might require expansion to a large number of generators. A simple 
and practical convergence condition is as follows. In practice, convolutional codes do not act on 
an infinite stream of qubits but instead act on a finite number of qubits. It may be that there 
are codes for which this procedure either does not converge or must be expanded until the frame 
size of the expanded code exceeds the number of qubits that the code acts on. In this case, we 
would not employ this procedure, and instead would treat the code as a block code, where we 
could employ the methods from Refs. [34, 35] for encoding and decoding. It is unclear if this 
practical convergence condition will ever really be necessary. This procedure does converge for a 
large number of useful codes, so that the frame size of the expanded code is much less than the 
number of qubits that the code acts on, and we have not found an example where this procedure 
fails. 

6.3 Encoding and Decoding Circuits 

This section proves the main theorem of this chapter. The theorem assumes that we have al- 
ready processed an arbitrary check matrix with the Gram-Schmidt algorithm and that the shifted 
symplectic product matrix corresponding to the processed check matrix has the standard form in 
(3.7). The theorem shows how to encode a set of information qubits, ancilla qubits, and halves 
of ebits into a code that has equivalent error-correcting properties to those of a desired set of 
convolutional generators. The theorem uses both finite-depth and infinite-depth operations in the 
encoding circuit and finite-depth operations in the decoding circuit. 
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Theorem 6.3.1. Suppose we have a set of quantum convolutional generators in the n — k x 2n- 
dimensional matrix H (D) where 



H{D) 



Z{D) X{D) 



Its shifted symplectic product matrix fl {D) is as follows: 



Q{D) = Z {D) (£)-!) + X {D) Z^ 



Suppose check matrix H (D) is the matrix resulting from processing with the polynomial symplectic 
Gram-Schmidt orthogonalization procedure. Therefore, Q{D) has the standard form in (3.7) with 

parameters c and a = n — k — 2c. Then there exists an online encoding circuit for the code that 
uses finite-depth and infinite- depth operations in the shift-invariant Clifford group and there exists 
an online decoding circuit for the code that uses finite-depth operations. The code encodes /c + c 
information qubits per frame with the help of c ebits and a = n — k — 2c ancilla qubits. 

Proof. We prove the theorem by giving an algorithm to compute both the encoding circuit and 
the decoding circuit. The shifted symplectic product matrix fl (D) for check matrix H (D) is in 
standard form so that the first 2c rows have the commutation relations of c halves of ebits and the 
last a rows have the commutation relations of a ancilla qubits. Perform the algorithm outlined in 
Refs. [46, 45] on the last a generators that correspond to the ancilla qubits. The algorithm uses 
finite-depth CNOT gates, Hadamard gates, and phase gates. The resulting check matrix has the 
following form; 



(6.8) 



Z"{D) Z'{D) X'{D) 
T{D) 

where the matrix Z" (D) and the null matrix at the top left of the "X" matrix each have dimension 
2cx a, the matrices Z' (D) and X' (D) each have dimension 2cx n— a, and all matrices in the second 
set of rows each have a rows. The matrix F (D) may have entries that are rational polynomials. 
If so, replace each of these entries with a "1" so that the resulting matrix has the following form: 



Z" {D) Z' {D) 
I 



X' (D) 
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This replacement is equivalent to taking a subcode of the original that has equivalent error- 
correcting properties and rate [45] . We can also think of the replacement merely as row operations 
with rational polynomials. We then perform row operations from the last a rows to the first 2c 
rows to obtain the following check matrix: 



Z' (D) 

1 



X'{D) 




The shifted symplectic product matrix still has the standard form in (3.7) because these last row 
operations do not change its entries. We now focus exclusively on the first 2c rows because the 
previous steps decoupled the last a rows from the first 2c rows. Consider the following submatrix: 



H'iD) 



Z{D) X{D) 



where we have reset variable labels so that Z (D) = Z' (D) and X (D) = X' (D). Perform row 
permutations on the above matrix so that the shifted symplectic product matrix for H' (D) changes 
from 



to become 



i=l 

I 

1 



(6.9) 



where each identity and null matrix in the above matrix are c x c-dimensional. We can employ 
the algorithm from Ref. [46] on the first c generators because the first c rows of the resulting 
check matrix form a commuting set (we do not use the row operations in that algorithm). The 
algorithm employs finite-depth CNOT gates, Hadamard gates, and phase gates and reduces the 
check matrix to have the following form; 





U{D) Z2{D) 



L{D) 
X,{D) X2{D) 
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where L (D) is a c x c lower triangular matrix and U [D] is a c x c upper triangular matrix. The 
z"^ diagonal entry uu (D) of U (D) is equal to l/lu (D^^) where la (D) is the i*"^ diagonal entry 
of L[D). This relationship holds because of the shifted symplectic relations in (6.9) and because 
gates in the shift-invariant Clifford group do not affect the shifted symplectic relations. We now 
employ several row operations whose net effect is to preserve the shifted symplectic relations 
in (6.9) — we can therefore include them as a part of the original polynomial symplectic Gram- 
Schmidt orthogonalization procedure. Multiply row i of the above check matrix by l/lu {D) and 
multiply row i + c by l/uu {D) for all i G {1, . . . ,c}. Then use row operations to cancel all the 
off-diagonal entries in both L (D) and U {D). The resulting check matrix has the following form: 





1 Z'^{D) 



I 
X[{D) X'^{D) 



where the primed matrices result from all the row operations. One can check that the shifted 
symplectic relations of the above matrix are equivalent to those in (6.9). Perform Hadamard gates 
on the first c qubits. The check matrix becomes 



I 







(6.10) 



We show how to encode the above matrix starting from c ebits and k + c information qubits. The 
following matrix stabilizes a set of c ebits: 



I 


I 























I 


I 






(6.11) 



where each identity matrix is c x c and the last column of zeros in each matrix is c x (/c + c). The 
receiver Bob possesses the first c qubits and the sender Alice possesses the last k + 2c qubits. The 
following matrix is the information-qubit: 









/ 


























I 
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where each identity matrix is (fc + c) x (fc + c) and each cohinin of zeros is (fc + c) x c. It is 
important to track the information-qubit matrix throughout encoding and decoding so that we 
can determine at the end of the process if we have truly decoded the information qubits. Perform 
finite-depth CNOT operations from the first c ebits to the last k+c qubits to encode the numerators 
of the entries in matrix Z2 {D) . Let Z2 jv denote this matrix of the numerators of the entries 
in Z'2 (D) . The stabilizer matrix becomes 



I/O 


and the information-qubit matrix becomes 







Z^^^iD) I 





0/ 



where ^2,Ar (-^) matrix that results on the "Z" side after performing the CNOT operations 

corresponding to the entries in Z2 {D). Perform Hadamard gates on the last k + c qubits. The 
stabilizer matrix becomes 



/ / 
Z'2^^{D) 

and the information-qubit matrix becomes 




7/0 



Z^,,{D) 
0/ 



0/ 




Let X2 jv (D) denote the matrix whose entries are the numerators of the entries in X2 (D). Perform 
CNOT gates from the first c qubits to the last k+c qubits corresponding to the entries in X2 {D). 
The stabilizer matrix becomes 



/ / 
A{D) Z'{D) 




/ / X'{D) 
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where A (D) = -2^2 w (-^) -^2 n i-^) ^^"^ -^2 n i-^) ^hc matrix that results on the "Z" side after 
performing the CNOT operations on the "X" side corresponding to the entries in X2 jv The 
information-qubit matrix becomes 



Z^j^iD) 
In{D) I 



X'^ 



0/ 





Let r (D) be a diagonal matrix whose i^^ diagonal entry is the denominator of the i*"^ row of 
Z2 (-D) and X2 {D). We perform infinite-depth operations corresponding to the entries in F (D). 
The stabilizer matrix becomes 



I r-i (£)-!) 
AiD) Z'2^^{D) 

and the information-qubit matrix becomes 





1 T{D) X'{D) 



Z'l^{D)T-^{D-^) 
X'l^{D)T-^{D-^) I 



0/ 




The above stabilizer matrix is equivalent to the desired one in (6.10) by several row operations. 
We premultiply the first set of rows by F (/?~^) and multiply the second set of rows by (D). 
We can also use the resulting identity matrix in the first set of rows to perform row operations 
from the first set of rows to the second set of rows to realize the matrix X'l (D) . The operators 
that Bob would really measure need to have finite weight so he would measure the operators 
corresponding to the entries in the following stabilizer matrix: 



r / 

A{D) Z'2^^{D) 





/ T{D) X'{D) 



(6.12) 



We are done with the encoding algorithm. Alice begins with a set of ebits and performs the 
encoding operations detailed in (6.11-6.12) and then performs the finite-depth operations detailed 
in (6.8-6.10) in reverse order. We now detail the steps of the decoding algorithm. Perform CNOT 
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gates corresponding to the entries in X2 jq {D) from the first set of c qubits to the last set of fc + c 
qubits. The stabihzer matrix becomes 



I r-i 

B{D) A{D) Z'{D) 





/ T{D) 



where B (D) = Z'2 ^ (D) X2 (D). The information-qubit matrix becomes 



Z'^^^{D)T-^D-^) 
X'i^iD) X!{^^{D)T-^D-') I 



1 




Perform Hadamard gates on the last set of A; + c qubits. The stabilizer matrix becomes 



I T-^{D-^) 
B (D) A (D) 

and the information-qubit matrix becomes 







Z-^^{D)T-^D-^) I 
X'iNi.D) X'l^{D)T-'{D-') 




0/ 



Perform CNOT gates from the first set of c qubits to the last set of fc + c qubits. These CNOT 
gates correspond to the entries in Z2 (-D). The stabilizer matrix becomes 



B (D) A (£>) 




/ T{D) 



and the information-qubit matrix becomes 



X'l^iD) X'l^{D)T-^{D-^) 




7 
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Row operations from the first set of rows of the stabihzer to each of the two sets of rows in the 
information-qubit matrix reduce the information-qubit matrix to the following form: 









I 


























I 



Then we perform the finite-depth operations detailed in (6.8-6.10). We have now finished the 
algorithm for the decoding circuit because the logical operators for the information qubits appear 
in their original form. □ 



6.3.1 Discussion 

Similar practical issues arise in these circuits as we discussed in the previous chapter. Encoding 
circuits with infinite-depth operations are acceptable if we assume that noiseless encoding is pos- 
sible. Otherwise, infinite-depth operations could lead to catastrophic propagation of uncorrected 
errors. Noiseless encoding is difficult to achieve in practice but we may be able to come close to 
it by concatenation of codes at the encoder. 

There is a dichotomy of these codes similar to that in the previous chapter. Some of the codes 
may have a simpler form in which the encoding circuit requires finite-depth operations only. These 
codes fall within the first class of codes discussed in the previous chapter and will be more useful 
in practice because they propagate errors in the encoding circuit to a finite number of qubits only. 
The remaining codes that do not have this structure fall within the second class of codes whose 
encoding circuits have both finite-depth and infinite-depth operations and whose decoding circuits 
have finite-depth operations only. 

6.3.2 Importing Classical Convolutional Codes over GF (4) 

One benefit of the new entanglement-assisted quantum convolutional codes is that we can produce 
one from an arbitrary classical convolutional code over GF (4). The error-correcting properties of 
the classical convolutional code translate to the resulting quantum convolutional code. It is less 
clear how the rate translates because we use the expansion technique. We know that the term 
{2k — ri) jn lower bounds the "entanglement-assisted" rate where n and k are the parameters from 



120 



the imported classical code. The rate should get a significant boost from entanglement — the rate 
boosts by the number of ebits that the code requires. 

The construction for importing an [n, k] classical convolutional code over GF (4) is as follows. 
Suppose the check matrix for the classical code is an n — A: x n-dimensional matrix H (D) whose 
entries are polynomials over GF (4). We construct the quantum check matrix H (D) according to 
the following formula: 



where 7 denotes the isomorphism between elements of GF(4) and symplectic binary vectors de- 
tailed in (3.38). We use this construction for the example in the next section. 



We take the convolutional generators from Ref. [51] as our example. Ref. [34] originally used these 
generators in a block code. We import the following classical convolutional code over GF (4): 



We produce two quantum convolutional generators by multiplying the above generator by uj and 
CO and applying the map in (3.38). The resulting quantum convolutional generators are as follows: 




6.4 Example 



(•••|0000|lwlO|1101|0000|---)- 



(6.13) 



(• • • \IIII\ZXZI\ZZIZ\IIII\ ■■■), 



(• • • \IIII\XYXI\XXIX\IIII\ ■■■). 



These generators have the following representation in the polynomial formalism: 



1 + D DID 







1 











10 



1 + D 1 + D 1 D 



The shifted symplectic product matrix Q {D) for the above code is as follows: 



D + D-^ 



D-^ 



fl{D) = 



D 



D + D-^ 
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The above matrix is not reducible to the standard form by any row operations. We therefore 
expand the code by a factor of two to give four generators with a frame size of eight. The two- 
expanded check matrix H (D) is as follows: 



H{D) 



1 0101101 

1 

13 D D 1 1 

1 



1000000 

1 110 110 1 
1 
DDODlllO 



The shifted symplectic product matrix ^2 (D) of the two-expanded check matrix is as follows: 



^^2 {D) 






1 + D 
D 





1 

1 + D 



1 + D-i 

1 1 + D-i 












We proceed with the Gram-Schmidt procedure because this matrix satisfies its initial requirements. 
We swap generators two and three to be the first and second generators of the check matrix because 
they have the commutation relations of half of an ebit. The shifted symplectic product matrix 
becomes 

" 1 l + £>-i 

1 1 + D 

1 + D-i 
1 + D D 

Multiply generator two by 1 + _D and add to generator four. Multiply generator one by 1 + D~^ 
and add to generator three. The shifted symplectic matrix becomes 






1 








1 

















14 


-D-^ + D-^ 





IH 


^D + D^ 
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Wc finally divide generator four by 1 + D + and the shifted symplectic product matrix then 
becomes 



so that it has the commutation relations of halves of two ebits. The check matrix resulting from 
these operations is as follows: 



H2iD) 



Z2 {D) X2 (D) 



(6.14) 



where 



Z2{D) 






1 




















D 


D 





D 


1 





1 





1 




1 





1 


1 





1 


d'^+d 


D^+D 





d'^+d 


D+1 


1 


D+1 















D^+D+1 



X2 (D) = 



D 



1 



; + 
D 



1 



D 
D 



D 



L D^+D+1 








D 



D 



D 



1 

1 

+ 
D 



D^+D+1 D^+D+1 D^+D+1 D^+D+1 



D 



The error-correcting properties of the above check matrix are equivalent to the error-correcting 
properties of the original two generators. 

This code sends six information qubits and consumes two ebits per eight channel uses. The 
rate pair for this code is therefore (3/4, 1/4). 

We can now apply the algorithm in Theorem 6.3.1 to determine the encoding and decoding 
circuits for this code. The encoding circuit begins from a set of two ebits and eight information 
qubits per frame with the following stabilizer matrix: 



Ho{D) 



1010000000 

0101000000 
0000000000 
0000000000 



0000000000 

0000000000 
1010000000 
0101000000 
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Wc label the eight qubits on the right side of each matrix above as 1, . . . , 8. We label Bob's two 
qubits on the left as i?l and B2. Perform the following finite-depth operations (in order from left 
to right and then top to bottom): 

C (1, 4, £> + 1)2) C (1, 5, 1 + D^) C (1, 6, 1) C (1, 7,1 + D)C (2, 4, D) C (2, 5, 1 + D) C (2, 6, 1) 
H (3, . . . , 8) C (1, 4, L> + £»2 ^ L>4) C (l, 5, D^) C (1, 6, 1 + D) C (l, 7, D^) C (2, 4, £) + D^) 
C (2, 5, 1 + £>) C (2, 6, 1) C (2, 7, 1 + £>) 

where we use the notation C {qi,q2,f {D)) to represent a finite-depth CNOT gate from qubit one 
to qubit two that implements the polynomial / (£>), H {qi, . . . ,qj) is a sequence of Hadamard gates 
applied to qubits qi through qj in each frame, P {q) is a phase gate applied to qubit q in each 
frame, and C {q, l/f{D)) is an infinite-depth CNOT gate implementing the rational polynomial 
l/f{D) on qubit q. Alice performs the following infinite-depth operations: 

She then finishes the encoding circuit with the following finite-depth operations: 

H (1, 2) C (2, 3, 1) C (2, 5, 1) C (2, 6, 1) C (2, 8, 1) P {2) H {3, . . . ,8) S (2, 3) (6.15) 
C (1, 2, 1) C (1, 3, 1) C (1, 5, 1) C (1, 6, 1) C (1, 8, 1) P (2) 

The code she encodes has equivalent error-correcting properties to the check matrix in (6.14). 

Bob performs the following operations in the decoding circuit. He first performs the operations 
in (6.15) in reverse order. He then performs the following finite-depth operations: 

C {Bl, A,D + D'^ + £)") c (Bl, 5, D^) C{B1,6,1 + D)C {Bl, 7, D^) C {B2, 4,D + D'^) 
C {B2, 5, 1 + £>) C (B2, 6, 1) C {B2, 7,1 + D) H {3, . . . ,8) 
C {Bl, 4,0 + D^) C {Bl, 5, 1 + D'^)C {Bl, 6, 1) C{B1,7,1 + D)C {B2, 4, D) 
C{B2,5,1 + D)C{B2,6,1) 
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The information qubits then appear at the output of this onhne decoding circuit. 



6.5 Optimal Entanglement Formulas 

We conjecture two formulas for the optimal number of ebits that a general (non-CSS) entanglement- 
assisted quantum convolutional code or one imported from a classical quaternary convolutional 
code need to consume per frame of operation. We show that this conjecture holds for a particular 
example. 

Conjecture 6.5.1. The optimal number c of ebits necessary per frame for an entanglement- 
assisted quantum convolutional code is 

c = rank {Hx{D)H^{D-^) + Hz{D)H]^{D-^)) /2 (6.16) 

where H{D) = Hz{D) Hx{D) represents the parity check matrix for a set of quantum 
convolutional generators that do not necessarily form a commuting set. 

It is clear that the above formula holds after we have expanded an original set of generators 
and their commutation relations are reducible by the polynomial symplectic Gram-Schmidt or- 
thogonalization procedure to the standard form in (3.7). The proof technique follows from the 
proof technique outlined in Section 3.6. But we are not sure how to apply this formula to an 
initial set of unexpanded generators. 

Conjecture 6.5.2. The optimal num,ber c of ebits required per frame for an entanglement- assisted 
quantum convolutional code imported from a classical quaternary convolutional code with parity 
check matrix H [D) is 

c = lank {H{D)H\D-^)) . (6.17) 

We know the number of ebits required for a CSS entanglement-assisted quantum convolutional 
code is rank (i?i(D)ilJ(D~^)) where Hi{D) and H2{D) correspond to the classical binary con- 
volutional codes that we import to correct respective bit and phase flips. Comparing the CSS 
convolutional formula, the CSS block formula in Corollary 3.6.1, and the general block formula in 
Theorem 3.6.1, the above conjecture seems natural. 
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Wc finally provide an example of the above conjecture. It is a slight modification of the code 
presented in Ref. [53]. 

Example 6.5.1. Consider the quantum convolutional code with quantum check matrix as follows: 


















h{D) 





D 


1 


h{D 


hiD) 


D 





1 


hiD) 























D 


D 


D 





1 





1 


1 





1 

D 


1 


1 

D 























1 

D 

















1 









where h (D) ~ 1 + D. This code requires two cbits and one ancilla qubit for quantum redundancy 
and encodes two information qubits. The shifted symplectic product matrix Hx{D)H^(D~^) + 
Hz{D)H'^{D-^) [51, 52] for this code is as follows: 

10 

1 

(6.18) 

1 

1 

The rank of the above matrix is four and the code requires two ebits. Therefore, the conjecture 
holds for this example. 



6.6 Closing Remarks 

There are several differences between the methods used for general, non-CSS codes discussed in 
this chapter and the CSS codes used in the previous chapter. It was more straightforward to 
determine how to use ebits efficiently in CSS entanglement-assisted quantum convolutional codes, 
but we have had to introduce the expansion technique in Section 6.1 in order to determine how to 
use ebits efficiently for codes in this chapter. There was also no need for an explicit Gram-Schmidt 
orthogonalization procedure in the previous chapter. The Smith algorithm implicitly produced 
the row operations necessary for symplectic orthogonalization. 
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Wc do have some methods in the next chapter that do not require expansion of a check 
matrix or an explicit Gram-Schmidt procedure but these methods do not make efficient use of 
entanglement and have a lower rate of quantum communication and higher rate of entanglement 
consumption than the codes discussed in this chapter. Nonetheless, we have determined ways to 
make these other codes more useful by encoding classical information in the extra entanglement 
with a superdense-coding-like effect [13]. These other codes are grandfather codes in the sense of 
Refs. [53, 84] because they consume entanglement to send both quantum and classical information. 
We discuss these techniques in the next chapter. 

One negative implication of the expansion of a code is that the expanded code requires more 
qubits per frame. The expanded code then requires a larger buffer at both the sender's and 
receiver's local stations. The increased buffer will be a concern right now because it is difficult to 
build large quantum memories. This issue will become less of a concern as quantum technology 
advances. The entanglement-inefficient codes in the next chapter have the advantage that they 
do not require expansion and thus require smaller buffers. It therefore should be of interest 
to find solutions in between the entanglement-efficient codes discussed in this chapter and the 
entanglement-inefficient codes discussed in the previous paragraph. 

Some outstanding issues remain. We have not proven the convergence of the polynomial sym- 
plectic Gram-Schmidt orthogonalization procedure and have instead provided a practical stopping 
condition. We have a conjecture for how to proceed with proving convergence. Suppose that we 
would like to construct a code consisting of one generator that does not commute with shifts of 
itself. We have found for many examples that the correct expansion factor for the generator is 
equal to the period of the inverse polynomial of the generator's shifted symplectic product. We 
do not have a proof that this factor is the correct one and we are not sure what the expansion 
factor should be when we would like to construct a code starting from more than one generator. 

The techniques developed in this chapter represent a useful way for encoding quantum infor- 
mation. This framework will most likely lead to codes with reasonable performance but they will 
most likely not come close to achieving the quantum capacity of an entanglement-assisted channel. 
The next step should be to combine the theory in this chapter with Poulin et al.'s recent theory of 
quantum serial-turbo coding [58] . This combination might lead to entanglement-assisted quantum 
turbo codes that come close to achieving capacity. 
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Chapter 7 



Entanglement- Assisted Quantum Convolutional 
Coding: Free Entanglement 



Sometimes we 've a free cater, 
Perhaps on the house we've a waiter, 
But free entanglement you say? 

Not a price to pay? 
Oh assume it for now and pay later. 

Wc mentioned in the previous chapter that the expansion of a set of quantum convolutional 
generators increases the frame size of a code. This increase implies that each round of transmission 
in the protocol sends a larger number of encoded qubits and requires a larger quantum memory 
for its operation. Thus increasing the frame size is somewhat undesirable. 

We offer the "free-entanglement" codes in this chapter as an alternative option to the codes 
in the previous chapter. The free-entanglement codes have a basic generator set that includes 
an arbitrary set of Pauli sequences. The free-entanglement codes consume more entanglement 
than is necessary and have a lower quantum information rate, but the benefit is that they do not 
require a larger quantum memory for each round of transmission. It is up to the quantum coding 
engineer to decide which option is more desirable for the particular coding application: using 
less entanglement and having a higher rate of quantum transmission or using a smaller quantum 
memory to send fewer qubits for each round of transmission. 
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Wc make the additional assumption that shared entanglement is freely available. Quantum 
information theorists make the "free entanglement" assumption when deriving the entanglement- 
assisted capacity of a quantum channel [40, 41]. This model makes sense when the sender and 
receiver can use the noisy channel and entanglement distillation protocols [65, 66] during off-peak 
times to generate a large amount of noiseless entanglement. 

The "entanglement-assisted" rate (Section 3.3) of a quantum code is the ratio of the number 
of encoded information qubits to the number of physical qubits assuming that entanglement is 
available for free. The entanglement-assisted rate of our codes is at least k/n. Parameter k is the 
number of information qubits that the code encodes and parameter n is the number of qubits used 
for encoding. A basic set of n — A; generators specifies all of our entanglement-assisted quantum 
convolutional codes and yields the above entanglement-assisted rate. 

We say that the entanglement-assisted rate is at least k/n because the example in Section 7.2 
discusses a method to boost the rate. This method encodes some classical information into the 
extra entanglement along the lines of the "grandfather" coding technique discussed in the next 
chapter and in Ref. [84]. The sender can use this classical information and the entanglement to 
teleport [3] additional information qubits and increase the rate of quantum transmission. 

7.1 Construction 

The proof of our main theorem below outlines how to encode a stream of information qubits, 
ancilla qubits, and shared ebits so that the encoded qubits have the error-correcting properties of 
an arbitrary set of Paulis. The receiver may employ an error estimation algorithm such as Viterbi 
decoding [2] to determine the most likely errors that the noisy quantum communication channel 
induces on the encoded stream. We then show how to decode the encoded qubit stream so that 
the information qubits become available at the receiving end of the channel. 

The encoding circuits in the proof of our theorem employ both finite-depth and infinite-depth 
operations. The decoding circuits employ finite-depth operations only. Infinite-depth operations 
can lead to catastrophic error propagation as discussed in previous chapters. In our proof below, 
we restrict infinite-depth operations to act on qubits before sending them over the noisy channel. 
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Catastrophic error propagation does not occur under the ideal circumstance when the operations 
in the encoding circuit are noiseless. 

Our theorem below begins with a "riuantum check matrix" that consists of a set of Pauli 
sequences with desirable error-correcting properties. This quantum check matrix does not neces- 
sarily correspond to a commuting stabilizer. The proof of the theorem shows how to incorporate 
ebits so that the sender realizes the same quantum check matrix for her qubits and the sender and 
receiver's set of generators form a valid commuting stabilizer. 

Theorem 7.1.1. Suppose the following quantum check matrix 



S{D) 



Z{D) X{D) 



e F2 [D] 



(n— fc)x2n 



where S {D) is of full rank and does not necessarily form a commuting stabilizer and ¥2 [D] is 
the field of binary polynomials. Then an entanglement-assisted quantum convolutional code exists 
that has the same error- correcting properties as the above quantum check matrix S{D). The 
entanglement- assisted rate of the code is at least k/n. 

Proof. Suppose that the Smith form [80] of X (D) is as follows 



X{D) = A{D) 



Fi (D) 
T2 {D) 




B{D), 



(7.1) 



where A (D) is (n — k) x (n — k), B [D] in n x n, Ti [D] is an ,s x ,s diagonal matrix whose entries 
are powers of D, and r2 {D) is a (c — s) x (c — s) diagonal matrix whose entries are arbitrary 
polynomials. Perform the row operations in (D) and the column operations in (D) on 
S{D). The quantum check matrix S {D) becomes 



E{D) 



Fi [D) 
F2 (£>) 




(7.2) 
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where E [D) = A^^ (L>) Z [D) [D-^). Suppose Ei [D) is the first c cohinins of E [D) and 
E2 (-D) is the next n — c columns of E (D) so that the quantum check matrix is as follows: 



Ei{D) E2{D) 



Ti (D) 
r2 {D) 




(7.3) 



Perform Hadamard gates on the last n — c qubits so that the quantum check matrix becomes 



El {D) 



Ti (£)) E2,i {D) 
Fa {D) E2,2 {D) 
E2,3 {D) 



(7.4) 



where 



E2A (D) 

E2{D)^ E2.2{D) ■ (7.5) 
£^2,3 {D) 

Perform CNOT operations from the first s qubits to the last n — c qubits to clear the entries in 
-£^2,1 {D). The quantum check matrix becomes 



El {D) 



Ti (D) 
T2 {D) E2,2 (£>) 
E2,3 (£>) 



(7.6) 



The Smith form of £^2,3 (-D) is as follows 



E2,3{D) = Ae{D) r{D) Be{D), 



(7.7) 



where Ae (D) is {n — k — c) x {n — k — c), Be (D) is (n — c) x (n — c), and T (D) is a (n — fc — c) x 
{n ~ k — c) diagonal matrix whose entries are polynomials. The Smith form of i?2,3 {D) is full 
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rank because the original quantum check matrix S [D] is fuU rank. Perform the row operations 
in [D) and the column operations in {D). The quantum check matrix becomes 



E[{D) 



Ti{D) 

T{D) 



(7. 



where 



E[ {D) 



I 

A-^'{D) 

{D) = E2,2 {D)B-'iD) = 



El {D) , 



(7.9) 
(7.10) 



Perform a modified version of the Smith algorithm to reduce the (c — s) x (n — c) matrix E'^ i^) 
to a lower triangular form (discussed in Chapter 5). This modified algorithm uses only column 
operations to transform 



E2,2biD)^ L{D) 



(7.11) 



where L (D) is (c — s) x (c — s) and the null matrix is (c — s) x (n + s — 2c). The quantum check 
matrix becomes 



E[{D) 



Fi (£») 

r2{D) e'.^^ad) l{d) 

T{D) 



(7.12) 



We have now completed the decomposition of the quantum check matrix with column and row 

operations. 

We turn to showing how to encode and decode a certain quantum check matrix that proves 
to be useful in encoding the above quantum check matrix. Consider the following quantum check 
matrix 

7 
r2{D) L{D) 



(7.13) 
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where r2 (D) and L (D) are from the matrix in (7.12) and each of them, the identity matrix, and 
the null matrices have dimension (c — s) x (c — s). We use a method for encoding and decoding 
the quantum check matrix in (7.13) similar to the method in Chapter 5 for the second class of 
CSS entanglement-assisted quantum convolutional codes. We begin with a set of c — s ebits and 
c — s information qubits. The following matrix stabilizes the ebits 



(7.14) 



where Bob possesses the c — s qubits on the "left," Alice possesses the 2 (c — s) qubits on the 
"right," and each matrix is (c — s) x (c — s). The following matrix represents the information 
qubits: 



/ 


/ 























/ 


/ 












I 


























I 



(7.15) 



The above information-qubit matrix represents the logical operators for the information qubits 
and gives a useful way of tracking the information qubits while processing them. Tracking the 
information-qubit matrix helps to confirm that the information qubits decode properly at the 
receiver's end [52]. We track both the above stabilizer and the information-qubit matrix as they 
progress through some encoding operations. Alice then performs CNOT gates from her first c — s 
qubits to her next c — s qubits. These gates multiply the middle c — s columns of the "X" matrix 
by L (D) and add the result to the last c — s columns and multiply the last c — s columns of the 
"Z" matrix by and add the result to the last c — s columns. The stabilizer becomes 



I I 


and the information-qubit matrix becomes 





1 I L{D) 



(7.16) 



L^(i?-i) / 





0/ 



(7.17) 
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Alice performs infinite-depth operations on her first c — s qubits corresponding to the rational 
polynomials ■ • -i l2c-s (^~^) ™ -^^2^^ {^^^) ■ These operations multiply the middle 

c — s columns of the "Z" matrix by Fj^^ (-^^^) '^^d multiply the middle c — s columns of the "X" 
matrix by r2 {D). The stabilizer matrix becomes 



/ T^'{D~^) 




and the information-qubit matrix becomes 




/ T2{D) L{D) 



(7.18) 



L'^ {D-')T^^ {D--") I 





0/ 



(7.19) 



Alice's part of the above stabilizer matrix is equivalent to the quantum check matrix in (7.13) by 
row operations (premultiplying the first set of rows by r2 (D~^).) 

We now illustrate a way to decode the encoded stabilizer in (7.18) and information-qubit 
matrix in (7.19) so that the information qubits appear at the output of the decoding circuit. Bob 
performs CNOT gates from the first set of qubits to the third set of qubits corresponding to the 
entries in L (D). The stabilizer becomes 



I r^'{D-^) 
000 



and the information-qubit matrix becomes 





/ r2(£') 








7 



(7.20) 



(7.21) 



Bob finishes decoding at this point because we can equivalently express the information-qubit 
matrix as follows 

n n r n n n 

(7.22) 









I 


























I 
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by multiplying the first c — s rows of the stabilizer by L?" (-D^^) and adding to the first c — s rows 
of the inforniation-qubit matrix. The information qubits are available at the receiving end of the 
channel because the above information-qubit matrix is equivalent to the original one in (7.15). 

We show how to encode the quantum check matrix in (7.12) using ebits, ancilla qubits, and 
information qubits. We employ the encoding technique for the submatrix listed above and use 
some other techniques as well. Suppose that we have the following matrix that stabilizes a set of 
c ebits per frame, n — k — c ancilla qubits per frame, and k information qubits per frame: 



/ / 
0/0/000 













/ / 

7 / 

/ 



(7.23) 



The first and third sets of rows have s rows and correspond to s ebits per frame, the second and 
fourth sets of rows have c — s rows and correspond to c — s ebits per frame, and the last set of 
n—k — c rows corresponds to n — fc — c ancilla qubits per frame. The above matrix has n + c 
columns on both the "Z" and "X" side so that the above matrix stabilizes k information qubits 
per frame. Bob possesses the first c qubits and Alice possesses the next n qubits. Alice performs 
the encoding operations in (7.14-7.19) to get the following stabilizer: 



/ 





/ 






































/ 





(/?-!) 





















































/ 





/ 






































/ 





r2 {D) 





L{D) 






































/ 









(7.24) 



We perform several row operations to get the quantum check matrix in (7.12). Premultiply the 
middle set of rows by Fi (/?). Premultiply the last set of rows by E2 2a i^) ''^'^'^ the result to 
the set of rows above the last set. Premultiply the last set of rows by F (£)). Finally, premultiply 
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the first two sets of rows by E'{ [D) = E[ [D) [l r2 (-D"^)] and add the result to the last three 
sets of rows. The quantum check matrix becomes 



/ 7 

7 Fa ^ (7)-i) 



E'{{D) E[{D) 









TiiD) TiiD) 

7 












^2{D) E,^,JD) L{D) 
T{D) 



(7.25) 

□ 

Alice's part of the above quantum check matrix and the last three sets of rows are equivalent 
to the quantum check matrix in (7.12). Alice then performs all finite-depth encoding operations 
(column operations) in (7.1-7.12) in reverse order to obtain the desired quantum check matrix in 
the statement of the theorem. Decoding consists of performing all the operations in (7.1-7.12) and 
then applying the decoding operations in (7.20-7.22). The entanglement-assisted rate of the above 
code is k/n because the code uses a noisy quantum communication channel n times per frame to 
send k information qubits per frame. 



7.2 Example 

We now present an example that begins with the same generators as those in the example from 
the previous chapter. We begin with the following two Pauli generators: 

(• • • \IIII\ZXZI\ZZIZ\IIII\ ■■■), {■■■ \IIII\XYXI\XXIX\IIII\ ■■■). 



We write the above two generators as a quantum check matrix: 



l + D DID 
10 



10 

l + D l + D 1 D 



(7.26) 
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We encode two information qubits per frame with the help of two ebits. The stabihzer matrix for 
the unencoded qubit stream is as follows: 



110 

10 10 











110 

10 10 



(7.27) 



Rows one and three correspond to one ebit and rows two and four correspond to the other. 
Multiply row one by D and add the result to row three, multiply row one by 1 + + D"^ and 
add the result to row four, and multiply row two by 1 + D~^ and add the result to row four. These 
row operations give the following equivalent stabilizer: 



1 

D 



1 

D 





1 




l + D-2 l + £)-i + £)2 l + £)-i + £)2 1 + D- 











110 

10 10 



(7.28) 



Figure 7.1 illustrates the operations that transform the unencoded stabilizer to the encoded one 
in an online encoding circuit. The final stabilizer is as follows 



1 

D 



10 

1 + D DID 







10 





1 







1 



1 1 + D 1 + D 1 D 



Compare Alice's Paulis in the last two rows of the above matrix to the quantum check matrix 
in (7.26). We have constructed a code with the same error-correcting properties because these 
two matrices are equivalent. The entanglement-assisted rate of the above code is 1/2 because it 
encodes two information qubits for every four uses of the noisy quantum channel. 
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Figure 7.1: (Color online) An online encoding circuit for an entanglement-assisted quantum con- 
volutional code. The receiver Bob possesses the first two qubits in the two ebits and the sender 
Alice possesses the second two qubits in the two ebits. The sender encodes two information qubits 
per frame with the help of her half of the two ebits. 
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Consider the following two operators; 



D-^ 1 + D + D-^ 
00 l + £>2 00 



1 
1 



The first row anticommutes with the first row in (7.28) and commutes with all other rows in 
(7.28). The second row anticommutes with the second row in (7.28) and commutes with all 
other rows in (7.28). These commutation relations imply that the above operators are useful for 
encoding classical information in a superdense-coding-like fashion. These operators encode two 
classical bits into the code and make use of the first two rows in (7.28) instead of just 'Vasting" 
them. Measuring the first two rows in (7.28) reveals the values of the two classical bits. We can 
determine the encoded versions of these "classical-information-encoding" operators by tracing how 
the operators change in the Heisenberg picture through the rest of the encoding circuit. We can 
use these two classical bits and consume one ebit to tclcport an additional information qubit. This 
technique boosts the cntangiemcnt-assistcd rate of this code from 1/2 to 3/4 because wc can now 
encode three information qubits for every four uses of the noisy quantum communication channel. 



7.3 Closing Remarks 

The "free entanglement" approach of this chapter uses entanglement less efficiently than the 
protocol in the previous chapter. It does not require expanding a set of generators and therefore 
does not require a heuristic convergence argument as do the codes of the previous chapter. The 
free entanglement method always works and results in encoding and decoding circuits that act on 
smaller numbers of qubits than the circuits in the previous chapter. 
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Unified Quantum Convolutional Coding 



Chapter 8 



"Never chase after a bus, a girl, or a unifying theory. 
There'll always he another one coming along." 
—John Archibald Wheeler (1911-2008) 

In this chapter, wc design a framework for "grandfatlier'' quantum convolutional codes. Our 
grandfather codes are useful for the simultaneous transmission of classical and quantum informa- 
tion. Rather than using block codes for this purpose, we design quantum convolutional codes.* 
Our technique incorporates many of the known techniques for quantum coding: subsystem codes 
[27, 29], entanglement-assisted codes [34], convolutional codes [44, 45, 1], and classical-quantum 
mixed coding [63, 85, 84]. The goal of our technique is to provide a formalism for designing codes 
that approach the optimal triple trade-off rates in the grandfather resource inequality in (1.5). 

We structure this chapter as follows. Section 8.1 details our "grandfather" quantum convolu- 
tional codes. We explicitly show how to encode a stream of classical-quantum information using 
finite-depth operations and discuss the error-correcting properties of our codes. We end with an 
example of a grandfather quantum convolutional code. We discuss which errors the code corrects 
actively and others that it corrects passively. 

*Kremsky, Hsieh, and Brun address the formulation of grandfather block codes in Ref. [84]. 
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8.1 Grandfather Quantum Convolutional Codes 

We detail the stabilizer formalism for our grandfather quantum convolutional codes and describe 
how they operate. This formalism is a significant extension of the entanglement-assisted formalism. 

An [n,k,l;r,c] grandfather quantum convolutional code encodes k information qubits and I 
information classical bits with the help of c ebits, a = n— k — I — c— r ancilla qubits, and r gauge 
qubits. Each input frame includes the following: 

(i) . Alice's half of c ebits in the state ]$"'"). 

(ii) . a = n — k — c— I — r ancilla qubits in the state |0). 

(iii) . r gauge qubits (which can be in any arbitrary state a). 

(iv) . I classical information bits ■ ■ ■ x\ given by a computational basis state \x) = • • • (8) 

X^' |0)®'. 

(v) . k information qubits in a state |'^).'^ 

The left side of Figure ?? shows an example initial qubit stream before an encoding circuit 
operates on it. 

The stabilizer matrix {D) for the initial qubit stream is as follows: 



^0 {D) = 



/ / 

/ 




/ / 




(8.1) 



where all identity matrices in the first two sets of rows are c x c, the identity matrix in the last 
row is a X a, the three columns of all zeros in both the "Z" and "X" matrices are respectively 

(a + 2c) X r, (a + 2c) x I, and (a + 2c) x k. The first two sets of rows stabilize a set of c ebits and 
the last set of rows stabilize a set of a ancilla qubits. The first c columns of both the "Z" and "X" 
matrix correspond to halves of ebits that Bob possesses and the last n columns in both matrices 
correspond to the qubits that Alice possesses. 



tThis statement is not entirely true because tiie information qubits can be entangled across multiple frames, or 
with an external system, but we use it to illustrate the idea. 
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Different generators for the grandfather code are important in active error correction, in passive 
error correction, and for the identification of the I classical information bits. We first write the 
unencoded generators that act on the initial qubit stream. The first subgroup of generators is the 
entanglement subgroup Sea with the following generators: 



Se,o {D) = 



/ 





/ 



(8.2) 



The above generators are equivalent to the first two sets of rows in (8.1) acting on Alice's n qubits. 
The next subgroup is the isotropic subgroup <S/^o with the following generators: 



Sifi {D) 



0/000 00000 



(8.3) 



The above generators are equivalent to the last set of rows in (8.1) acting on Alice's n qubits. 
The encoded versions of both of the above two matrices are important in the active correction of 
errors. The next subgroup is the gauge subgroup Sg,o whose generators are as follows: 



Scfi {D) 



0/00 





0/00 



(8.4) 



The generators in <Sg,o correspond to quantum operations that have no effect on the encoded 
quantum information and therefore represent a set of errors to which the code is immune. The 
last subgroup is the classical subgroup <Sc,o with generators 



Scfi {D)= 000/0 00000. 



(8.5) 



The grandfather code passively corrects errors corresponding to the encoded version of the above 
generators because the initial qubit stream is immune to the action of operators in <Sc,o (up to a 
global phase). Alice could measure the generators in Scfl to determine the classical information in 
each frame. Unlike quantum information, it is possible to measure classical information without 
disturbing it. 
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Alice performs an encoding circuit with finite-depth operations to encode her stream of qubits 
before sending them over the noisy quantum channel. The encoding circuit transforms the initial 
stabilizer (D) to the encoded stabilizer S (D) as follows: 



S{D) 



I Zei{D) 
Ze2{D) 
Zj (D) 



Xei{D) 
Xi {D) 



(8.6) 



where Zei {D), Xei {D), Ze2 {D) and Xe2 {D) are c x n-dimensional, Zj (D) and Xj (D) are 
a X n-dimensional. The encoding circuit affects only the rightmost n entries in both the "Z" and 
"X" matrix of (D) because these are the qubits in Alice's possession. It transforms Se,o {D), 
Si^o (D), Sg,o (D), and ^co (D) as follows: 



Se {D) = 


Zei (D) 


Xei (D) 


, Si{D) = 




_ Ze2{D) 


Xe2 (D) _ 




Sg (D) = 


Zgi (D) 


Xgi (D) 


, Sc{D) = 




Zg2 {D) 


Xg2 {D) 





Zj {D) Xj {D) 
Zc {D) Xc {D) 



(8.7) 
(8.8) 



where Zgi {D), Xgi (-D), Zg2 (-D) and Xg2 {D) are r x n-dimensional and Zc {D) and Xc {D) 
are I x n-dimensional. The above polynomial matrices have the same commutation relations 

as their corresponding unencoded polynomial matrices in (8.2-8.5) and respectively generate the 
entanglement subgroup Se, the isotropic subgroup Si, the gauge subgroup Sg, and the classical 
subgroup Sc- 

A grandfather quantum convolutional code operates as follows. 

(i) . Alice begins with an initial qubit stream as above. She performs the finite-depth encoding 

operations corresponding to a specific grandfather quantum convolutional code. 

(ii) . She sends the encoded qubits online over the noisy quantum communication channel. The 

code passively protects against errors in {Si,Sg,Sc)- 

(iii) . Bob combines the received qubits with his half of the ebits in each frame. He obtains the 

error syndrome by measuring the generators in (8.6). He processes these syndrome bits with 
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a claswical error estimation algorithm to diagnose errors and applies recovery operations to 
reverse the errors. 

(iv). He then performs the inverse of the encoding circuit to recover the initial qubit stream 
with the information qubits and the classical information bits. He recovers the classical 
information bits either by measuring the generators in Sc before decoding or the generators 
in Sc,o after decoding. 

A grandfather quantum convolutional code corrects errors in a Pauli error set £ that obey one 
of the following conditions VEa, Eh G £: 

3ge {Si,Se) : {9,ElEb} = or eIe^ S {Si,Sg,Sc) • 

It corrects errors that anticommute with generators in {Si,Se) by employing a classical error 
estimation algorithm. The code passively protects against errors in the group (<S7,<Sg,<Sc). 

Our scheme for quantum convolutional coding incorporates many of the known techniques for 
quantum error correction. It can take full advantage of the benefits of these different techniques. 

8.2 Example 

We present an example of a grandfather quantum convolutional code in this section. The code 
protects one information qubit and classical bit with the help of an ebit, an ancilla qubit, and a 
gauge qubit. The first frame of input qubits has the state 

po = |$+) ($+1 |0) (0| ® (70 ® |a;o) (a;o| |Vo) (V'ol , (8.9) 

where l^"*") is the ebit, |0) is the ancilla qubit, ctq is an arbitrary state for the gauge qubit, |a;o) is a 
classical bit represented by state |0) or and |V'o) is one information qubit equal to ao \0)+(3o 
The states of the other input frames have a similar form though recall that information qubits 
can be entangled across multiple frames. 
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The initial stabilizer for the code is as follows: 



So (D) = 



1 1 

1 





1 1 




The first two rows stabilize the ebit shared between Alice and Bob. Bob possesses the half of the 
ebit in column one and Alice possesses the half of the ebit in column two in both the left and right 
matrix. The third row stabilizes the ancilla qubit. We name Alice's qubits one through five (they 
are actually two through six in the above matrix from the left to the right) . 

The generators for the initial entanglement subgroup Se,o, isotropic subgroup «S/,o, gauge 
subgroup Scfi, and classical subgroup Sc,o are respectively as follows: 



Sea (D) 
Si,o (D) 
Sg,o {D) 
Sc,o {D) 



1 



10 

10 







1 





10 



00010 00000 



The sender performs the following finite-depth operations (order is from left to right and top to 
bottom): 



H{2) C{2,3,D) C(2,4,l + £)) C(2,5,£>) H {3,A,5) 
C(2,3,D) C{2,5,D) H {2) C{l,2,D) C(l,4,l + D) 
C(l,5,l + L') F(l,2,3,4,5) C(1,3,L>) C(l,4,l + L») 
C(l,5,l + £') 5(1,4). 
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where the notation for the above encoding operations was estabhshed in Section 4.5. The initial 
stabiUzer matrix Sq (D) transforms to S (D) under these encoding operations, where 



SiD) 



1 

1 + D D 1 1 + D 
D D D 



1 + D D 1 1 + D 

1 
10 1 1 



B.IO) 



The generators for the different subgroups transform respectively as follows: 



Se (D) = 
Si{D) = 
Sg (D) = 
Sc{D) = 




1 + D D 1 1 + D 



1 + D D 1 1 + D 




D D D 



i) 1 i) 

0^000 



10 11 


10 



1 1 + D-^ 1 + D-i 0- 



The code actively protects against an arbitrary single-qubit error in every other frame. One can 
check that the syndromes of the stabilizer in S {D) satisfy this property. Consider the Pauli gen- 
erators corresponding to the generators in the entanglement subgroup and the isotropic subgroup: 



X I I X X 
Z I I Z Z 
I X I X X 



X I X I X 
Z Z I I Z 
I I Z Z Z 



(8.11) 



where all other entries in the left and right directions are tensor products of the identity. We 
can use a table-lookup syndrome-based algorithm to determine the error-correcting capability of 
the code. The method is similar to the technique originally outlined in detail in Ref. [1]. The 
syndrome vector s consists of six bits where s = Si - ■■ Se- The first bit Si is one if the error 
anticommutes with the operator XIIXX in the first part of the first generator above and zero 
otherwise. The second bit S2 is one if the error anticommutes with the operator XIXIX in the 
delayed part of the first generator above and zero otherwise. The third through sixth bits follow 
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Error 


Syndrome 


Error 


Syndrome 


Error 


Syndrome 




001100 




000001 


X5 


001101 


Yi 


111100 




010001 


n 


mill 


Zi 


110000 


Z3 


010000 


Z5 


110010 


X2 


000100 


X4 


001001 






Y2 


000110 


Y4 


101011 






Z2 


000010 


Z4 


100010 







Table 8.1: A list of possible single-qubit errors in a particular frame and the corresponding syn- 
drome vector. The syndrome corresponding to any single-qubit error is unique. The code therefore 
corrects an arbitrary single-qubit error in every other frame. 

a similar pattern for the second and third generators above. Table 8.1 lists all single-qubit errors 
over five qubits and their corresponding syndromes. The code corrects an arbitrary single-qubit 
error in every other frame using this algorithm because the syndromes are all unique. A syndrome- 
based Viterbi algorithm might achieve better performance than the simple syndrome table-lookup 
algorithm outlined above. 

This code also has passive protection against errors in (<S/,<Sg,«Sc)- The Pauli form of the 
errors in this group span over three frames and are as follows: 



mil 


IXIXX 


IIZZZ 


IZIZI 


IIZII 


mil 


IZIII 


IIXII 


mil 


IZIZI 


ZZIZI 


mil 



(8.12) 



The smallest weight errors in this group have weight two and three. The code passively corrects 
the above errors or any product of them or any five-qubit shift of them. 

There is a trade-off between passive error correction and the ability to encode quantum in- 
formation as discussed in Ref. [31]. One can encode more quantum information by dropping the 
gauge group and instead encoding extra information qubits. The gauge generators then become 
logical X and Z operators for the extra encoded qubits. One can also turn classical bits into 
information qubits by dropping the generators in the classical subgroup. These generators then 
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become logical Z operators for the extra encoded qubits. By making the above replacements, the 
code loses some of its ability to correct passively, but we gain a higher quantum information rate. 
On the other hand, if we replace a gauge qubit with an ancilla qubit or an ebit, we gain the ability 
to correct extra errors. The trade-off now is that replacing gauge qubits with ebits or ancillas 
enhances the active error-correcting capability of the code, but increases the overall complexity of 
error correction. 



8.3 Closing Remarks 

We have presented a framework and a representative example for grandfather quantum convolu- 
tional codes. We have explicitly shown how these codes operate, and how to encode and decode a 
classical-quantum information stream by using ebits, ancilla qubits, and gauge qubits for quantum 
redundancy. The ultimate goal for this theory is to find quantum convolutional codes that might 
play an integral part in larger quantum codes that approach the grandfather capacity [64]. One 
useful line of investigation may be to combine this theory with the recent quantum turbo-coding 
theory [58]. 
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Convolutional Entanglement Distillation 



Chapter 9 



A chap at the "Entanglement Distillery, " 
Was drunk so they gave him the pillory, 
They'd foul dirty ebits. 
He said, "Convolutional circuits!" 
So they augmented the quantum artillery. 

The goal of entanglement distillation resembles the goal of quantum error correction [65, 66]. An 
entanglement distillation protocol extracts noiseless, maximally-entangled ebits from a larger set 
of noisy ebits. A sender and receiver can use these noiseless ebits as a resource for several quantum 
communication protocols [13, 3]. 

Bennett et al. showed that a strong connection exists between quantum error-correcting codes 
and entanglement distillation and demonstrated a method for converting an arbitrary quantum 
error-correcting code into a one-way entanglement distillation protocol [66] . A one-way entangle- 
ment distillation protocol utilizes one-way classical communication between sender and receiver 
to carry out the distillation procedure. Shor and Preskill improved upon Bennett et al.'s method 
by avoiding the use of ancilla qubits and gave a simpler method for converting an arbitrary CSS 
quantum error-correcting code into an entanglement distillation protocol [86]. Nielsen and Chuang 
showed how to convert a stabilizer quantum error-correcting code into a stabilizer entanglement 
distillation protocol [32] . Luo and Devetak then incorporated shared entanglement to demonstrate 
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how to convert an cntanglcmcnt-assistcd stabilizer code into an entanglement-assisted entangle- 
ment distillation protocol [76]. All of the above constructions exploit the relationship between 
quantum error correction and entanglement distillation we further exploit the connection in this 
chapter by forming a convolutional entanglement distillation protocol. 

In this last chapter, our main contribution is a theory of convolutional entanglement distillation. 
Our theory allows us to import the entirety of classical convolutional coding theory for use in 
entanglement distillation. The task of finding a good convolutional entanglement distillation 
protocol now becomes the well-established task of finding a good classical convolutional code. 

We begin in Section 9.3 by showing how to construct a convolutional entanglement distillation 
protocol from an arbitrary quantum convolutional code. We translate earlier protocols [86, 32] for 
entanglement distillation of a block of noisy ebits to the convolutional setting. A convolutional 
entanglement distillation protocol has the benefit of distilling entanglement "online." This online 
property is useful because the sender and receiver can distill entanglement "on the fly" as they 
obtain more noisy ebits. This translation from a quantum convolutional code to an entanglement 
distillation protocol is useful because it paves the way for our major contribution. 

Our major advance is a method for constructing a convolutional entanglement distillation 
protocol when the sender and receiver initially share some noiseless ebits. As stated previously, 
prior quantum convolutional work requires the code to satisfy the restrictive self-orthogonality 
constraint, and authors performed specialized searches for classical convolutional codes that meet 
this constraint [43, 44, 48, 1]. We lift this constraint by allowing shared noiseless entanglement. 
The benefit of convolutional entanglement distillation with entanglement assistance is that we can 
import an arbitrary classical binary or quaternary convolutional code for use in a convolutional 
entanglement distillation protocol. The error-correcting properties for the convolutional entangle- 
ment distillation protocol follow directly from the properties of the imported classical code. Thus 
we can apply the decades of research on classical convolutional coding theory with many of the 
benefits of the convolutional structure carrying over to the quantum domain. 

We organize this chapter as follows. We review stabilizer entanglement distillation in Section 
9.1 and entanglement-assisted entanglement distillation in Section 9.2. In Section 9.3, we show how 
to convert an arbitrary quantum convolutional code into a convolutional entanglement distillation 
protocol. In Section 9.4, we provide several methods and examples for constructing convolutional 
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J 

Figure 9.1: An example of a block entanglement distillation protocol. A sender creates a set of 
noisy ebits by sending half of a set of Bell states through a noisy quantum channel. Both sender 
and receiver perform multi-qubit measurements to diagnose channel error. The sender transmits 
her measurement results to the receiver over a classical communication channel. Both perform 
recovery and decoding operations to obtain a set of noiseless ebits. 

entanglement distillation protocols where two parties possess a few initial noiseless ebits. These 
initial noiseless ebits act as a catalyst for the convolutional distillation protocol. The constructions 
in Section 9.4 make it possible to import an arbitrary classical binary or quaternary convolutional 
code for use in convolutional entanglement distillation. 

9.1 Stabilizer Entanglement Distillation without Entangle- 
ment Assistance 

The purpose of an [n, k] entanglement distillation protocol is to distill k pure ebits from n noisy 
ebits where < k < n [65, 66]. The yield of such a protocol is k/n. Two parties can then use 
the noiseless ebits for quantum communication protocols. Figure 9.1 illustrates the operation of 
a block entanglement distillation protocol. 

The two parties establish a set of shared noisy ebits in the following way. The sender Alice 
first prepares n Bell states |cf)+)'^" locally. She sends the second qubit of each pair over a noisy 
quantum channel to a receiver Bob. Let |<I>+) be the state 1$+)'^" rearranged so that all of Alice's 
qubits are on the left and all of Bob's qubits are on the right. The noisy channel applies a Pauli 
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error in the error set £ C H" to the set of n qubits sent over the channeh The sender and receiver 
then share a set of n noisy ebits of the form (I (g) A) |<i>+) where the identity I acts on Alice's 
qubits and A is some Pauli operator in £ acting on Bob's qubits. 

A one-way stabilizer entanglement distillation protocol uses a stabilizer code for the distillation 
procedure. Figure 9.1 highlights the main features of a stabilizer entanglement distillation protocol. 
Suppose the stabilizer <S for an [n, k] quantum error-correcting code has generators gi,.. ■ ,gn-k- 
The distillation procedure begins with Alice measuring the n — k generators in <S. Let {Pi} be 
the set of the 2""*' projectors that project onto the 2""'^ orthogonal subspaces corresponding to 
the generators in <S. The measurement projects |$+) randomly onto one of the i subspaces. Each 
Pj commutes with the noisy operator A on Bob's side so that 

(P, ® I) (I ® A) |$+) = (1 A) (P, I) |$+) . (9.1) 

The following important "Bell-state matrix identity" holds for an arbitrary matrix M: 

(M0I)|$+) = (I®M^)|$+). (9.2) 

Then (9.1) is equal to the following: 

(I ® A) (P, ® I) |$+) = (I ® A) (P2 ® I) |$+) = (I ® A) (P, ® Pf ) |$+) . 

Therefore each of Alice's projectors P^ projects Bob's qubits onto a subspace Pf corresponding 
to Alice's projected subspace Pj. This operation is one of the "weird" properties of entanglement 
because it projects the set of noisy ebits onto the codespace effectively "before" the noise acts on it. 
Alice restores her qubits to the simultaneous +l-eigenspa€e of the generators in <S. She sends her 
measurement results to Bob. Bob measures the generators in <S. Bob combines his measurements 
with Alice's to determine a syndrome for the error. He performs a recovery operation on his qubits 
to reverse the error. He restores his qubits to the simultaneous +l-eigenspace of the generators 
in <S. Alice and Bob both perform the decoding unitary corresponding to stabilizer <S to convert 
their k logical ebits to k physical ebits. 
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9.2 Stabilizer Entanglement Distillation with Entanglement 
Assistance 



Luo and Devetak provided a straightforward extension of the above protocol [76]. Their method 
converts an entanglement-assisted stabilizer code into an entanglement-assisted entanglement dis- 
tillation protocol. 

Luo and Devetak form an entanglement distillation protocol that has entanglement assistance 
from a few noiseless ebits. The crucial assumption for an entanglement-assisted entanglement 
distillation protocol is that Alice and Bob possess c noiseless ebits in addition to their n noisy 
ebits. The total state of the noisy and noiseless ebits is 

(I^® (A® 1)^)1$++,) (9.3) 

where I"* is the 2"+"^ x 2"+'= identity matrix acting on Alice's qubits and the noisy Pauli operator 
(A (g) I)^ affects Bob's first n qubits only. Thus the last c ebits are noiseless, and Alice and Bob 
have to correct for errors on the first n ebits only. 

The protocol proceeds exactly as outlined in the previous section. The only difference is that 
Alice and Bob measure the generators in an entanglement-assisted stabilizer code. Each generator 
spans over n + c qubits where the last c qubits are noiseless. 

We comment on the yield of this entanglement-assisted entanglement distillation protocol. 
An entanglement-assisted code has n — k generators that each have n + c Pauli entries. These 
parameters imply that the entanglement distillation protocol produces k+c ebits. But the protocol 
consumes c initial noiseless ebits as a catalyst for distillation. Therefore the yield of this protocol 
is k/n. 

In Section 9.4, we exploit this same idea of using a few noiseless ebits as a catalyst for distil- 
lation. The idea is similar in spirit to that developed in this section, but the mathematics and 
construction are different because we perform distillation in a convolutional manner. 
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Figure 9.2: An example of a convolutional entanglement distillation protocol taken from the 
quantum convolutional code in Ref. [1]. The code in Ref. [1] has rate 1/3 and can correct for 
single-qubit errors in every other frame. Alice and Bob first measure the operators in the stabilizer 
for the quantum convolutional code. Alice performs conditional unitaries on her qubits to restore 
them to the +1 eigenspace of the stabilizer code. Alice forwards her measurement results to Bob. 
Bob performs a maximum- likelihood decoding procedure such as Viterbi decoding [2] to determine 
the qubit errors. He corrects for these errors. He restores his qubits to the +1 eigenspace of the 
stabilizer code. Alice and Bob both perform online decoding to obtain ebits with yield 1/3. 

9.3 Convolutional Entanglement Distillation without En- 
tanglement Assistance 

We now show how to convert an arbitrary quantum convolutional code into a convolutional en- 
tanglement distillation protocol. Figure 9.2 illustrates an example of a yield-1/3 convolutional 
entanglement distillation protocol. The protocol has the same benefits as a quantum convolutional 
code: an online decoder with less decoding complexity than a block protocol, good error-correcting 
properties, and higher ebit yield than a block protocol. The protocol we develop in this section is 
useful for our major contribution presented in the next section. 

We can think of our protocol in two ways. Our protocol applies when a sender Alice and a 
receiver Bob possess a countably infinite number of noisy ebits. Our protocol also applies as an 
online protocol when Alice and Bob begin with a finite number of noisy ebits and establish more 
as time passes. The countably infinite and online protocols are equivalent. We would actually 
implement the entanglement distillation protocol in the online manner, but we formulate the 
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forthcoming mathematics with the countably infinite description. Each step in the protocol does 
not need to wait for the completion of its preceding step if Alice and Bob employ the protocol 
online. 

The protocol begins with Alice and Bob establishing a set of noisy ebits. Alice prepares a 
countably infinite number of Bell states |$+) locally. She sends one half of each Bell state through 
a noisy quantum channel. Alice and Bob then possess a state p"^^ that is a countably infinite 
number of noisy ebits pf^ where 

P^^ = (8)pf^. (9.4) 

i=l 

The state p"^^ is equivalent to the following ensemble 

{ft,i*+>r}- (9-5) 

In the above, pi is the probability that the state is |$+)^^, where 

|$+)f^^(I®A,)|$i)^^, (9.6) 

and I'J'Jj) is the state I |<1>+) J rearranged so that all of Alice's qubits are on the left and 
all of Bob's are on the right. € 11^^ is a Pauli sequence of errors acting on Bob's side. These 
errors result from the noisy quantum channel. I is a sequence of identity matrices acting on Alice's 
side indicating that the noisy channel does not affect her qubits. Alice and Bob need to correct 
for a particular error set in order to distill noiseless ebits. 

Alice and Bob employ the following strategy to distill noiseless ebits. Alice measures the 
n — k generators in the basic set Go- The measurement operation projects the first n{u+l) ebits 
{v is the constraint length) randomly onto one of 2""*^ orthogonal subspaces. Alice places the 
measurement outcomes in an (n — fc)-dimensional classical bit vector ao- She restores her half of 
the noisy ebits to the simultaneous +l-eigenspace of the generators in Qo if ao differs from the 
all-zero vector. She sends ao to Bob over a classical communication channel. Bob measures the 
generators in Qo and stores the measurement outcomes in a classical bit vector bo. Bob compares 
bo to ao by calculating an error vector eo = ao®bo. He corrects for any errors that eo can identify. 
He may have to wait to receive later error vectors before determining the full error syndrome. He 
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restores his half of the noisy ebits to the simultaneous +l-eigenspace of the generators in Qq if 
the bit vector bo indicates that his logical ebits are not in the +l-space. Alice and Bob repeat 
the above procedure for all shifts D (Go), (^o)i • • • of the basic generators in Q(f. Bob obtains 
a set £ of classical error vectors e^: £ = {e^ : i e Z^'}. Bob uses a maximum- likelihood decoding 
technique such as Viterbi decoding [2] or a table-lookup on the error set £ to determine which 
errors occur. This error determination process is a purely classical computation. He reverses the 
estimated errors after determining the syndrome. 

The states that Alice and Bob possess after the above procedure are encoded logical ebits. 
They can extract physical ebits from these logical ebits by each performing the online decoding 
circuit for the code Q. The algorithm outlined in Ref. [45] gives a method for determining the 
online decoding circuit. 

Example 9.3.1. We use the rate-1/3 quantum convolutional code in Example 4.1.1 to produce 
a yield- 1/3 convolutional entanglement distillation protocol. Alice measures the generators in 
the stabilizer in (4.8) for every noisy ebit she shares with Bob. Alice communicates the result 
of her measurement of the first two generators to Bob. Alice restores the qubits on her side to 
be in the simultaneous -t-l-eigenspace of the first two generators. Bob measures the same first 
two generators. Alice measures the next two generators, communicates her results, etc. Bob 
compares his results to Alice's to determine the error bit vectors. Bob performs Viterbi decoding 
on the measurement results and corrects for errors. He rotates his states to the simultaneous 
+l-eigenspace of the generators. Alice and Bob perform the above procedure in an online manner 
according to Figure 9.2. Alice and Bob can decode the first six qubits after measuring the second 
two generators. They can decode because there is no overlap between the first two generators 
and any two generators after the second two generators. They use the circuit from [45] in reverse 
order to decode physical ebits from logical ebits. They distill ebits with yield 1/3 by using this 
convolutional entanglement distillation protocol. The ebit yield of 1/3 follows directly from the 
code rate of 1/3. 
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9.4 Convolutional Entanglement Distillation with Entan- 
glement Assistance 

The convolutional entanglement distillation protocol that we develop in this section operates iden- 
tically to the one developed in the previous section. The measurements, classical communication, 
and recovery and decoding operations proceed exactly as Figure 9.2 indicates. 

The difference between the protocol in this section and the previous one is that we now assume 
the sender and receiver share a few initial noiseless ebits. They use these initial ebits as a catalyst 
to get the protocol started. The sender and receiver require noiseless ebits for each round of the 
convolutional entanglement distillation protocol. They can use the noiseless ebits generated by 
earlier rounds for consumption in later rounds. It is possible to distill noiseless ebits in this way 
by catalyzing the process with a few noiseless ebits. The protocol we develop in this section is a 
more powerful generalization of the previous section's protocol. 

The construction in this section allows sender and receiver to use an arbitrary set of Paulis for 
the distillation protocol. The set does not necessarily have to be a commuting set of Paulis. 

The implication of the construction in this section is that we can import an arbitrary binary 
or quaternary classical convolutional code for use as a quantum convolutional code. We explicitly 
give some examples to highlight the technique for importing. The error-correcting properties and 
yield translate directly from the properties of the classical convolutional code. Thus the problem 
of finding a good convolutional entanglement distillation protocol reduces to that of finding a good 
classical convolutional code. 

9.4.1 Yield (n-l)/n Convolutional Entanglement Distillation 

We present our first method for constructing a convolutional entanglement distillation protocol 
that uses entanglement assistance. The shifted symplectic product from Section 4.3 is a crucial 
component of our formulation. 

Suppose Alice and Bob use one generator N (u (D)) for an entanglement distillation protocol 
where 

n{D) = [z{D)\^{D)] = 
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Zi{D) ••• Zn{D) \ Xi{D) ••• Xn{D) 



and where z\ (D), .... z„ (D), xi (D), . . . , x„ (D) are binary polynomials. We do not impose a 
commuting constraint on generator N(u (£>)). Alice and Bob choose generator N(u(£))) solely 
for its error-correcting capability. 

The shifted symplectic product helps to produce a commuting generator from a noncommuting 
one. The shifted symplectic product of u {D) is 



The coefficient (u u)q for zero shifts is equal to zero because every tensor product of Pauli 
operators commutes with itself: 



Recall that u{D) is self-time-reversal symmetric (4.32). We adopt the following notation for a 
polynomial that includes the positive-index or negative-index coefficients of the shifted symplectic 
product (u © u) {D): 




(9.7) 



(u0u)o = O. 



(9.8) 




(9.9) 



The following identity holds: 



(u©u) {D)+ = (u©u) (D-i) 



(9.10) 



Consider the following vector of polynomials: 



a(D) 



(u©u) {D) 




(9.11) 



Its relations under the shifted symplectic product are the same as u{D): 



(a © a) (D) = (u © u) {D)+ + (u © u) (D) 



(u©u) (D). 



158 





Zi 




X2 


Z2 


Y2 


1 





1 


1 





1 














1 


1 





1 


1 


1 





1 


1 





1 












Table 9.1: The convolutional entanglement distillation protocol for Example 9.4.1 corrects for a 
single-qubit error in every fourth frame. Here we list the syndromes corresponding to errors Xi, 
Yi, and Zi on the first qubit and to errors X2, Y2, and Z2 on the second qubit. The syndromes 
are unique so that the receiver can identify which error occurs. 

The vector a (D) provides a straightforward way to make N (u (D)) commute with all of its shifts. 
We augment u (D) with a (D). The augmented generator u' (D) is as follows: 



2{D) (u©u)(I?)+ I x(D) 1 



(9.12) 



The augmented generator u' (D) has vanishing symplectic product, (u' u') (D) = 0, because the 
shifted symplectic product of a (D) nulls the shifted symplectic product of u (D). The augmented 
generator N (u' (D)) commutes with itself for every shift and is therefore useful for convolutional 
entanglement distillation as outlined in Section 9.3. 

We can construct an entanglement distillation protocol using an augmented generator of this 
form. The first n Pauli entries for every frame of generator N (u' (D)) correct errors. Entry 
n + 1 for every frame of N (u' (D)) makes N (u' (D)) commute with every one of its shifts. The 
error-correcting properties of the code do not include errors on the last (extra) ebit of each frame; 
therefore, this ebit must be noiseless. It is necessary to catalyze the distillation procedure with 
nu noiseless ebits where n is the frame size and v is the constraint length. The distillation 
protocol requires this particular amount because it does not correct errors and generate noiseless 
ebits until it has finished processing the first basic set of generators and u — 1 of its shifts. 
Later frames can use the noiseless ebits generated from previous frames. Therefore these initial 
noiseless ebits are negligible when calculating the yield. This construction allows us to exploit the 
error-correcting properties of an arbitrary set of Pauli matrices for a convolutional entanglement 
distillation protocol. 
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Wc discuss the yield of such a protocol in more detail. Our construction employs one generator 
with n + 1 qubits per frame. The protocol generates n noiseless ebits for every frame. But it also 
consumes a noiseless ebit for every frame. Every frame thus produces a net of n — 1 noiseless ebits, 
and the yield of the protocol is (n — 1) /n. 

This yield of (n — 1) /n is superior to the yield of an entanglement distillation protocol taken 
from the quantum convolutional codes of Forney et al. [1]. Our construction should also give 
entanglement distillation protocols with superior error-correcting properties because we have no 
self-orthogonality constraint on the Paulis in the stabilizer. 

It is possible to construct an online decoding circuit for the generator u' {D) by the methods 
given in [45]. A circuit satisfies the noncatastrophic property if the polynomial entries of all of 
the code generators have a greatest common divisor that is a power of the delay operator D [45] . 
The online decoding circuit for this construction obeys the noncatastrophicity property because 
the augmented generator u' (£>) contains 1 as one of its entries. 

Example 9.4.1. Suppose we have the following generator 



N (u (£))) = (• • • \II\ZZ\IX\XZ\ZI\II\ ■■■), 



where 



up) 



D 



The above generator corrects for an arbitrary single-qubit error in every fourth frame. Table 9.4.1 
lists the unique syndromes for errors in a single frame. The generator anticommutes with a shift of 
itself by one or two to the left or right. The shifted symplectic product confirms these commutation 
relations: 

(u Qu){D)^D + D'^ + D-^ + D-^. 

Let us follow the prescription in (9.12) for augmenting generator N(u(Z))). The following poly- 
nomial 



a(£>)= (u0u)(£')+ 



D + D^ 
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Figure 9.3: The above figure illustrates Bob's side of the convolutional entanglement distillation 
protocol that uses entanglement assistance. The noise affects the first and second of every three 
ebits that Bob shares with Alice. Every third ebit that Alice and Bob share are noiseless. The 
measurements correspond to those in Example 9.4.1. 
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has the same commutation relations as u (D): 



(a©a) (D) = (u0u) (D) . 



(9.13) 



We augment u (D) as follows: 



The overall generator now looks as follows in the Pauli representation: 



N (u' (£))) = (• • • \III\ZZX\IXZ\XZZ\ZII\III\ ■■■). 



The yield of a protocol using the above construction is 1/2. Figure 9.3 illustrates Bob's side of 
the protocol. It shows which of Bob's half of the ebits are noisy and noiseless, and it gives the 
measurements that Bob performs. 

9.4.2 Yield (n-m)/n Convolutional Entanglement Distillation 

The construction in the above section uses only one generator for distillation. We generalize the 
above construction to a code with an arbitrary number of generators. We give an example that 
illustrates how to convert an arbitrary classical quaternary convolutional code into a convolutional 
entanglement distillation protocol. 

Suppose we have the following m generators 



{N {ui (£))) : 1 < i < m} , 



where 



ui (D) 




Zip) 


XI (D) 


U2 (D) 




Z2 (D) 


X2 (D) 


Um {D)_ 






x„ (D) _ 



(9.14) 
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Wc make no assumption about the commutation relations of the above generators. We choose 
them solely for their error-correcting properties. 

We again utilize the shifted symplectic product to design a cxmvolutional entanglement distil- 
lation protocol with multiple generators. Let us adopt the following shorthand for the auto and 
cross shifted symplectic products of generators Ui {D) , Um {D): 



u+ = {ui © Ui) {D)+ , 
Uij = {ui Uj) {D) . 



Consider the following matrix: 



ai (D) 




a2 {D) 




a™ {D)_ 





Ui Ui,2 

u+ 



Ul,„ 
U2,r, 



(9.15) 
(9.16) 



(9.17) 



••• u+ 

The symplectic relations of the entries {D) are the same as the original Uj (D) : 

(ai © a^) (£>) = (ui © Uj) {D) Vz, j €{!,..., m}, 

or equivalently, if Cl^ (D) is the shifted symplectic product matrix for the generators in (9.17) and 
(£)) is the shifted symplectic product matrix for the "u" generators, then 



(D) = (D) . 

We mention that the following matrix also has the same symplectic relations: 











U2,l 



Um,l Um,2 



(9.18) 
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Let us rewrite (9.17) as follows: 



(9.19) 

The above matrix provides a straightforward way to make the original generators commute with all 
of their shifts. We augment the generators in (9.14) by the generators ci^ (D) to get the following 
m X 2 (n + m) matrix: 

zi(£)) z[{D) XI (D) x'i(D) 
Z2(£>) Z^(£>) X2(D) x'2(D) 

(D) {D) {D) ^'^ {D) _ 

Every row of the augmented matrix U' {D) has vanishing symplectic product with itself and any 
other row. This condition is equivalent to the following matrix condition for shifted symplectic 
orthogonality [44]: 

Z (D) (D-^) - X (£)) Z'^ (D-^) = 0. (9.20) 

The construction gives a commuting set of generators for arbitrary shifts and thus forms a valid 
stabilizer. 

We can readily develop a convolutional entanglement distillation protocol using the above 
formulation. The generators in the augmented matrix U' {D) correct for errors on the first n 
ebits. The last m ebits are noiseless ebits that help to obtain a commuting stabilizer. It is 
necessary to catalyze the distillation protocol with {n + m)u noiseless ebits. Later frames can 
use the noiseless ebits generated from previous frames. These initial noiseless ebits are negligible 
when calculating the yield. 

We comment more on the yield of the protocol. The protocol requires a set of m generators 
with n + m, Pauli entries. It generates n ebits for every frame. But it consumes m noiseless ebits 
per frame. The net yield of a protocol using the above construction is thus (n — m) / n. 



ai {D) 




z[{D) 


x'l {D) 


a2 {D) 




4(^) 




am {D) 









Z{D) X{D) 
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The key benefit of the above construction is that we can use an arbitrary set of Paulis for 
distilling noiseless ebits. This arbitrariness in the Paulis implies that we can import an arbi- 
trary classical convolutional binary or quaternary code for use in a convolutional entanglement 
distillation protocol. 

It is again straightforward to develop a noncatastrophic decoding circuit using previous tech- 
niques [45] . Every augmented generator in U' {D) has "1" as an entry so that it satisfies the 
property required for noncatastrophicity. 

Example 9.4.2. We begin with a classical quaternary convolutional code with entries from GF (4): 



(• • • |0000|la)10|1101|0000| • • • ) • (9.21) 

The above code is a convolutional version of the classical quaternary block code from Ref. [34]. 
We multiply the above generator by co and oj as prescribed in Refs. [19, 1] and use the map in 
(3.38) to obtain the following Pauli generators 



N (ui (£))) = (• • • \IIII\ZXZI\ZZIZ\IIII\ ■■■), 

N (U2 (£>)) = (• • • \IIII\XYXI\XXIX\IIII\ ■■■). (9.22) 



We determine binary polynomials corresponding to the above Pauli generators: 



" ui (D) ' 




l + DID 


10 


_ U2 (£>) _ 




10 


1 + D 1 + D 1 D 



The first generator anticommutes with itself shifted by one to the left or right, the second generator 
anticommutes with itself shifted by one to the left or right, and the first generator anticommutes 
with the second shifted by one to the left. The following shifted symplectic products confirm the 
above commutation relations: 



(ui ui) (£)) = £>-i + £>, (u2 0U2)(r>) = + {m Q U2) (D) = D-\ (9.24) 
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Consider the following two generators; 



" ai (D) 




D 


1 


_ a2 (D) _ 




D D 


1 



(9.25) 

Their relations under the shifted symplectic product are the same as those in (9.24). 

(ai ai) (£>) = (ui ui) (D) , (aa as) (£>) = (u2 U2) (£>) , (ai as) (£>) = (ui U2) (£») . 

We augment the generators ui (£>) and U2 (-D) to generators u'^ {D) and u'2 {D) respectively as 
follows. The augmented matrix U' {D) is 



\J'{D) = 



1 + D D I D D 
1 D D 



1 10 

1 + D 1 + D 1 D 1 



(9.26) 



The first row of U' (D) is generator u'l (D) and the second row is U2 (-D). The augmented gener- 
ators have the following Pauli representation. 



N (u'l {D)) = {■■■ \IIIIII\ZXZIXI\ZZIZZI\IIIIII\ ■■■), 
N (U2 (£>)) = (• • • \IIIIII\XYXIIX\XXIXZZ\IIIIII\ ■■■). 



(9.27) 



The original block code from Ref. [34] corrects for an arbitrary single-qubit error. The above 
entanglement distillation protocol corrects for a single-qubit error in eight qubits — two frames. 
This error-correcting capability follows from the capability of the block code. The yield of a 
protocol using the above stabilizer is again 1/2. 



9.4.3 CSS-Like Construction for Convolutional Entanglement Distilla- 
tion 

Wc finally present a construction that allows us to import two arbitrary binary classical codes for 
use in a convolutional entanglement distillation protocol. The construction is similar to a CSS 
code because one code corrects for bit flips and the other corrects for phase flips. 
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Wc could simply use the technique from the previous section to construct a convolutional 
entanglement-distillation protocol. We could represent both classical codes as codes over GF (4). 
We could multiply the bit-flip code by uj and the phase-flip code by Cu and use the map in (3.38) 
from GF (4) to the Paulis. We could then use the above method for augmentation and obtain a 
valid quantum code for entanglement distillation. But there is a better method that exploits the 
structure of a CSS code to minimize the number of initial catalytic noiseless ebits. 

Our algorithm below uses a Gram-Schmidt like orthogonalization procedure to minimize the 
number of initial noiseless ebits. The procedure is similar to the algorithm in [35] with some key 
differences. 

Suppose we have m generators {N (wj (£>)) : 1 < i < m} where 



(9.28) 



and each vector Wj {D) has length 2n. The above matrix could come from two binary classical 
codes. The vectors zi (£*),. . . ,Zp (D) could come from one code, and the vectors xi {D),. . . ,x„_p (D) 
could come from another code. The following orthogonality relations hold for the above vectors: 

V 1 < i, j < P : (wi w,) (£)) = 0, (9.29) 
V p+l<i',j' <m: (wj/ wy) {D) = 0. (9.30) 

We exploit the above orthogonality relations in the algorithm below. 

We can perform a Gram-Schmidt process on the above set of vectors. This process orthogo- 
nalizes the vectors with respect to the shifted symplectic product. The procedure does not change 

the error-correcting properties of the original codes because all operations are linear. 

The algorithm breaks the set of vectors above into pairs. Each pair consists of two vectors 
which are symplectically nonorthogonal to each other, but which are symplectically orthogonal 



wi {D) 




" zi (D) 





Wp{D) 









Wp+i {D) 







XI (£>) 


_ w„ (D) _ 







^m—p (-^) 
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to all other pairs. Any remaining vectors that are syniplectically orthogonal to all other vectors 
are collected into a separate set, which we call the set of isotropic vectors. This idea is similar to 
the decomposition of a vector space into an isotropic and symplectic part. We cannot label the 
decomposition as such because the shifted symplectic product is not a true symplectic product. 

We detail the initialization of the algorithm. Set parameters i = 0, c=0, 1 = 0. The index i 
labels the total number of vectors processed, c gives the number of pairs, and I labels the number 
of vectors with no partner. Initialize sets U and V to be null: U = V = 9. U keeps track of the 
pairs and V keeps track of the vectors with no partner. 

The algorithm proceeds as follows. While i < m, let j > 2c + / + 2 be the smallest index 
for a Wj (£>) for which (w2c+i+i Wj) {D) ^ 0. Increment I and i by one, add i to V, and 
proceed to the next round if no such pair exists. Otherwise, swap Wj (£>) with W2c+i+2 {D). For 
r e {2c + Z + 3, . . . , m}, perform 

Wr- {D) = (W2c+;+2 W2c+;+l) {D) Wr {D) + {Wr W2c+!+2) (-D"^) W2c+;+l (£>) , 

if Wr {D) has a purely z component. Perform 

Wr (-D) = (W2c+;+l W2c+;+2) {D) Wr {D) + {Wr W2c+;+l) (-D"^) W2c+(+2 {D) , 

if Wr (-D) has a purely x component. Divide every element in (-D) by the greatest common 
factor if the GCF is not equal to one. Then 

[Wr W2c+(+l) (£>) = [Wr W2c+;+2) {D) = 0. (9.31) 

Increment c by one, increment i by one, add i to U, and increment i by one. Proceed to the next 
round. 

We now give the method for augmenting the above generators so that they form a commuting 
stabilizer. At the end of the algorithm, the sets U and V have the following sizes: \U\ = c and 
|V| = I. Let us relabel the vectors Wj {D) for all 1 < i < 2c+ Z. We relabel all pairs: call the first 
Uj {D) and call its partner Vj {D) for all 1 < i < c. Call any vector without a partner Uc+i {D) for 
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all 1 < i < I. The relabeled vectors have the following shifted symplectic product relations after 
the Gram-Schmidt procedure: 



(u, V,) (D) = /, (D) V I, J e {!,..., c}, 
{uiQUj){D)=0 V i,j 

{viQvj){D)=0 V i,j G{l,...,c}, 



(9.32) 



where /» {D) is an arbitrary polynomial. Let us arrange the above generators in a matrix as 
follows: 



uJiD) vJ{D) 



u 



c+l 



(9.33) 



We augment the above generators with the following matrix so that all vectors are orthogonal to 
each other: 

... 









f2{D) 









Ocxl 
Olxl 



Ocxl 
Olxl 



Olxc 
Olxc 



Olxc 
Icxc 



(9.34) 



fc{D) 

Ocxl 
■■■ Olxl 

The yield of a protocol using the above construction is (n — m) /n. Suppose we use an [n, ki] 
classical binary convolutional code for the bit flips and an [n, classical binary convolutional 
code for the phase flips. Then the convolutional entanglement distillation protocol has yield 
(fci + k2-n) /n. 

Example 9.4.3. Consider a binary classical convolutional code with the following parity check 
matrix: 

1 + D D ll . (9.35) 
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We can use the above parity check matrix to correct both bit and phase flip errors in an entan- 
glement distillation protocol. Our initial quantum paxity check matrix is 



1 + D D 1 





1 + D D 1 



(9.36) 



The shifted symplectic product for the first and second row is D ^ + D. We therefore augment 
the above matrix as follows: 



1 + D D 1 D-'^ + D 





1 + D D 1 1 



(9.37) 



The above matrix gives a valid stabilizer for use in an entanglement distillation protocol. The 
yield of a protocol using the above stabilizer is 1/3. 



9.5 Closing Remarks 

We constructed a theory of convolutional entanglement distillation. The entanglement-assisted 
protocol assumes that the sender and receiver have some noiseless ebits to use as a catalyst for 
distilling more ebits. These protocols have the benefit of lifting the self-orthogonality constraint. 
Thus we are able to import an arbitrary classical convolutional code for use in a convolutional 
entanglement distillation protocol. The error-correcting properties and rate of the classical code 
translate to the quantum case. Brun, Devetak, and Hsieh first constructed the method for im- 
porting an arbitrary classical block code in their work on entanglement-assisted codes [35, 34]. 
Our theory of convolutional entanglement distillation paves the way for exploring protocols that 
approach the optimal distillable entanglement by using the well-established theory of classical 
convolutional coding. 

Convolutional entanglement distillation protocols also hold some key advantages over block 
entanglement distillation protocols. They have a higher yield of ebits, lower decoding complexity, 
and are an online protocol that a sender and receiver can employ as they acquire more noisy ebits. 

We suggest that convolutional entanglement distillation protocols may bear some advantages 
for distillation of a secret key because of the strong connection between distillation and privacy 
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[86] . We are currently investigating whether convolutional entanglement distillation protocols can 
improve the secret key rate for quantum key distribution. 
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Chapter 10 



Conclusion 



Fly Photon! Trap Ion! Stay Spin! 

Wliicli of you we'll implement in? 

For quantum coherence 
Demands perseverance, 
We'll try and God only knows when. 

Quantum error correction theory plays a fundamental role in quantum computing and com- 
munication. Without error-correcting protocols, quantum computing and communication devices 
will fall prey to the hands of decoherence. It is crucial for theorists to continue developing tech- 
niques to protect quantum information because a fundamental discovery in the theory might bring 
quantum computing closer to reality. 

We have augmented the theory of quantum error correction by contributing a theory of 
entanglement-assisted quantum convolutional coding. With the ability to import arbitrary clas- 
sical convolutional codes, we can now construct quantum convolutional codes that inherit the de- 
sirable characteristics of their ancestral classical convolutional codes. Our entanglement-assisted 
quantum convolutional coding theory should be useful for future quantum communication engi- 
neers if they would like to have codes with a good performance/complexity trade-off. 

We have said and again stress that the next important line of investigation is to combine the 
quantum turbo coding theory [58] with this theory. Convolutional codes form the constituent codes 
of a quantum turbo code and it would be interesting to investigate if we can enhance performance 
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with cntanglcmcnt-assistcd codes as the constituent codes. This investigation might produce 
quantum codes that come close to achieving the entanglement-assisted or "father" capacity. 

There are also other avenues to pursue these avenues include any scenario where entanglement 
assistance might help. It is useful to inspect the results of quantum Shannon theory to determine 
whether we can construct a coding scenario that fits with the constructions in the asymptotic 
scenario. We have conducted little analysis of the performance of our codes beyond stating that 
they inherit the properties of the imported classical codes. It would be interesting to observe 
the performance of these codes in a realistic noisy quantum channel when using syndrome-based 
Viterbi processing for correction of quantum error. 

We have seen much creativity in the field of quantum error correction in the past decade because 
of the many strange resources available in quantum theory. Experimentalists are increasingly using 
an array of quantum error correction techniques with the goal of bringing us closer to having qubits 
with good quantum coherence. One can only imagine what resources future quantum coding 
theorists will exploit to protect their valuable quantum information. 
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